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Chapter 1 
Introduction 



The quest for the fundamental building blocks and laws of the world surrounding 
us has been a driving force to mankind since its early days. The idea that nature 
consists of small, invisible constituents was first expressed by the ancient Greek 
Democritus in the fifth century BC. It was not until the nineteenth century AD 
that this idea was picked up again and embedded in a scientific context. Over 
time experiments discovered ever smaller substructures, from atoms to electrons 
and hadrons, and thereon to quarks. From a theoretical point of view, the aim 
is to embed these experimental results in a model which is based on as few 
assumptions as possible and can explain all other physical effects. 

The currently established model which performs this task is the Standard 
Model of elementary particle physics [1, 2]. It is one of the best-tested theories 
of contemporary physics. All known elementary particles are accommodated in 
this model. Solely the scalar Higgs boson [3] is included in the theory, but could 
not be found in experiments so far [4]. It is this particle which is assumed to be 
responsible for the masses of the fermions and weak gauge bosons. 

In spite of its success, the Standard Model also has its insufficiencies, and 
new theories are searched for, which might provide an even better description 
of nature. One of the most popular ones is supersymmetry [5]. It extends the 
two, fundamental symmetries of the Standard Model, the Poincare group and 
the non-Abehan gauge group SU{3)c ® SU{2)l ® U{1)y of strong, weak and 
electromagnetic interactions, by an anticommuting operator which induces an 
equal number of bosonic and fermionic states. 

The search for supersymmetry and the Higgs boson are main tasks of the Large 
Hadron Collider (LHC) at CERN. It will start operation in mid-2007 and provide 
a wealth of data. To verify or falsify theories and to relate this data to parameters 
of a model, it is necessary to calculate precise theoretical predictions, which match 
the accuracy which LHC will be able to obtain. As both the Standard Model and 
its supersymmetric extension are defined as perturbative theories with a series 
expansion in Planck's constant h, the inclusion of effects beyond leading order is 
often necessary. 
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In this thesis production processes for Higgs bosons in the Standard Model 
and its supersymmetric extension, the Minimal Supersymmetric Standard Model 
(MSSM) [6], at hadron colliders are considered. The calculations are performed 
at the one-loop level and include the SUSY-QCD corrections, i.e. corrections with 
squarks and gluinos running in the loop, for the MSSM Higgs bosons. 

The outline of this thesis is as following. First, a short introduction to the 
Standard Model (SM) is given in chapter 2. Special emphasis is put on the Higgs 
sector of the SM. Here also a possible extension including higher-order operators 
is discussed. Despite being a well-tested theory, the Standard Model also has its 
shortcomings, which are mentioned in the last section of this chapter. 

Out of the possible extensions of the Standard Model which aim to solve these 
deficencies, supersymmetry is the most popular one, as it is appealing from both 
an experimental and a theoretical point of view. Its discussion in chapter 3 of 
this dissertation starts with the basic principles of the theory. After the necessary 
ingredients to build a phenomenologically viable model are investigated, the focus 
is put on the simplest supersymmetric extension of the Standard Model, the 
Minimal Supersymmetric Standard Model (MSSM) [6]. The Lagrangian of the 
MSSM after supersymmetry breaking is written down and the particle content of 
the model is explained. 

Chapter 4 is concerned with the methods of regularization and renormaliza- 
tion. The first one is necessary to cancel the divergences which appear in the 
calculation of one-loop cross sections, and renders the amplitudes finite. Renor- 
malization then restores the physical meaning of the calculated cross sections. 
After a general introduction to the concepts, the renormalization of the strong 
coupling constant a^ in the way it is used in this thesis, is presented. The chap- 
ter concludes with a discussion of the bottom-quark Yukawa coupling. Here the 
mass counter term introduces large one- loop corrections to the cross section [7, 8]. 
They are universal, so they can be included in an effective tree-level coupling. 
Additionally, they are a one-loop exact quantity, so a resummation to all orders 
in perturbation theory is possible. 

The next chapter deals with the calculation of hadronic cross sections. The 
underlying theory, QCD and the parton model, is briefly introduced. Then ex- 
plicit formulae for the calculation of integrated and differential hadronic cross 
sections are given. An important technique to improve the cross-section ratio 
of signal over background processes and to enable the reconstruction of particu- 
lar event-types in the detector is the application of cuts to final-state particles. 
The implementation of these formulae in computer code is done in a program, 
called HadCalc, which is developed by the author of this thesis and which is 
lastly presented. It is based on the tools FeynArts [9, 10], FormCalc [11, 12, 13] 
and LoopTools [11, 14, 15]. The latter is extended to include now the five-point 
loop integrals, such that a complete one-loop calculation of 2 — > 3 processes is 
possible. HadCalc completes the tool set to provide a largely automated way of 
calculating hadronic cross sections. 



In the subsequent chapters, this program is apphed to the calculation of pro- 
cesses which contain supersymmetric Higgs bosons in the final state. The full 
one-loop SUSY-QCD corrections, i.e. corrections with squarks and gluinos run- 
ning in the loop, are included in the numerical results. 

The associated production of a charged Higgs boson H^ and a W boson via 
bottom quark-ant i- quark annihilation is studied in chapter 6. The discovery of a 
charged Higgs boson would be a clear sign of physics beyond the Standard Model. 
The above-mentioned universal corrections to the bottom-quark Yukawa coupling 
are expected to yield a numerically large and dominant contribution for certain 
regions of the MSSM parameter space, but the size of the SUSY-QCD corrections 
in the other regions is not known and requires a full one-loop calculation, which 
is presented in this thesis. 

In chapter 7 the production of the lighter CP-even neutral Higgs boson h^ 
via vector-boson fusion is investigated. This process has a clear final state of two 
jets in the forward region of the detector and forms an important /i^-pro duct ion 
mode with small theoretical uncertainties. For the corresponding Standard Model 
process with a Standard Model Higgs boson H in the final state, the Standard- 
QCD corrections are already known. They are the same as for /i^-production 
in the MSSM up to the replacement of the Higgs coupling. In the MSSM case 
additional SUSY-QCD corrections appear. In this thesis the complete one-loop 
SUSY-QCD corrections are calculated and their effect on the total cross section 
is discussed. In the last section of this chapter a background to the vector-boson- 
fusion process, /i°-pro duct ion with two outgoing jets and one or two gluons in 
the initial state, is considered and its numerical impact studied. 

The SUSY-QCD corrections to /i°-pro duct ion in association with heavy, i.e. 
bottom or top, quarks are presented in chapter 8. Besides being additional dis- 
covery channels for the Higgs boson, these processes can also be used to extract 
the respective quark Yukawa couplings from the data. This task can only be per- 
formed if the theoretical uncertainty of the cross section is small. The Standard- 
QCD corrections to these processes are available in the literature and greatly 
reduce the dependence on the renormalization and factorization scale. Addi- 
tionally, there are SUSY-QCD corrections which can also yield large corrections 
and must be taken into account. Therefore, a full calculation of the one- loop 
SUSY-QCD corrections is necessary, which is presented in this dissertation. 

Lastly, the possibility to measure the quartic Higgs coupling at hadron collid- 
ers is analyzed in chapter 9. For this purpose triple-Higgs production via gluon 
fusion is studied at the leading one-loop order. In this chapter not the MSSM 
is used as the underlying model, but an effective theory based on the Standard 
Model where the trilinear and quartic Higgs self-couplings are left as free param- 
eters. 

In appendix A the numerical values of the Standard Model parameters, which 
were kept fixed for all calculations in this thesis, and of the MSSM parameters 
for the reference point SPSla' [16] are noted. Appendix B contains the defini- 
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tions of mathematical quantities which are used throughout the dissertation, and 
appendix C the parametrization of the phase space for two- and three-particle 
final states. 

In loop calculations integrals over the loop momentum appear which can be 
solved analytically. The definition of these integrals is given in appendix D. Spe- 
cial attention is paid to the five-point integrals which have not been implemented 
in the package Loop Tools [11, 14, 15] before. The numerical method of Gaussian 
elimination, which is used to further improve the stability of the loop-integral 
calculation, is presented in appendix E. 

Finally, the complete user manual of HadCalc is attached in appendix F. The 
program itself can be obtained from the author^. 



"'^email: mraucMmppmu . mpg . de 



Chapter 2 



Standard Model 



2.1 Structure of the Standard Model 

The Standard Model (SM) of elementary particle physics [1, 2] is one of the best 
tested theories in physics. It consists of an outer symmetry of the Poincare group 
of space-time transformations and a non-Abelian gauge group of the inner direct 
product SU (3) c® SU (2) l(E)U {1)y ■ SU{3)c is the color gauge group and describes 
the strong interactions by the theory of QCD. The product SU{2)l^U{1)y spec- 
ifies the electroweak interactions which unify the electromagnetic and weak in- 
teractions. The Higgs mechanism, which will be described in chapter 2.2, breaks 
this symmetry spontaneously, thereby leaving a U{1)q symmetry of electromag- 
netic interactions which is described by QED. The one remaining interaction, 
gravitational interaction, is beyond the scope of the SM. In fact, a consistent 
theory which formulates general relativity in terms of a quantum field theory 
is not known until today. At the center-of-mass energies used at present or at 
planned future colliders, which are maximally of the order of a few hundred TeV, 
the effects due to gravitational interactions are negligibly small. The Standard 
Model therefore provides an excellent approximation to describe collider physics. 
The fermionic sector of the SM consists of spin-| leptons (z/e,t'^,z/T-,e,/i,r) and 
quarks {u,c,t,d,s,b) which appear in three different generations. The particles of 
each generation have the same quantum numbers but a different coupling to the 
Higgs field which will be introduced below. Left-handed fermions transform as a 
doublet under SU{2)l where the upper component forms the neutrinos (i^e,t'/i,'^r) 
and up-type quarks {u,c,t), respectively, and the lower component the electron- 
type leptons {e,fi,T) and the down-type quarks {d,s,b). Right-handed fermions 
transform as a singlet under SU{2)l the only exception being that there are no 
right-handed neutrinos at all. For each group generator a spin-1 gauge boson 
exists which transforms under the adjoint representation of the respective group. 
Consequently there are eight gauge bosons for SU{3)c, the gluons, three gauge 
bosons for SU(2)l, the W bosons, and one for U{1)y, called B. 
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Experiments show that not all gauge bosons are massless [17]. Adding an ex- 
plicit mass term for these gauge bosons is not possible for renormalizable quantum 
field theories. Such terms are forbidden due to the postulate that the Lagrangian 
should be invariant under gauge transformations. Otherwise the resulting theory 
would be non- renormalizable. For this reason another way of giving masses to 
the gauge bosons is needed. This is achieved by the Higgs mechanism which will 
be described in the next chapter. 



2.2 Higgs mechanism 

2.2.1 Standard Model Higgs sector 

As mentioned above, it is a difficult task to construct a gauge theory which is 
renormalizable and has massive gauge bosons. In the Standard Model this prob- 
lem is solved by the Higgs mechanism [3]. The idea is to add additional terms to 
the Lagrangian, such that the Lagrangian is invariant under the SU{2)i^U{1)y 
gauge transformations with a ground state which does not share this invariance. 
To realize this idea one introduces a new complex scalar field, the Higgs field 
$, which behaves like a doublet under SU{2)i gauge transformations and has 
hypercharge Y = +\. Its ground state acquires a vacuum expectation value w, 
such that a U{1)q symmetry of electromagnetic interactions is preserved. The 
electromagnetic charge is defined as Q = /s + y , where /s is the quantum number 
of the third component of the weak isospin operator. Therefore only the lower 
component of the doublet can have a vacuum expectation value, as assigning a 
vacuum expectation value to the upper component would also break the U{1)q. 
The Higgs field can be parametrized as 

where G^ is a complex and H and G" are two real scalar fields. The Higgs 
potential, i.e. the non-kinematic part of the SM Lagrangian which contains only 
Higgs fields, can be written as 

y($)=-!|i($t$) + ||($t$)2_ (2.2) 

The breaking of a continuous global symmetry leads to massless scalar par- 
ticles, the Goldstone bosons [18]. One Goldstone boson occurs for each broken 
generator of the symmetry group. In case of a broken continuous local symme- 
try, like a gauge symmetry, these Goldstone bosons are unphysical. They can be 
eliminated by an appropriate choice of gauge, the unitary gauge. Their degrees of 
freedom are "eaten up" by the gauge bosons which become massive. Once "eaten 
up", the Goldstone bosons form the longitudinal modes of the gauge bosons. 
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For electroweak symmetry breaking there are three broken generators leading to 
three "would-be" Goldstone bosons G^ and G^. Only the field H in eq. (2.1) is 
physical. It is the field of the Higgs boson which has not been discovered yet. Its 
mass m,H is a free parameter of the theory. It is bounded from below by exper- 
imental searches m^ > 114.4 GeV [19] and from above by electroweak precision 
data where a best fit yields mn = 114^45 GeV [20]. 

After electroweak symmetry breaking the gauge boson triplet ly* z = 1 . . . 3, 
of SU{2)l and the gauge boson B^ (f/(l)y) no longer form the mass eigenstates 
of the theory. The mass eigenstates are obtained by rotations 

Wi~[WlTiWl), Z,=cwWl-swB,, A,=swWl + cwB,. (2.3) 

sw and cw denote the sine and cosine of the electroweak mixing angle, the Wein- 
berg angle. The photon field A^ stays massless and can be interpreted as the 
gauge boson of the remaining U{1)q symmetry of electromagnetic interactions. 
The electromagnetically neutral Z and the charged W bosons receive a mass, 
which is proportional to the vacuum expectation value of the Higgs field: 

6 6 

mz =- V, mw =- — V (2.4) 

ZSwCw ZSw 

where e is the electromagnetic unit charge. As W and Z have already been found 
in experimental searches these equations determine the Weinberg angle and the 
scale of electroweak symmetry breaking v = 247 GeV. 

The Goldstone bosons G^ and G^ of eq. (2.1) are absorbed by the W and 
Z bosons, respectively. In this thesis the 't Hooft-Feynman gauge is used which 
has technical advantages for loop calculations since the gauge boson propagators 
in this gauge take a simpler form. In the 't Hooft-Feynman gauge the Goldstone 
bosons appear explicitly as internal propagators with a mass equal to that of the 
associated gauge boson. For external propagators their contribution is accounted 
for in the longitudinal component of the polarization vector of the respective 
gauge boson. 

In analogy to the inclusion of massive gauge bosons into a renormalizable 
quantum field theory there is no possibility to introduce fermion mass terms 
directly. To generate fermion masses one introduces Yukawa interactions which 
couple the fermions to the Higgs field 

>CYukawa = " XI'Lj^cr^j - A^^Q.^^M^^j - X'/Qj^dn,j + h.c. (2.5) 

with 

$^ = ^a2^* = (_^J+*) (2.6) 

which is also an SU{2)l doublet but has hypercharge Y = —1. The vacuum 
expectation value v in the decomposition of $ (eq. (2.1)) leads to terms which 
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are bilinear in the fermion fields, i.e. to mass terms for the fermions. The Xy 
are 3x3 Yukawa coupling matrices. They parameterize the masses of the quarks 
and further mixing effects in the quark sector. 

2.2.2 Higher-dimensional operators 

The realization of the Higgs sector in the SM is minimal in the sense that it 
contains just enough additional parameters and fields to give a consistent the- 
ory of the particles known nowadays. In extensions of the SM additional terms 
are possible, which lead to the following general parameterization of the Higgs 
potential with one doublet $ [21, 22]: 

^(*) = E:^(***-f)"" = ^»(***-f)' + '^(ij)- '") 

The expansion for n = on the right-hand side is identical to the SM Higgs 
potential eq. (2.2) with Aq = 2:J- up to the constant term which is not a physical 
observable and leaves the equations of motion unchanged. The additional terms 
for n > contain operators of mass dimension 6 and higher. Such terms are non- 
renormalizable but can be considered as effective terms of an extended theory. 
They are suppressed by the scale A which is the scale where new physics sets 
in. The only requirement eq. (2.7) has to fulfill is that its highest no n- vanishing 
coefficient Aj is positive so that the potential is bounded from below. 



2.3 Problems of the Standard Model 

Despite its large success there are both experimental and theoretical hints that 
the SM is only the low-energy limit of a more general theory. 

An experimental clue is the measured value of the anomalous magnetic mo- 
ment of the muon [23] . This observable is known to an extremely high precision 
from both experiment and theory, where the uncertainty stems from unknown 
higher-loop contributions and experimental errors on the input parameters. The 
deviation from the SM prediction is about 0.7-3.26 standard deviations [24]. 

Another evidence comes from the dark matter problem in the universe [25]. 
Looking at the rotation of galaxies as a function of the distance from the cen- 
ter shows that for large distances the circular velocity is constant, whereas the 
observed radiating matter would result in a decrease of the velocity with the 
distance. This implies that there is some fraction of matter which is contribut- 
ing to the overall mass density of the galaxy, but not emitting electromagnetic 
radiation, hence the name dark matter. Precision measurements of the cosmic 
microwave background [26] yield an average density of the universe that is very 
close to the so-called critical density, where the curvature of the universe van- 
ishes. Combining these data with our current understanding how the universe 
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emerged and evolves requires that the total matter content of the universe which 
contributes to this density is about 27%. The rest is some form of energy, so- 
called dark energy. Of these 27% of matter content, only about 4% of the total 
matter content consist of the usual baryonic matter, i.e. of matter built up of 
protons and neutrons. The remaining 23% must be made of non-baryonic, only 
weakly-interacting matter. The only particles in the SM which fulfill this require- 
ment are the neutrinos. Current upper limits on their masses [17] imply however 
that they cannot account for the whole required dark matter density. 

One of the theoretical clues is the unification of coupling constants in Grand 
Unified Theories (GUT), where all three SM gauge groups merge in a single 
gauge group. Possible GUT gauge groups are SU{5) [27], which is experimentally 
not viable due to a too large proton decay rate [28] or 5*0(10) [29]. Via the 
renormalization group equations the coupling constants of the three SM gauge 
groups can be written as running coupling constants which depend on the energy. 
GUT theories predict that at a high energy scale, typically of the order Mqut ~ 
10^^ GeV, all three gauge couplings unify. Such a unification does not occur in 
the SM, even if one takes into account that new particles at the GUT scale might 
slightly modify the running. 

Another hint is the so-called hierarchy problem. If one considers one-loop cor- 
rections to the mass of the Higgs boson quadratic divergences appear [30]. These 
divergences can be erased by renormalization. One finds that the corrections are 
of the order of the largest mass in the loop. If the SM is indeed the ultimate the- 
ory up to arbitrary high energies, this heaviest particle is the top quark and the 
corrections are well under control. But if the SM is replaced by a new theory at 
higher energies, like a Grand Unified Theory which unifies the electroweak with 
the strong interactions or a quantum theory which includes gravity, new particles 
with masses of the order of this new theory will appear, typically with masses of 
the Planck scale Mpianck ~ 10^^ GeV. In such new models extreme fine-tuning is 
necessary to get a Higgs mass of the order of the electroweak scale, as is predicted 
by electroweak precision data [20]. In particular there is no symmetry, neither 
conserved nor broken, which would explain such a fine-tuning in a natural way. 

The last problem concerns the neutrino sector. Neutrinos are assumed to be 
massless in the SM. It is known from the observation of neutrino oscillations [31] 
that neutrinos possess a tiny mass. There is no conceptual problem to introduce 
such a mass in the SM. As neutrinos are not of importance for the work presented 
in this dissertation the exact formulation of the neutrino sector can be ignored. 

To solve the problems mentioned above various models have been proposed. 
The model widely believed to be the most promising candidate is supersymmetry. 
This extension of the Standard Model was studied in this thesis and will be 
introduced in the following. 
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Supersymmetry 



3.1 Basic principles 

It was shown by Coleman and Mandula [32] that combining space-time and in- 
ternal symmetries is only possible in a trivial way. In the proof of this theorem 
only general assumptions on the analyticity of scattering amplitudes and the as- 
sumption that the S-matrix is invariant under Lorentz transformations are made. 
Later it was realized [33] that besides of Lie-algebras, which are defined via 
commutation relations, one can also use so-called superalgebras, which also con- 
tain anticommutators. Then a new type of operators Q is allowed which has the 
following properties [6, 34, 35]: 

\Qa ,Ql3 } = {QaA,Qi3B} =0 

[^M'Qal = [Pf^^QaA] =0 (3.1) 

The supersymmetry generators Q and Q carry Weyl spinor indices a, a, (3 and 
/3 which run from 1 to 2, where the undotted indices transform under the (0, |) 
representation of the Poincare group and the dotted ones under the (|,0) con- 
jugated representation. The indices A and B refer to an internal space and run 
from 1 to a number A^ > 1. For A^ > 1 chiral fermions are not allowed [36]. 
These are necessary to construct the observed parity violation via SU{2)l, where 
left- and right-handed fermions carry different quantum numbers. Therefore only 
(A^ = l)-supersymmetries are relevant for phenomenologically interesting energy 
ranges and in the following only such supersymmetries will be considered. P^ 
denotes the generator of Lorentz translations, the energy-momentum operator, 
and a'^: = (1, a* l) is the four-dimensional generalisation of the Pauli matrices. 
The first line of eq. (3.1) shows the entanglement of space-time symmetry and 
the internal symmetry. The last line indicates the invariance of supersymmetry 
under Lorentz transformations. 

11 
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As the operators anticommute with themselves, they have half-integer spin 
according to the spin-statistics theorem. A detailed calculation shows that their 
spin is always |. Therefore we have 

Q |boson) = |fermion) Q |fermion) = |boson) . (3-2) 

The one-particle states belong to irreducible representations of the supersym- 
metry algebra, the so-called supermultiplets. Each supermultiplet includes both 
bosonic and fermionic states which are called superpartners to each other. They 
can be transformed into each other by applying Q and Q. 

Each supermultiplet contains the same number of bosonic and fermionic de- 
grees of freedom. For example the simplest supermultiplet incorporates a Weyl 
fermion with two helicity states, hence two degrees of freedom. Its bosonic part- 
ners are two real scalars each with one degree of freedom, which can also be 
combined into one complex scalar field. This is called the scalar or chiral super- 
multiplet. 

The next possibility is a spin-1 vector boson. To guarantee the renormalis- 
ability of the theory this has to be a gauge boson which is massless and contains 
two degrees of freedom. It follows that the partner is a massless Weyl fermion. 
A spin- 1 fermion would render the theory non-renormalisable, so it must be a 
spin-| fermion. This is called a gauge or vector supermultiplet. 

From eq. (3.1) follows 

[P^P^Q,] = [P^P^g<i] =0 . (3.3) 

Pfj^P^ is just M^, the squared mass of a state in the supermultiplet. Applying 
the supersymmetry operator therefore does not change the mass of the state and 
all states in a supermultiplet have the same mass if supersymmetry is unbroken. 
This will be important later on when the Lagrangian is constructed. 



3.2 Superfields 

Starting from the supermultiplets one can construct superfields. To simplify the 
notation Grassmann variables are introduced. These are anticommuting numbers 
whose properties are defined in chapter B.3. The superalgebra can now be written 
in terms of commutators 



0"Qa,QpO^ =2r(T^e^P^ (3.4) 

[e''Qo.,9fQp] = \Qj'',Q^9'']=0 . (3.5) 



In general a finite supersymmetric transformation is given by the group ele- 
ment 

G {x„ r , e^) = eH-^^''+''"Q"+^"<^<3n , (3.6) 
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in complete analogy to a general non-Abelian gauge transformation e**^""^" with 
the group generators T*^. P'^, Q^ and Q" are the generators of the supersym- 
metry group. The coordinates can be combined into a tuple which represents a 
superspace coordinate z = {x^,6°',9a). The set of all possible coordinates spans 
the superspace. 

The fields on which these generators operate must then also be a function of 
9 and 6 besides x^. These are the so-called superfields $(a;^, 6, 6). In superspace 
one can obtain an explicit representation of Qa and Q" as differential operators. 
For that purpose one considers a supersymmetry transformation of $ 

G{y,,i,l)^{x,e,e). (3.7) 

Taking the parameters as infinitesimal and performing a Taylor expansion one 
obtains the following explicit representation of the supersymmetric generators 

<?"--J. + '«^-S«^ (3.9) 



d_ 



P. = ^l^- (3-10) 



For the further treatment it is sufficient to consider only infinitesimal super- 
symmetric transformations which have the following form 



Sg{^,O^x,,0, 



4+4-^(w^"-^-.a ^ 



36 ^de ^^ ^ '^^' dx^. 



$(x^,^,^), (3.11) 



where ^ and ^ are also Grassmann variables. Contracted indices which are 
summed over have been suppressed in this equation. 

Analogously to the covariant derivative in gauge theories one also introduces 
covariant derivatives D^ and D^ with respect to the supersymmetry generators. 
These derivatives have to be invariant under Q and Q, which is equivalent to the 
postulate 

{D,, Qc?i = {Da, Qa} = {Da, Qa} = {Da, Qa} = 0. (3.12) 

Thus the covariant derivatives are 

D^ = ^ + ia\e('^ (3.13) 

Da = -T^ - ^^^^Lt^- (3-14) 

QQa /3a ^3,/. V J 

From eqs. (3.8, 3.9, 3.13, 3.14) one can also deduce that the Grassmann variables 
6 and ^ have spin — |, while D and Q have spin +|. 
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Superfields can be expanded into component fields. The general expansion of 
superfields in terms of Grassmann variables is 

$(a;, e, 6) =f{x) + e(t){x) + Oxix) 

+ eem{x) + een{x) + ea^ev^{x) (3.15) 

+ eee\{x) + eeei){x) + eeeed{x). 

Due to the anticommuting properties of Grassmann variables this expansion is 
complete, i.e. it truncates with the last shown term. 

Up to now all expressions have been written out for general superfields. To 
construct a supersymmetric Lagrangian only two special types of superfields are 
needed. They are irreducible representations of the supersymmetry algebra. One 
obtains them by imposing covariant restrictions on a general superfield. In this 
way they still span a representation space of the algebra but have less components. 

3.2.1 Chiral Superfields 

One possibility are chiral superfields. They are defined by applying the covariant 
derivative Da on the scalar superfield $ as defined in eq. (3.15) 

Da^{z) = 0. (3.16) 

The solution of this differential equation leads to a chiral superfield which can be 
expressed in general component fields as 

$ =A{x) + iOa^'ed^A^x) + -66668 ^d^'A^x) 

(3.17) 
+ V26i){x) - -=66d^,i){x)a^'6 + 66F{x). 

y4 is a complex scalar field, (f) a complex Weyl spinor and F an auxiliary complex 
scalar field which has mass dimension two. It transforms under supersymmetry 
transformations into a total space-time derivative and therefore does not represent 
a physical, propagating degree of freedom. The product of chiral superfields is 
again a chiral superfield. For two chiral superfields $1 and $2 this follows directly 
from the product rule for derivatives 

Da ($1$2) = {Da^l) $2 + $1 (5^$2) = 0. (3.18) 

Analogously one can define an antichiral superfield \l/ by the equation 

Da^{z) = 0. (3.19) 

In particular the hermitian conjugate <l>^ of a chiral superfield $ is an antichiral 
superfield. 
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3.2.2 Vector Superfields 

The second special type of superfields are vector superfields. They are derived 
from a general scalar superfield V by demanding it to be real: 

V\z) = V{z). (3.20) 

The name vector superfield stems from the fact that in the expansion a real vector 
field appears as a component field and that these fields are used as generalized 
gauge fields when supersymmetric gauge theories are constructed. 
The complete expansion in terms of component fields is 

V{x, e, 6) =C{x) + i9x{x) - i6x{x) + Uo [M{x) + iN{x)] - ^66 [M{x) - zN{x)] 



ea^9vJx)+ie99 



^M 



A(x) + \cj>^d,x{x) 



ieee 



\{x) + '-a'^d.xix) 



-eeee 

2 



D{x) + \d^d^C{x) 



(3.21) 

C, D, M and A^ are scalar fields and v^ is the vector field which gives the name 
to this type of superfields. They all have to be real to fulfill eq. (3.20). A and x 
are Weyl spinors. 

For the vector superfield we can now define a supersymmetric gauge transfor- 
mation which is in the general non-Abelian case described by 

e9y -^ ^-ig^^gv^ig^ ^ (3 22) 

where $ denotes again a chiral superfield. This simplifies in the Abelian case to 

V ^V + i{^-<^^). (3.23) 

Using this gauge transformation we can simplify eq. (3.21) and choose 

X(x) = C{x) = M{x) = N{x) = (3.24) 

thereby eliminating unphysical degrees of freedom. This choice of gauge is called 
Wess-Zumino gauge [5]. As we have used only three of the four bosonic degrees 
of freedom in $ the "ordinary" gauge freedom of an Abelian gauge group is still 
present and the Wess-Zumino gauge is compatible with the usual gauges. 
The vector superfield is now simplified to 

V = -ea^ev^{x) + ieeex{x) - me\{x) + -eeeeD{x) (3.25) 

with the scalar auxiliary field D with mass dimension two. As in the case of chiral 
superfields this auxiliary field turns into a total derivative under supersymmetry 
transformations and does not contribute to the propagating degrees of freedom. 
Now we have all building blocks to construct a supersymmetric extension of 
the Standard Model. 
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3.3 A Supersymmetric Lagrangian 

A supersymmetric Lagrangian requires the action to remain unchanged under 
supersymmetry transformations 

Sg I d^xC{x) = . (3.26) 

This requirement is fulfilled if the Lagrangian C turns into a total space-time 
derivative under supersymmetry transformations. A comparison with the trans- 
formation properties of chiral and vector superfields shows that the F and D 
terms of eq. (3.17) and (3.21) show exactly this behavior. Schematically the 
Lagrangian can be written simply as 

£= [d^9CF+ [ d^ed^eCo . (3.27) 



As was noted already in the previous chapter the product of two chiral super- 
fields is again a chiral superfield. Explicit multiplication of the component fields 
yields a term proportional to ipiipj which has the form of a fermion mass term. The 
product of three chiral superfields which is by induction also a chiral superfield 
contains terms of the type ipiipjAk which describe Yukawa-like couplings between 
two fermions and a scalar. Products of four or more chiral superfields would lead 
to terms with a mass dimension greater than four and yield a Lagrangian which 
is no longer renormalizable. Thus the terms which can contribute to a supersym- 
metric Lagrangian can be written in a compact way with the superpotential 

iy($,) = A,$, + -rriy^.^j + -gijk^i^j^k ■ (3.28) 

The product $$''' of a chiral superfield with its hermitian conjugate is self- 
conjugate. Therefore it is a vector superfield according to the definition eq. (3.20) 
and a possible candidate for a contribution to Co'- 

I d^^d^*^"^ = FF* - Ad^d^'A^ - i^a^d^'i: . (3.29) 

The expression contains terms for the kinetic energy of both the scalar and the 
fermionic component. The auxiliary fields F do not have any kinematic terms so 
they can be integrated out. 

Gauge interactions are introduced by a supersymmetric generalization of the 
"minimal coupling" $"I'<1> -^ $te2(yy,|, ^j^]^ g^ vector superfield V with V = T^V^, 
where T" are the generators of the gauge group. Written in component fields one 
can replace the ordinary derivatives by covariant derivatives D^ = d^ + igv'^Ta. 

The terms for the kinetic energy of the gauge fields can also be expressed in 
terms of a superpotential 

W^ = -^{DD)e-^'^Dy3^ . (3.30) 
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The product WaW"" is gauge invariant and also a chiral superfield, so its ^^-term 
can appear in the supersymmetric Lagrangian. Again only the gauge bosons and 
their superpartners, the gauginos, obtain kinetic terms, but not the auxiliary 
fields, so we can eliminate them. 

Therefore the most general form of a supersymmetric Lagrangian has the 
following form: 



■^susY — I d 6 



^ -yV'^W'' + PF ($) ) + h.c. 



dHdH ($te23^$) 



16^2 

(3.31) 
As the two auxiliary fields F and D do not have any terms for the kinetic 
energy, their equations of motion have a simple form 

Solving these equations for the F and D fields 

'dW{A. 



dA, 



Da = -gAirjAj (3.33) 



and inserting these expressions into eq. (3.31) the Lagrangian can be completely 
expressed in terms of physical fields. 

3.4 Supersymmetry breaking 

As shown in eq. (3.3) all members of a supermultiplet must have the same mass. 
This means if the Standard Model was supersymmetrized by just replacing the 
fields with their respective superfields there would exist for example a supersym- 
metric partner to the electron with a mass of 511 keV/c^. This partner particle 
is a boson with spin 0, but otherwise with the same quantum numbers as the 
electron, i.e. particularly with a charge of one negative elementary charge. Such 
a particle would have been discovered experimentally a long time ago. 

This problem can be circumvented by requiring that supersymmetry is bro- 
ken. In this way one can give a mass to the supersymmetric partners which is 
beyond the current experimental limits. In analogy to spontaneous symmetry 
breaking in the electroweak sector the Lagrangian itself should be invariant un- 
der supersymmetry transformations, but have a vacuum expectation value which 
is not invariant under such transformations. For supersymmetry this problem is 
somewhat complicated because additional constraints appear which have to be 
fulfilled simultaneously. Such a constraint follows immediately from the definition 
of the supersymmetry algebra eq. (3.1) which implies 

H = P'> = ^ {QiQi + QiQi + Q2Q2 + Q2Q2) > 0. (3.34) 
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Applying the Hamiltonian H onto a state |\l/) leads to the result that supersym- 
metry is broken if neither the D nor the F term can be made zero simultaneously. 

The Fayet-Iliopoulos mechanism [37] achieves supersymmetry breaking by 
adding a D term to the Lagrangian which is linear in the auxiliary field, while 
O'Raifeartaigh models [38] do this via chiral supermultiplets and a superpotential 
such that not all auxiliary fields F can be made zero at the same time. Both mech- 
anisms are phenomenologically not viable because they can lead to color breaking 
or the breaking of electromagnetism, or need an unacceptable fine-tuning [39] . 

Hence one expects that supersymmetry is not broken directly by renormal- 
izable tree-level couplings, but indirectly or radiatively. For these purposes one 
introduces a hidden sector of particles in which supersymmetry is broken and 
which has only small or no direct couplings at all to the normal visible sector. 
The two sectors however share some common interaction which mediates the 
breaking from the hidden to the visible sector and leads to additional super- 
symmetry breaking terms. Two possible scenarios for this mediation are widely 
discussed in the literature [40] . The first one is gravity-mediated supersymmetry 
breaking. At the Planck scale gravity is anticipated to become comparable in 
size to the gauge interactions. The mediating interaction is associated with the 
new gravitational interactions which enter at this scale. Because of the fiavor 
blindness of gravity these gravitational interactions are expected to be fiavor- 
blind as well. A second possibility is that the mediating interactions are the 
ordinary QCD and electroweak gauge interactions. They connect the visible and 
the hidden sector via loop diagrams involving messenger particles. This scenario 
is called gauge-mediated supersymmetry breaking. 

For a phenomenological analysis it is often not relevant what the exact way 
of supersymmetry breaking is but only which additional terms in the Lagrangian 
are generated. Thereby the cancellation of quadratic divergences should remain 
valid, such that the solution of the naturalness problem of the Standard Model 
is not lost. Terms which do not spoil the cancellation are called soft supersym- 
metry breaking terms. It was shown [41] that only the following terms are soft 
supersymmetry breaking up to all orders in perturbation theory: 

• scalar mass terms m1jA*Aj 

• trilinear scalar interactions tijkAiAjAk + h.c. 

• mass terms for gauge particles -niiXiXi 

• bilinear terms bijAiAj + h.c. 

• linear terms hAi 

Now all building blocks are in place and we can turn to building a supersym- 
metric version of the Standard Model. 
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3.5 Minimal Supersymmetric Standard Model 

The simplest possibility of a supersymmetric extension of the Standard Model is 
called Minimal Supersymmetric Standard Model (MSSM). The underlying alge- 
bra is an (N=l)-supersymmetry with soft supersymmetry breaking. As in the 
Standard Model the MSSM shall have a local gauge symmetry with respect to 
the gauge group SU{3)c ® SU{2)l ® t/(l)y, which describe the strong, weak 
and electromagnetic interactions. Its particle content is obtained by replacing all 
fields with their corresponding superfields. 

Each matter field is assigned a chiral superfield. Its fermionic part describes 
the usual fermions of the Standard Model and its bosonic part contains the "scalar 
fermions" , the sfermions, as superpartners. For each gauge group a vector super- 
field is introduced whose vector bosons form the usual gauge bosons of the Stan- 
dard Model, and the fermionic superpartners are two-component Weyl spinors, 
in general called gauginos. The nomenclature of the new particles usually follows 
the convention that the bosonic superpartners carry the name of the fermion with 
a prefix "s", which is short for "scalar", and the fermionic superpartners carry 
the name of the gauge boson with the suffix "-ino" . 

In the Higgs sector of the MSSM it is not sufficient to replace the scalar field 
by a vector superfield. One would need both the field H and its hermitian con- 
jugate H* to give mass to both up- and down-type quarks. This is forbidden 
by the requirement that the superpotential must be analytic and so one needs a 
second Higgs doublet with negative hypercharge. Additionally the fermion which 
emerges from the single Higgs superfield would carry a non- vanishing hypercharge 
Y. This hypercharge contributes to the chiral anomaly [42] which is not com- 
pensated by other particles. The quantized version of such a theory would be 
inconsistent. In the MSSM the two fermions, one from each Higgs doublet, have 
opposite hypercharge and their contribution to the anomaly exactly cancels. 

Table (3.1) gives an overview of the particle content of the MSSM in the 
interaction basis. For the gauge superfields we have the following field strengths 
in the MSSM 

Wc"^ = -J (DD) e-^'^^D^e^3.G (3 35^, 

WLi = ~ {DD) e-23»^D«e23-^ (3.36) 

Wva = -\ (DD) e-'^'y^Dy^yB ^ _9^i)^jj^B . (3.37) 

Additionally the superpotential must be fixed. In the MSSM it is defined as 

(3.38) 
where Ae, A„ and A^ are 3x3 Yukawa couphng matrices and I and J denote the 
generation index. 
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Table 3.1: Superfields and particle content of the MSSM in the interaction basis. 
Superfields are denoted with a hat and the superpartners all carry a tilde. The 
generation index I of the quarks and leptons runs from 1 to 3. For the gauge 
fields the color index a runs from 1 to 8 and the weak isospin index i from 1 
to 3. The bold numbers in the group representation of the non-Abelian groups 
SU{'i)c and SU{2)l denote the dimension of the representation, where 1 is the 
trivial representation and the gauge bosons are in the adjoint representation of 
the group. The number for the Abelian group f/(l)y denotes the hypercharge of 
the particle. 
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Inserting eqs. (3.35)-(3.38) into eq. (3.31) and adding the F terms to the 
Lagrangian yields the supersymmetric part of the MSSM Lagrangian 

ajjr aa , -•- tjt i jjr ia 






-^susY = d 9 

+ T7r^WYcWY'" + VTmssm ) + h.c. 



^■■^"■■^ ™y 






lii^2g^W+2gyBl^^^^^2gyB^ ^ ' > 



j^ Q\ ^2gsG+2g^W+2gyB Q _^ fj] ^-2gsG'^ +2gyBfj _^ jj] ^-2gsG'^ +2gyB f^ 
j^fj\^2g^W+2gyBjj^ _^ jj] ^2g^W+2gyB ^^ 

Supersymmetry in the MSSM is broken exphcitly by soft supersymmetry 
breaking terms, i.e. only the terms mentioned at the end of chapter 3.4 are al- 
lowed. This leads to the following contributions to the MSSM Lagrangian: 

• Majorana mass terms for all gauginos 

Aoft.majoranamass =1^ (MsA" A^ + MsA^At^ + MiA^Ab) + H.C. (3.40) 

• mass terms for all scalar superpartners of the Standard Model fermions and 
for the scalar Higgs fields 

'^soft.scalarmass = — M^ \^L,I^L,I + ^L,I^L,l) — Me^*Rj(^R,I 

- M? (uIjUlj + d*LjdL,ij - M^u*jijURj - Mld*j^jdRj 
-ml\Hi\'^ -ml\H2\'^ (3.41) 

• bilinear term which couples the two scalar Higgs fields 

-Csoft.bilinear =''^12 y^ijH{Hl + h.C.) (3.42) 

• trilinear interaction terms for the scalar superpartners of the Standard 
Model fermions 



+ h.C. . (3.43) 
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In the general case the Yukawa couphngs Ae, A„ and A^^ as well as the trilinear 
couplings Af,, Au and A^ are complex 3x3 matrices. The scalar mass parameters 
Ml, M^, Mq, Mjj, Mfj, are hermitian 3x3 matrices. The scalar Higgs mass 
parameters mi and m2 are real numbers, and the gaugino mass parameters Mi, 
M2 and M3 as well as the bilinear Higgs coupling mu are complex numbers. 

There is an additional possibility for soft-breaking trilinear couplings [43] 
which has the form 

>Csoft,tri2 = {a!J'hI*V'EJ - A'J'hI*Q''UJ + A'J'W^*Q''DJ^ + h.c. . (3.44) 

This expression involves charge-conjugated Higgs fields which, in contrast to the 
superpotential, are possible for soft supersymmetry-breaking terms. However, it 
turns out that in most scenarios of supersymmetry breaking such terms are not 
generated. Therefore they are normally not considered and will also be neglected 
in this thesis. 

The complete soft supersymmetry breaking Lagrangian is given by 

''-soft ''-'soft,majoranamass ' •'-softjScalarmass ~r •'-soft, bilinear ' •'-soft, trilinear yO.^Oj 

As next step gauge fixing terms must be added to the Lagrangian. This is 
required so that all Green functions are still calculable. In this dissertation the 
R^- or 't Hooft gauge is used 

/:gauge-fixing = " ^ (d'^G;)' - ^ (^d'^W^ + -^m^^ (G+ - 0')^ 

- ^ (d'^W^, - ^rn^^ (G+ + G-))' (3-46) 

- ^ {d'K + ^wmziG'Y - ^ [d^B^ - swniziG^f . 

G^ and G° are the Goldstone bosons which were already described for the Stan- 
dard Model case in chapter 2.2 and appear in the MSSM in the same way. 

Setting ^ = 1 yields the 't Hooft-Feynman gauge which is advantageous for 
one-loop calculations, because the propagators take a very simple shape, while 
the Goldstone bosons appear explicitly in the calculation. This kind of gauge is 
used throughout this thesis. 

Finally unphysical modes which were introduced by the gauge-fixing terms 
are compensated by Faddeev- Popov ghost terms £ghost [44]. 

Adding up all contributions gives the complete Lagrangian of the MSSM 

-^MSSM = -^SUSY + -^soft + -^gauge-fixing + -^ghost- (3.47) 

Additional terms could be added to the superpotential in eq. (3.38) which are 
also gauge-invariant and analytic in the superfields, but violate lepton or baryon 
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number conservation which has not been observed experimentally so far. Such 
terms include the coupling of three lepton or quark superfields or the coupling 
of lepton to quark superfields. The strictest limits on lepton and baryon number 
violation are obtained by searching for a possible decay of the proton which vio- 
lates each baryon and lepton number by one unit. Experiments have established 
a lower limit on the proton lifetime of 10^^ years [17] while general violating 
terms predict a decay time in the order of minutes or hours. Thus a mechanism 
must exist which forbids or at least heavily suppresses these terms. The simplest 
possibility is to postulate a conservation of baryon and lepton number. Such a 
postulate would be a regression with respect to the Standard Model. There the 
conservation is fulfilled automatically and a consequence of the fact that there 
are no renormalizable lepton and baryon number violating terms. Furthermore, 
postulating lepton and baryon number conservation as a fundamental principle 
of nature is generally not viable. It is known that there are non-perturbative 
effects in the electroweak sector which do violate lepton and baryon number 
conservation, although their effect is negligible for the energy ranges of current 
experiments. 

Instead a symmetry should be introduced which has the conservation of these 
quantum numbers as a natural consequence. So in the MSSM as a perturbative 
theory baryon and lepton number conservation is again guarantueed while the 
existence of non-perturbative effects is not contradicted by demanding a funda- 
mental symmetry. Such a symmetry is given by R-parity [45]. A new quantum 
number R is introduced and from that a so-called i?-parity Pr = (—1)^ is de- 
rived. It is induced by the generators of supersymmetry, stays intact after spon- 
taneous supersymmetry breaking and is multiplicatively conserved. -R = for 
all Standard Model particles and the additional Higgs scalars and R = 1 for all 
supersymmetric partners. The link to lepton and baryon number conservation is 
obvious if one writes the -R-parity quantum number in terms of baryon number 
B, lepton number L and spin s 

P^ = (_1)2^+3{B-L) , (3 48) 

S is -|-| for the left-handed chiral quark superfield Q/, — | for the right-handed 
quark superfields Uj and Dj, and for all other particles. Analogously L is +1 
for the left-handed lepton superfield Lj, —1 for the right-handed lepton superfield 
Ej, and for all other particles. Then all Standard Model particles and the Higgs 
scalars have Pr = +1 and the supersymmetric partners have an odd -R-parity of 
Pr = -1. 

An interesting consequence of this is that each interaction vertex connects an 
even number of supersymmetric particles. Therefore they can only be produced 
in pairs and the lightest supersymmetric particle (LSP) must be stable. 
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3.6 Particle content of the MSSM 

3.6.1 Higgs and Gauge bosons 

As in the Standard Model the SU{2)l(E)U{1)y symmetry is broken by the vacuum 
expectation values of the Higgs fields in such a way that a U{1)q symmetry of 
electromagnetic interactions remains. Its associated conserved quantum number 
is the usual electromagnetic charge. As shown before the Higgs sector of the 
MSSM must consist of two scalar isospin doublets 

with opposite hypercharge. (pi, 02, Xi ^^^ Xo, ^-^e real scalar fields, and (p^ and 
(j)2 are complex scalar fields. In eq. (3.49) an expansion around the vacuum 
expectation values has been performed, which satisfy the equation 

{//,> = (;) (»,) = («) . (3.50) 

Collecting all terms in the Lagrangian which contain only the Higgs fields we 
have contributions to the Higgs potential from the F terms in the superpotential, 
from the D terms and finally from the soft supersymmetry breaking terms 

^Higgs = l/^l {\Hl\ + \H2\ ) 

+ 1 [al + 9i) im' - WY + Igl \HlH2f (3.51) 

+ ml \Hif + m\ \E2'\ - mj {eijH{Hi + h.c) . 

This equation shows the close entanglement between supersymmetry breaking 
and electroweak symmetry breaking. Only including the soft breaking terms it 
is possible that the minimum of the Higgs potential is not at the origin and the 
fields acquire a vacuum expectation value. 

The mass matrices of the Higgs fields are obtained by differentiating twice 
with respect to the fields and x- This leads to four uncoupled real 2x2 
matrices. To obtain the mass eigenstates these matrices have to be diagonalized 
by unitary matrices. In the case of a real 2x2 matrix this is simply a rotation 
matrix. We obtain 

I) . (3.54) 



A" 
H" 
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C/3, sp, Ca and Sa is a short-hand notation for cos/3, sin/3, cos a and sin a, respec- 
tively. Similar abbreviations will also be used for the other angles in this thesis, 
as well as tf3 denoting tan/3. The mixing angle /3 is defined as the ratio of the 
two vacuum expectation values 

tf3 = — with < /3 < ^ . (3.55) 

t/3 is a free parameter of the MSSM. The mixing angle a is determined by the 
relation 

t2a = hr^-^^ with - ^ < a < . (3.56) 

m\ — m| 2 

The restriction on the given interval determines a uniquely and is chosen such 
that always m/jo < tjiho. By electroweak symmetry breaking three group genera- 
tors are broken and therefore as in the Standard Model three unphysical would-be 
Goldstone bosons G^ and G° emerge. The five remaining Higgs bosons are phys- 
ical ones. There are two electrically neutral CP-even Higgs bosons h^ and if*^, 
one CP-odd A^ and two electrically charged ones H^. The mass of the CP-odd 
Higgs boson mA is usually chosen to be the second free parameter of the MSSM 
Higgs sector. The masses of the other Higgs bosons at tree-level are then 

mhO,HO =2 ( "^A + "il =F Y ("^i + 'm\f - Am\m%clp \ (3.57) 

mH± =m\ + m^ . (3.58) 

These relations receive large corrections at higher orders which must be taken into 
account when one wants to obtain realistic predictions. The one-loop corrections 
are known completely [46, 47, 48, 49]. On the two-loop level the calculation 
of the supposedly dominant corrections in the Feynman diagrammatic approach 
[50] of 0{atas) [51, 52, 53, 54, 55], O {aj) [51, 56, 57], C («,,«,) [58, 59] and 
O {atah + a1) [60], a calculation in the effective potential approach [61] and the 
evaluation of momentum-dependent effects [62] have been performed. As these 
expressions are rather lengthy they are not written out here. For the numerical 
evaluation the expressions given in [63] have been used. 

As in the Standard Model, electroweak symmetry breaking turns the W^ and 
B gauge bosons into the mass eigenstates W^, Z and the photon 7. W and Z 
bosons acquire a mass, where the single vacuum expectation value of eq. (2.1) is 
replaced by w = \/v\ + f|. 

The gauge bosons of SU{?>)c are the eight massless gluons. Their mass eigen- 
states are identical to the interaction eigenstates g't = G'^ . 

3.6.2 Higgsinos and Gauginos 

All particles which have the same quantum numbers can mix with each other. As 
the SU{2)l ® U{1)y symmetry is broken, only the SU{?>)c and U{1)q quantum 
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numbers have to match. 

In the sector of non-colored, charged particles there are the Winos W^ and 
the charged Higgsinos H^ and H2 with 

M 



r± 



-^>^w 



W^=\ ^ 



i>^w 



Ht = \i^. 



Ht 




(3.59) 



As for the W bosons the relation 
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holds. 

These four two-component Weyl spinors combine into two four- component 
Dirac fermions called charginos. Their mass matrix is diagonalized by 
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M2 V2mwSi3\ 
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U and V are two unitary matrices which are chosen such that m 
positive and m-+ < Tn-+. The chargino mass eigenstates are given by 
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(3.61) 
are both 

(3.62) 



The uncolored neutral higgsinos and gauginos also mix among each other. We 
have the two neutral Higgsinos }i\ and H\^ the Zino Z and the Photino A 

■|) ^^^(i ^m ^4t 



E\ 



(3.63) 



The latter two are obtained, as in the case of Z and 7, by rotating A^ and A^ 
by the Weinberg angle 

^Z =CvyAjy — S-w\b A^ =SvyA^ -|- CvyAs . (3.64) 

The four Weyl spinors form four Majorana fermions, called neutralinos, whose 
mass matrix is also diagonalized by a unitary matrix N 
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(3.65) 
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Again the remaining freedom in the choice of A^ is used to order the neutrahno 
masses such that m^a < m^o < m^o < 171.^0. The neutrahno mass eigenstates are 

Xi — X2 — X3 — X4 ° 

given by 



xl 
\xV 



N 



Hi 

V m ) 



(3.66) 



The gauginos of 5't/(3)c7, the gluinos, do not mix with other particles as they 
are the only fermions which are subject to the strong interaction exclusively. 
There are eight gluinos with mass rUg = \M^\. Gluinos are Major ana particles 
and have the following form 



r 



-lie 



(3.67) 



3.6.3 Leptons and Quarks 

Leptons and quarks have similar properties as in the Standard Model. The Weyl 
spinors of left- and right-handed fermions can be combined into one Dirac spinor 
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(3.68) 



where / is again the generation index. The down-type quarks dj are not exact 
mass eigenstates. A rotation 

d'l = VS^KMdj (3.69) 

by a unitary matrix, the Cabibbo-Kobayashi-Maskawa(CKM)-matrix Vckm [64], 
is required to turn the flavor eigenstates dj into mass eigenstates d'j. As the 
CKM-matrix is close to a unity matrix and flavor-mixing effects do not play any 
role in the processes which are calculated in this thesis effects induced by the 
CKM-matrix will be neglected and the CKM-matrix is set to exactly the unity 
matrix. 

Leptons and quarks receive their masses via the Yukawa terms in the super- 
potential which are bilinear in the lepton and quark flelds: 
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(3.70) 



These equations are often rewritten such that the Yukawa couplings are expressed 
in terms of the fermion masses and the mass of the W boson 
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e denoting the elementary charge. 
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(3.71) 
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3.6.4 Sleptons and Squarks 

In the sfermion sector mixing between different interaction eigenstates can occur 
in the same way as for the gauginos. In general the 3x3 trilinear coupling ma- 
trices and mass matrices in the soft supersymmetry breaking part of the MSSM 
Lagrangian can be fully occupied, thus leading to mixing between the sfermions 
of different generations. Such mixing results in contributions to flavor changing 
neutral currents (FCNCs) besides the contribution of the CKM-matrix which is 
already present in the Standard Model. Experimental limits [17] show that such 
additional contributions have to be small [65]. Additionally, most popular mod- 
els of supersymmetry breaking mediate this breaking from the hidden sector by 
flavor-blind interactions. Therefore the soft breaking mass matrices and trilinear 
couplings are chosen purely diagonal. Then the mass matrices of the electron-like 
sleptons and the squarks decompose into 2x2 matrices where only the left- and 
right-handed fields of each generation mix. These can be written as 






Mf= \ f ;? 'V^ ^ (3-72) 



with 



r r r, / f o \ I M^ for left-handed sleptons 

Mf =™|(/3'-Q,4)c., + ^^ fo. left-handed =qua.kB ^'-''^ 

M| for right-handed electron-like sleptons 
Mf^ =m\ {QfSw) C2/3 + -j M? for right-handed up-like squarks (3.74) 

M? for right-handed down-like squarks 

r „ , 7- for up-like squarks 

' tp for electron-like sleptons and down-like squarks 



(3.75) 



Qf is the electromagnetic charge of the sfermion. J3 denotes the quantum number 
of the third component of the weak isospin operator T3 which is -|-^ for up-like 
squarks and — | for down-like squarks and electron-like sleptons. These mass 
matrices can again be diagonalized by a unitary matrix 

fm^~ \ 
Uf-MjUJ = I /^ _2 1 • (3.76) 

The fields then transform as 

^'^ - Uf f f ) . (3.77) 




f2 ' V/ 



'R 



3.6. Particle content of the MSSM 29 



In the sneutrino sector only left-handed fields exist. For this reason the mass ma- 
trix consists only of the M^^ element given in eq. (3.73). M^^ is therefore a free 
parameter of the theory which directly gives the squared mass of the sneutrinos 
according to 



m- 



2 _ 1_2^ , 7,^2 



-m',C2p + Ml . (3.78) 



ui 2 Z '^'^ L 



The interaction eigenstates uj are identical to the mass eigenstates. 



Chapter 4 

Regularization and 
Renormalization 



In general the Lagrangian of a model contains free parameters which are not 
fixed by the theory, but must be determined in experiments. On tree-level these 
parameters can be identified directly with physical observables like masses or 
coupling constants. If one goes to higher-order perturbation theory these rela- 
tions are modified by loop contributions. Additionally the integration over the 
loop momenta is generally divergent which further complicates the situation. To 
achieve a mathematically consistent treatment it is necessary to regularize the 
theory before predictions can be made. This introduces a cutoff in the relations 
between the parameters and the physical observables. As a consequence, the 
parameters appearing in the basic Lagrangian, the so-called "bare" parameters, 
have no longer a physical meaning. This physical meaning is then restored via 
renormalization. The renormalized parameters obtained in this way are again 
finite. Their value is fixed by renormalization conditions. 

The details of this procedure are described in the following sections. 

4.1 Regularization 

The ultra-violet divergences appearing in the integration over loop momenta must 
be treated via a regularization scheme. Therefore a regularization parameter 
is introduced into the theory which leads to finite expressions, but leaves the 
expressions dependent on the renormalization parameter. 

There exist different regularization schemes, three of which are shortly de- 
scribed in the following: 

Pauli-Villars Regularization 

This regularization scheme [66] is very simple. Originally the integration region 
over the four- dimensional loop momentum ranges from plus to minus infinity. In 
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this scheme it is restricted such that the absolute value of the loop momentum 
is below a certain finite value. This cutoff parameter must be much larger than 
any other mass scale appearing in the theory. Performing a regularization in this 
way usually destroys gauge symmetry, so it is not used for practical calculations 
and not further taken into account in this dissertation. 

Dimensional Regularization 

Loop integrals are divergent if the dimension of the integration is exactly four. 
Dimensional regularization [67] exploits this fact. If one shifts the dimension of 
the loop momentum by an infinitesimal value and performs the integration in 
Z? = 4 — 2e dimensions, the integral becomes finite. The divergences now appear 
as poles in the infinitesimal parameter e. Additionally, the dimensions of all fields 
are also set to D dimensions and the gauge couplings are multiplied by /i^'^. The 
parameter /i has the dimension of a mass and specifies the regularization scale. 
It is introduced to keep the coupling constants dimensionless. This scheme is 
normally used in Standard Model calculations as it preserves gauge symmetry. 
It does, however, not preserve supersymmetry. As the fields are D- dimensional, 
additional degrees of freedom are introduced so that the number of fermionic 
degrees of freedom no longer equals the number of bosonic degrees of freedom 
and therefore supersymmetry is broken. 

Dimensional Reduction 

This scheme [68, 69] is similar to dimensional regularization in the respect that the 
loop integration is performed in D dimensions and the divergences are recovered 
as poles in e. In this scheme the fields are kept four-dimensional in order to avoid 
explicit supersymmetry breaking. The mathematical consistency of dimensional 
reduction has long been questioned [70], but recently a consistent formulation [71] 
could be established. It could be shown that supersymmetry is conserved for 
matter fields at least up to the two-loop order. 

4.2 Renormalization 

The dependence on the unphysical scale fi can be removed via renormalization. 
It consists of a set of rules which consistently replaces the bare parameters in the 
Lagrangian by new finite ones. 

There exist different degrees of renormalizability. One possibility are super- 
renormalizable theories. They are characterized by the fact that the coupling has 
positive mass dimension. In these theories only a finite number of basic Feynman 
diagrams diverge. These divergences can, however, appear as subdiagrams at 
every order in perturbation theory. An example for such a theory is scalar 0"^- 
theory. Here apart from vacuum polarization graphs only the one- and two-loop 
tadpoles and the one-loop self-energy diagram are divergent. 
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In renormalizable theories only a finite number of amplitudes diverge, but 
tliese divergences occur at all orders of perturbation theory. In such theories 
there are also dimensionless couplings but none with a mass dimension smaller 
than zero. To cancel the divergences a finite set of rules is necessary. Non-Abelian 
gauge theories like the Standard Model and the MSSM belong to this category. 
Their renormalizability was first proven in ref. [72]. 

Finally a theory can be non-renormalizable. In this case all amplitudes are 
divergent if the order of perturbation theory is sufficiently high. The set of 
rules to absorb the divergences is infinite and new ones appear at each order of 
perturbation theory. This means that the theory loses its predictive power. It 
might at first sight look like such models would be useless, but this is not the case. 
Non-renormalizable models are often used for effective theories. Here operators 
of a mass dimension greater than four appear in the Lagrangian. As the final 
expression in the Lagrangian must be of mass dimension four, this is compensated 
by an appropriate power of a cut-off mass appearing in the denominator. This 
cut-off mass defines the energy scale up to which the effective theory is valid and 
above which it must be replaced by the full renormalizable theory. In the overlap 
region where both theories give a useful result, a matching between the two is 
performed, thereby fixing the renornialization conditions and allowing meaningful 
predictions. 

4.2.1 Counter terms 

One of the most popular renormalization approaches nowadays is multiplicative 
renornialization with counter terms. In this scheme the bare parameters go, i.e. 
couplings and masses appearing in the Lagrangian, are replaced by renormalized 
ones g, which are related to the bare ones via the renormalization constant Zg 

go = Zgg={l + 6ZJ^'^ + 6Zp + ...)g , (4.1) 

where on the right-hand side the renormalization constant has been expanded in 
orders of perturbation theory and the order is denoted by the superscript. The 
renormalized g have a finite value. The SZg absorb the divergences which ap- 
pear in the loop integrals and are parametrized in the regularization parameter. 
Therefore they remove the dependence on the unphysical regularization parame- 
ter. Additionally, finite parts can be absorbed in the renormalization constants, 
as the decomposition in eq. (4.1) is not unique. Which finite parts are absorbed 
in the renormalization constants depends on the chosen renormalization scheme, 
which will be discussed below. If one also adds the wave function renormalization 
of external particles, the renormalization of the parameters is sufficient to obtain 
finite S-matrix elements. To achieve the finiteness of off-shell Green functions, 
the fields must be renormalized as well. Therefore the bare fields $o are replaced 
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by the renormalized ones $ and the field renormahzation constant Zq, 

% = \fzl^ = A + ^(5Z« - ^<5zi^)' + \bzf + ...") $ . (4.2) 

Again on the right-hand side the field renormalization constant is written out as 
an expansion in orders of perturbation theory. Thereby, the term containing the 

squared of the one-loop renormalization constant (— |5Z^ ) is part of the two- 
loop contribution. Similarly, for higher orders the orders of all renormalization 
constants which appear in a term must be added up to give the loop order to 
which the term contributes. 

Using both parameter and field renormalization all Green functions are fi- 
nite. We can now insert the renormalized parameters and fields into the bare 
Lagrangian 

L {go, $o) = C (^Zgg, v^$) = C {g, $) + Cct {g, $, Zg, Z^) (4.3) 

and write it as a sum of the renormalized Lagrangian C {g, $) and the counter- 
term part which can be expanded in terms of the loop order 

Cct (g, $, Zg, Z^) =/:gi. [g, $, 5Zf\ ^4'^) + 

C% [g, $, 5ZW, 5Z«, 5Zf , 5Zf ) + . . . . ^^'^^ 

In this thesis only one-loop corrections to processes are considered. So only 
the one-loop counter terms 5Z^^^ enter the calculations, hence for simplicity the 
superscript (1) on the 6Z will be dropped from now on. 

4.2.2 Renormalization Schemes 

The finite part of the renormalization constants is not fixed by the divergences, 
but can be chosen in a suitable way. The definition of these finite parts together 
with an independent set of parameters comprises a renormalization scheme and 
therefore defines the relation between the observables and the parameters of the 
theory. If one adds up all orders of perturbation theory the result is indepen- 
dent of the chosen renormalization scheme. The value of the input parameters, 
however, still depends on the renormalization scheme and must be chosen appro- 
priately. For actual calculations only a finite number of orders can be taken into 
account. The resulting dependence on the renormalization scheme is then a mea- 
sure for the theoretical uncertainty which is induced by the missing higher-order 
terms. 

The simplest renormalization scheme is the minimal-subtraction scheme or 
short MS-scheme [73] . It is based on dimensional regularization as regularization 
scheme. In this scheme the counter terms absorb just the divergent --terms but 
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no finite contributions. This scheme is actually a whole set of schemes, as the 
scale fi, which was introduced in the regularization step, is still present. This 
scale is now taken as the renormalization scale fiR and for specifying a concrete 
renornialization scheme fiR must be fixed as well. 

A commonly used variant of the MS-scheme is the modified minimal subtrac- 
tion scheme or short MS scheme [74, 75, 76]. It is based on the observation that 
the i-terms are always associated with other constant terms that emerge from 
the continuation of the loop momentum in D dimensions and are denoted by A", 
where n is the loop order. At one-loop order it has the following explicit form 

A = -fE + lnAn , (4.5) 

where --fE denotes the Euler-Mascheroni constant. The absorption of the numer- 
ical constants 7^ and In Att corresponds to a redefinition of the renormalization 
scale 

f,l''':=^,y-'----. (4.6) 

If dimensional reduction is used as the regularization scheme, the renormal- 
ization scheme is called DR. Apart from that the procedure is identical to the 
MS scheme. The A" terms are subtracted by the renormalization constants but 
no other finite parts. As before, this corresponds to a redefinition of the renor- 
malization scale /i|. . On the one-loop level the counter terms are identical, 
while at higher orders they differ because the two regularization schemes induce 
different finite parts. 

Another, distinct possibility is the on-shell scheme (OS scheme) [77, 78]. The 
expression on-shell means that the renormalization conditions are set for particles 
which are on their mass shell. The mass of a particle which is on-shell is given by 
the real part of the pole of the propagator and can be interpreted as its physical 
mass. In the OS scheme the real parts of all loop contributions to the propagator 
pole and consequently to this mass are absorbed in the mass counter terms. 
Hence the counter terms in the OS scheme also have a non-vanishing finite part 
and the dependence on the regularization scale /i is completely eliminated in this 
scheme. Coupling constants are renormalized in the OS scheme by demanding 
that the coupling constants stay unchanged if all particles coupling to the vertex 
are on-shell. This means that all corrections to the vertex are compensated by 
the counter term of the coupling constant. For the on-shell renormalization of 
fields one demands that the propagators are correctly normalized, i.e. the residue 
of the renormalized on-shell propagator is equal to one. 

The renormalization of tp, the ratio of the two Higgs vacuum expectation val- 
ues, is performed via DR also when otherwise the OS scheme is used [79]. As tp 
does not receive any SUSY-QCD corrections at one-loop order, its renormaliza- 
tion is not necessary for the calculations of this thesis. Also the strong coupling 
constant ag is always renormalized in the MS or DR scheme. The details of the 
renormalization of as are presented in the next section. 
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A complete expression of all Standard Model one-loop counter terms in the 
OS scheme was given in ref. [78]. Its extension to the MSSM was performed in 
ref. [80] . In this thesis the same conventions as in these two references are used. 



4.2.3 Renormalization of the strong coupling constant 

As every other parameter, the coupling constant Qs of the strong interaction re- 
ceives divergent loop corrections. These divergences must be removed by renor- 
malization. As shown in the previous section, g^ = ^/A7^a^ is renormalized mul- 
tiplicatively such that 

9! = Z,^9s'''' = '''il + 5Zjg, . (4.7) 



The explicit form of SZg^ depends on the renormalization scheme. Choosing 
the OS scheme for this task, however, is not possible. If the renormalization 
condition for Qs is formulated completely analogous to the renormalization of the 
electromagnetic coupling constant, one must demand that the corrections to the 
gluon-quark-anti- quark vertex vanish in the limit of zero-momentum transfer. 
To formulate this condition the value of Qs would be needed in a region which 
is below the QCD scale Aqcd- Coming from values above, gs formally reaches 
infinity at this scale and perturbative methods are no longer defined. As the 
OS scheme is based on the validity of perturbation theory this would lead to a 
self-contradiction. 

Instead another renormalization scheme must be used, which avoids the de- 
pendence on Qs at zero-momentum transfer. The MS and DR schemes share this 
property. In these schemes the counter term 6Zg^ is fixed by the condition that 
the gluon-quark-anti-quark vertex is finite. Due to a Slavnov- Taylor identity, 
which guarantees the universality of Qs, this automatically results in finite three- 
and four-gluon vertices. The counter term has the following explicit form 



6Zg^ = -^in-^nf-2--nf]A , (4.^ 



where the contributions to the sum originate from gluons, quarks, gluinos and 
squarks. rif = Q denotes the number of quark flavors. The last two terms originate 
from the supersymmetric particles and are not present in the Standard Model. 

The behavior of Qs with respect to higher-order corrections can be improved by 
the use of renormalization group equations (RGE). The one- loop RGE sum up all 
leading-log contributions which have the form g^"" {fin) (In/i/j)". Their application 
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leads to the following expression for the strong coupling constant^ [16]: 

a^M = -^S^ . (4.9) 



1 - f aDR (m^) In ^ 



The experimental value [17] for a^ is given in the MS scheme at the scale mz and 
using the Standard Model RGE for extracting a^^ from the data. This must be 
converted to a^^ via 

ar(™z) = f^ (4.10) 

1 - Aa^ 

with 



2n \ 2 3 mz mz 6 ^^-^ V rra^ m^z 



(4.11) 



The In ^^ terms in the last equation decouple the particles heavier than m^ from 
the running of a^. 

Also the finite part of the one-loop contribution to the gluon-quark-anti- 
quark vertex depends on the renormalization scale fiR. It should best be chosen 
in a way that the error, which is induced by missing higher-order corrections, is 
as small as possible [81, 82]. Since i?-parity is conserved in the MSSM, the one- 
loop diagrams decompose into two distinct sets, where the loop either consists 
solely of SUSY particles or does not contain any supersymmetric particles at all. 
The latter ones form the corrections which also appear in the Standard Model. 
Except for the top quark, which is decoupled, they take part in the running of 
as- For these contributions the same renormalization scale fi^ should be used as 
in eq. (4.9) which is typically of the order of the energy scale of the considered 
process. 

For the additional SUSY contributions another, special value fiR is chosen [81, 
82]. This is possible because the two sets of diagrams are distinct and all super- 
symmetric particles are decoupled from the running of a^. The scale is chosen 
such that the contribution of these diagrams vanishes at zero-momentum transfer. 
Under this condition gs is taken at the scale fiR, so the procedure is well-defined. 
It is fulfilled if 

^^ squarks ^ ^^R h^R / 



^Here the formula for Us is quoted as this constant is normally used in calculations and also 
the experimental value of the coupling constant is given in terms of as- The corresponding 

2 

expression for gs can simply be derived from the relation as = |^. 
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Figure 4.1: One-loop SUSY-QCD diagram mediating an effective coupling be- 
tween the bottom quark and H^. The cross in the gluino line here represents 
a mass insertion, i.e. the rrig term is chosen when computing the trace over the 
fermion line. The subscripts a and b of the sbottom particles take the values 1 
and 2. 



Solving for fifi yields 



IJ'R = Vrn^ Y\_ {^qi^'12)^'' ■ (4-13) 

squarks 



This procedure reduces the numerical value of the one- loop corrections and there- 
fore makes the calculation more stable against the theoretical uncertainty from 
missing higher-order terms. 



4.3 Bottom-quark Yukawa Coupling 

The mass of the bottom quark and its Yukawa coupling to the Higgs particles 
are intimately related. They originate from the same term in the unbroken La- 
grangian. After the Higgs fields have acquired a vacuum expectation value, the 
vacuum-expectation-value component yields the mass term of the bottom quark 
in the Lagrangian and the other components describe the Yukawa coupling of the 
bottom quark to the various Higgs particles. This relation can be modified by 
loop corrections, and it turns out that these are very large in the case of bottom 
quarks [7, 8]. A resummation of the leading corrections to all orders in perturba- 
tion theory can be performed, which greatly reduces the theoretical uncertainty 
originating from unknown higher-order corrections. 

At tree-level the bottom quark only couples to the first Higgs doublet Hi as 
can be seen from the superpotential eq. (3.38). A coupling to the second one 
H2 is forbidden. Such a coupling can, however, be generated dynamically at the 
one-loop level. Taking into account only SUSY-QCD corrections, i.e. corrections 
with squarks and gluinos, this is done by the single diagram Fig. 4.1. Although 
this contribution is loop-suppressed, it can induce a potentially large shift in 
the tree-level relations, because it is enhanced by tj^. By electroweak symmetry 
breaking the Higgs field H^ acquires a vacuum expectation value V2 and firstly we 
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will consider only this part. On tree-level the bottom-quark mass and its Yukawa 
coupling Xjj are related via 

TUb = XbVi. (4.14) 

Adding the vacuum-expectation- value contribution from Fig. 4.1 changes this 
equation to 

nib = \bVi + AXbV2 = Vi{\b + AXbtf^) = \bVi{l + Arrib) ■ (4.15) 

As the numerical value of nib is fixed by experiments, this results in a change of 
the effective Yukawa coupling of the bottom quark 

nih 1 

Aft = —-, T • 4.16 

ui 1 + Amb ^ ^ 

Computing the diagram in Fig. 4.1 in the limit of vanishing external momen- 
tum yields the following explicit form for Anib- 

^fTT-b = —^nigfitfsl {nif^^,ni^^^,nig) (4.17) 



with 
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^ ("^6i , ^b2 ^^g) =- [ ^l rnl In -X + ^I'^l ^'^Z;i + ^l^l ^^ zJt 



mf "^ ^ m^ ^ "1 m~ 
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X -, ^-. ^- ^-. ^ (4.18) 



and /i denoting the MSSM parameter which couples the two Higgs doublets. In 
the limit where the squark and gluino masses have approximately the same value, 
denoted by a common SUSY mass ttisusy? the last equation simplifies to 

/("^SUSY,msUSY,rnsUSY) = 7^-^ • (4-19) 

If additionally /x is of comparable size, this results in 

Am6 = sign(/i) tp . (4.20) 

on 

So for large values of tp this effect can be of (9 (1) and does not vanish for heavy 
SUSY spectra. 

For computations up to one-loop order eq. (4.16) can be expanded so that it 
contains only corrections up to O (a^). The equation is then modified and reads 

A, = ^(l-Am,) . (4.21) 

Vl 
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So for large absolute values of Anih, which are phenomenologically very inter- 
esting, huge one-loop corrections appear. If Arrif, exceeds one, the standard way 
of computing one-loop cross sections by adding the interference term between 
tree-level and one-loop diagrams even yields negative total cross sections which 
are obviously wrong. One might even question if perturbation theory is still valid 
in this regime, but definitely higher-order calculations would be needed to reduce 
the theoretical uncertainty. 

This problem is solved by the observation that these corrections do not appear 
at higher orders. In ref. [7] it was proven that there are no contributions to Arrih 
of 

o((as-^^t^ ] (4.22) 

for n > 1. Higher-order corrections either lack the enhancement factor tf^ or are 
suppressed by a mass ratio ^"^'' . Therefore Arub is a one- loop exact quantity 
and including it as in eq. (4.16) contains the corrections to all orders in Og which 
have the form given in eq. (4.22). 

Using the resummed form eq. (4.16) is only useful when computing total cross 
sections. For a comparison with one-loop cross sections it is necessary to use 
eq. (4.21) so that the same order in ag is taken into account in both calculations. 
This will be explained in more detail in chapter 6, where this procedure is applied 
to a physical process. 

The Ami) corrections are universal. They occur in every coupling of the 
bottom quark to the different Higgs particles, both neutral and charged ones. 
They are also independent of the kinematic configuration. 

When the bottom quark couples to the physical Higgs fields an additional 
term occurs. It also originates from diagram Fig. 4.1, but now not the coupling 
to the vacuum expectation value but to the remaining neutral Higgs field (j)2 is 
considered. In addition to the tree-level coupling to the first Higgs doublet 

^bbcho = Xb = — (4.23) 

this induces another term 

r,s,. = A.AA. = ^^ , (4,24) 

' Vl tp 

After electroweak symmetry breaking the fields (f)\ and (f)\ must be rotated by 
the angle a to form the two CP-even mass eigenstates h^ and H^. Combining 
everything this leads to the following effective couplings of the bottom quark [7, 
83] 

Arr ib 
Arrib \' tfsta 

'^--rj,„,.^fl + ^) . (4.26) 



h;*. =rL. N . .„ I 1 - TT^ 1 (4.25) 
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where r° denotes the respective tree-level coupling. Expanding these equations 
up to the one-loop order yields 

rko =rLo (l - Am, (l + ^) ) (4.27) 

^IHO =r^5H0 (l - ^rn, (l - I) ) • (4.28) 

In the coupling of the top quark to the Higgs fields a similar effect occurs. On 
tree-level the top quark couples only to H2 and a coupling to the second doublet 
Hi is generated perturbatively. This results in a modified Yukawa coupling which 
is given by 

fTLf 1 , , 

A. = — ,^. , 4.29 

in complete analogy to eq. (4.16). The correction term Arrit has the form [84] 

Amt = —^rugfi—l (m^-^ , mi^ , Trig) . (4.30) 

OlT tf3 

In contrast to Am, this equation has a suppression factor of j-. Therefore its 
numerical impact is much smaller than the t/3-enhanced bottom-quark correction 
and it is largest for small values of tj^. The contribution of this correction is 
nevertheless significant and therefore it is justified to include its effect in the 
same way as for the bottom-quark correction. 

Also the coupling of the top quark to the physical Higgs particles gets an 
additional contribution from the coupling to the Hi doublet. In this case the 
modified couplings are 

^'"'' ^^°*"° 1 + Am, ^^ " ^"^*M«) (4-31) 

where r° denotes the respective tree-level coupling. An expansion up to the 
one-loop order yields 

r,\,o =r°f,„ (1 - Am, (1 + tpt^)) (4.33) 

^Imo =T'uHO (1 - Am, (1 - ^) ) • (4.34) 



Chapter 5 

Hadronic Cross Sections 



The cross sections which are obtained by applying the Feynman rules contain, 
amongst other particles, quarks and gluons. The leading interaction between 
these particles is the strong interaction, which is described by quantum-chromo 
dynamics (QCD). This theory possesses two characteristic properties: asymptotic 
freedom [85] and confinement. Asymptotic freedom describes the behavior of the 
theory at small distances. In this region the interaction is weak and the coupling 
constant gets smaller with decreasing distance or, equivalently, with rising energy. 
At large distances confinement appears, because the interaction becomes strong 
and binds the particles tightly together. If the space between them becomes 
even larger, it is energetically favorable to form new quark-anti-quark pairs. One 
consequence of this behavior is that quarks and gluons cannot be observed as free 
particles, but only as constituents of hadrons, i.e. mesons, which are quark-anti- 
quark pairs, and baryons, which are states of three quarks or three anti-quarks. 
An example for these hadrons are protons, which are the colliding particles at the 
LHC. To make theoretical predictions it is necessary to relate the interactions at 
the parton level to the interactions at the hadron level [86]. The basis for doing 
this is the parton model [87], which will be described in the next section. 



5.1 Parton Model 

The parton model describes the inner structure of hadrons in hard collisions. It 
starts from the assumption that every observable hadron consists of constituents, 
the so-called partons, which can be identified as quarks and gluons. Experimental 
evidence for this assumption comes from the observation of scaling [88] in deep 
inelastic electron-proton-scattering. If the hadron carries some momentum P'^, 
the partons which take part in the partonic subprocess have momentum xP'^ 
with a; G [0, 1]. As normally the mass of the hadrons is small compared to their 
kinetic energy one can assume P^ = 0. 

The interaction of an electron and a hadron or of two hadrons among each 
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other can be split into two parts. Because of Lorentz contraction and time dilation 
the interaction time of the two incoming particles in the laboratory frame is very 
short. Therefore effectively a static hadron is seen. For the hard scattering 
process interactions between partons of the same hadron need not be considered. 
Also the process of hadronization after the interaction happens on time scales 
which are much larger than the interaction itself. 

From this the theorem of factorization [89] follows immediately. It states that 
all diagrammatic contributions to the structure functions can be separated into 
a product of two functions C and /, which depend on two mass scales fiR and 
Hf- l^R is the renormalization scale which was already defined in chapter 4, ^p is 
the so-called factorization scale and separates the long-distance from the short- 
distance effects. Slightly simplifying one can say that every parton propagator 
which is off-shell by ^p or more contributes to C, while those which are below 
this value contribute to /. 

5.2 Integrated Hadronic Cross Sections 

The hard scattering process C therefore can be calculated in perturbation theory 
by Feynman rules, using partons as incoming particles. It is independent of 
long-distance effects and especially from the type of the colliding hadron. 

The parton distribution function (PDF) fi/h{x, ^p) contains the long-distance 
effects. It is independent of the underlying scattering process, but depends on ^p 
and the type of hadron h. It is normalized such that it can also be interpreted as a 
probability density, namely the probability of finding the parton i in the hadron h 
with a momentum xP^. Its behavior as a function of the parameters is determined 
by the Altarelli-Parisi integro-differential equations [90]. Its numerical value, 
however, cannot be calculated a priori from the theory. At a single reference 
point it must be determined by experiments. 

Therewith one obtains the expression [86] 



<ypp~*fin+X = 2Z / '^'^~I~^rnn-*f^n{JS,as{^^R)) (5.1) 

r .. .1 J Tn ^ 



"^ , d£ 

{m,n} "^ 

for an integrated hadronic cross section with the parton luminosity 
d£ r^ dx 1 



-"mn 



dr J^ X 1 + 6„ 

■ (jm/p {X, flp) fn/p (^-, fipj + fn/p {x, fJ^p) fm/p (^-, I^fJJ ■ (5.2) 

Here vo denotes the hadronic center-of-mass energy, i.e. the one of the two 
colliding protons, and dmn^fin the partonic cross sections of the subprocesses, 
where the two incoming partons m and n produce some final state, labeled fin. 
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The sum includes all possible parton combinations ni and n where the order of 
appearance is not taken into account. The integration variable r relates the par- 
tonic and hadronic center-of-mass energies with each other. More specifically, ^/T 
can be interpreted as the part of the hadronic center-of-mass energy which takes 
part in the partonic subprocess, as the partonic center-of-mass energy is given 
by vl = \^tS. The lower limit of the integral Tq is determined by the kinematic 
configuration. \/tqS is the minimal energy which is necessary to produce the 
final state fin and therefore denotes the production threshold. 

The formula given above is valid for processes with two or more particles in the 
final state. For hadronic cross sections it is also possible to calculate integrated 
cross sections for 2^1 processes. One first obtains for the partonic cross section 
of the process mn -^ f 

da_^; = "l \Mf, {mn ^ /) \'6 {pi + pi - p]) . (5.3) 

'iPf\/s\Pra\ 

Again m and n specify the incoming partons, / denotes the outgoing particle, 
ntf its mass, and p° the energy of the respective particle i. Pm indicates the 
three-momentum of particle m in the partonic center-of-mass system and Ai/i is 
the matrix element. 

When convoluting with the parton distribution functions the single remain- 
ing 5-function in the partonic cross sections solves the r integral in eq. (5.1) 
analytically. Thus one obtains for the integrated hadronic cross section 
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5.3 Differential Hadronic Cross Sections 

Additionally one can define hadronic cross sections that are differential in one 
or more parameters. For these parameters it is useful to take variables that are 
either invariant under Lorentz transformations or at least have very simple trans- 
formation properties. In this thesis three differential hadronic cross sections are 
presented which are also implemented in the HadCalc program that is described 
below in section 5.5. They are cross sections differential with respect to the in- 
variant mass of the final-state particles, the rapidity of one final-state particle 
and, thirdly, the transverse momentum. 

5.3.1 Invariant Mass 

The first differential hadronic cross section is the one with respect to the invariant 
mass of the final-state particles. The invariant mass of a process is equivalent to 
the partonic center-of-mass energy Vs = VtS of the process or, in other words. 
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the sum of the final-state momenta of the outgoing particles. The differential 
cross section takes the form 



{m,n\ s 

where fin again labeles a general final state. 

5.3.2 Rapidity 

The rapidity y of a particle is defined as 

y = artanh — = - In — (5.6) 

where pz = p ■ cg denotes the fraction of the particle's three-momentum p that 
goes in the direction of the beam axis, labeled z. The mass of the particle will 
later be referred to as m. Using the rapidity instead of directly taking the angle 
9 between the particle and the beam axis possesses some advantages because the 
rapidity of a particle has a few useful properties. Under a boost in the 2;-direction 
to a frame with a velocity j3, the rapidity transforms as y ^ y — artanh f3. Thus 
the shape of the rapidity distribution ^ stays unchanged. More generally, the 
sum of two rapidities when the momenta point into the same direction is given by 
the rapidity of the sum of the momenta, added via the formula for the relativistic 

addition of velocities: y {pi)+y (^2) = y ( ff^^^ ) . In experimental analyses often 
a slightly different measure, the pseudo-rapidity 77, is used. It is derived from the 
standard rapidity by taking the limit of a vanishing mass of the particle and is 
defined as 

r]=-\n- . (5.7) 

Z i — C0 

In the HadCalc program both normal rapidity and pseudo-rapidity are imple- 
mented. As conversion between both variables can be performed by the simple 
transformation 

y = artanh Wl 2 tanh?7 1 , (5.8) 

y y p +m? j 

in the following only the shorter expressions for the standard rapidity are given. 
The ones for pseudo-rapidity can then be deduced from them. 

Using the above-mentioned definition of the rapidity the differential hadronic 
cross section with respect to the rapidity for 2^2 processes then reads 

da _ /"^ dC d(T dcg 
dy y^g dr dc^ dy 
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The momenta and masses given in the formulae always refer to the particle for 
which the rapidity distribution is calculated. The angle Cq between the particle 
and the beam axis in the partonic center-of-mass system is fixed by the relation 



/ Tn^ I 1 X \ 

c,~= Jl + ^tanh(y + -ln — j (5.10) 

where the second term in the argument of tanh originates from the boost from the 
hadronic center-of-mass system, which is the laboratory frame, to the partonic 
one, in which the partonic subprocess is calculated. This leads to 

''' ■ - (5.11) 




dy y p^ cosh^ (y + I In ^) 



For processes with three or more particles in the final state the formula is 
very similar. Additional phase-space integrals appear for the further particles 
but otherwise eq. (5.9) stays unchanged. In the following equation the differential 
cross section for a 2 ^ 3 process is given 

dy ^ dr J 'J '7 ' dkldcgdk^^dr] dy 
The parametrization of the three-particle phase space is described in appendix C.2. 

5.3.3 Transverse Momentum 

The last implemented differential hadronic cross section is the one with respect 
to the transverse momentum pt = \/v\ + v\ of one of the final state particles. 
For 2-^2 processes it is defined as 

^^ ^'dr^^l^ (5.13) 



with 



dvT J TO dr dc^ dpT 
dcn 






(5.14) 



which follows from 



c^± = ±W1-t4 . (5.15) 

Here two possible solutions arise because of the sign ambiguity when taking the 
square root. In principle both solutions have to be taken into account and added 
up unless they are excluded by other constraints as shown below. The lower limit 



48 Chapter 5. Hadronic Cross Sections 

of the r-integral Tq must be adjusted such that Sg is always inside its co-domain 

[0;i] 



m\ +Pt + \/^h + Pt) 



fo = ^^^ ^ '- , (5.16) 

/i and /2 denoting the two final state particles. 

For 2 — i> 3 processes the extension to include the third final-state particle is 
straightforward. The lower limit for r in these processes is 



f ^J^j, + Pt + yi^.f2 + ^hf + Pr 1 



^0 = ^ ^ '— , (5.17) 

when the cross section is differential in the particle /i. Therefore the expression 
for the differential cross section reads 

^= [\r^ Ukl Ukl fdv ^^J'. ^ . (5.18) 

dpT 4 dry ^J V ' dkldc^dkldT] dpT 

5.4 Cuts 

In order to improve the ratio of the signal-process cross section to that of the 
background processes it can be useful to place appropriate cuts on the final-state 
particles. Also experimental techniques used in the reconstruction of events like 
jet-clustering algorithms can mandate the use of cuts in theoretical predictions, 
so that the behavior of these techniques is emulated there. 

In the HadCalc program cuts on three different properties of the final-state 
particles are implemented [91]. The first two are cuts on the rapidity and the 
transverse momentum of a particle. The definition of these two variables was 
already presented in the previous section. The third one is a mutual property of 
two particles, the jet separation ARij, which is defined as 



AR,, = J Ayf^ + A<Pl . (5.19) 



Ai/ij denotes the rapidity difference between the two particles i and j and Acpij 
the difference in the azimuthal angles of the two particles in the transverse plane. 
Its main use are exclusive hadronic cross sections where final-state jets shall be 
observed explicitly. It mimics the behavior of jet-clustering algorithms. There 
two jets, which are separated by a jet separation below a certain limit, are seen in 
the reconstruction as a single jet which has kinematic properties that are averaged 
over the two final-state partons. 

For the first two cut parameters, rapidity and transverse momentum, it is 
possible to translate these cuts into a limit on the integration parameters of 
the phase space. The most general case is assumed here that cuts on both the 
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rapidity i/cut and the transverse momentum PTcnt of ^ particle shall be applied. 
Using eq. (5.15) the transverse- momentum cut can be translated into a cut on Cg 
and one obtains 



c-^-^l-^<c,<^l-^ = c- . (5.20) 

Likewise, the cut on the rapidity can also be turned into a cut on c^ via eq. (5.10), 
yielding 



^e > W 1 + ?j tanh I -y^ut + - In ^ j = cj'^ 



P 
Ce<Jl + ^tanh U^, + 1 In ^"j = c^- . (5.21) 



To shorten the notation the abbreviation 



\ 



p 



1 + ^ 

p 



(5.22) 



is used in the following. Again the momenta and mass used in the equations all 
refer to the particle whose phase space should be constrained. 

Applying both cuts requires that the conditions on Cg are all fulfilled simul- 
taneously. This also restricts the integral on x which appears in the parton 
luminosity given in eq. (5.2). In total the x-interval divides into five different 
regions, which will be labeled by roman numbers. First the two cases where both 
cuts cannot be fulfilled simultaneously, are considered, because the lower limit of 
one cut lies above the upper limit of the other one: 
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Region V: 


/1 + r 
1 — r 


(5.24) 



These two regions are excluded and the cross section vanishes there. 

For specifying the other regions first two special cases are considered, where 
the lower limits on Cg and the upper limits, respectively, coincide. For these cases 
the according value of x is determined 

Oy Opj, 



max 






'l-r 


(5.25) 
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Using these two definitions the other intermediate regions can be specified. 
The ranges for Cg which are deduced from these following regions specify the 
allowed area where the cuts are fulfilled and therefore the cross section does not 
vanish. The next two regions handle the cases where the limits on Cq from rapidity 
and transverse momentum overlap and one limit is given by the rapidity cut and 
the other one by the transverse-momentum cut: 

Region II: cf^Kcf^Kc^Kcf'^^Kcf''^ => xj < x < mm{x^,^, x,^^,) 

(5.27) 
Region IV: cj"^ < cf" < c^ < cj^"" < cT"" => max(x^i,,x„,ax) < x < xy 

(5.28) 






Finally the definition of the last region is the case whether one cut gives a 
range on Cq that completely lies inside the other one. Depending on which cut 
this is, the limits on x are different: 



Region III a): cf" < cf ° < Cg < c|^^" < c^^" 



"pt "y " "y "pt 



(5.29) 



Region III b): cf^<cf^^<Cg<cf^^<cf^^ => x^ax < a; < a;„,„ . 

(5.30) 

In addition to those regions the original constraint for x for a hadronic cross 
section without cuts applies: 

r < X < 1 . (5.31) 

Combining the result of all regions one can see that no holes in the integration over 
X or Cg appear and the final borders of the integration routine can be simplified 
to 

max(r, xi) < x < min(xv, 1) (5.32) 

and 

max(cf ", cf °) < Cg < min(c^'^, c^") ■ (5.33) 

^ Bp^' By ' B ^ dp^' By ' ^ > 

For a cross section which is differential with respect to the rapidity of a final 
state particle the cut on the transverse momentum yields a restriction on Cg in 
the same way as in eq. (5.33) 

The constraint on x must then be adjusted such that Cg is always inside this 
allowed interval, yielding 



max(r, -yre"^-!/ ) < x < min(-\/re~^'i/ , 1) , (5.35) 

1 + r Vl — T 
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which corresponds to eq. (5.32) where the rapidity cut |/cut is replaced by its value 
y given as an argument to the cross section. 

Similarly, for cross sections that are differential in the transverse momentum 
of a final-state particle a cut on the rapidity puts a further constraint on the 
allowed interval for Cg^: 

cr < ce^ < cr (5-36) 



with 




(5.37) 
(5.38) 



Again this leads to a corresponding change in the limits of the x-integration which 
are given by 



' 1 fp I \ —I— f 

max(r, A/re~^""S/ z) <x < min(v/re^'^"S/ r, 1) (5.39) 

\/ 1 I ^ V 1 — r 



with 
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(5.40) 



p 

This again corresponds to eqs. (5.32) and (5.22) where instead of the cut on 
the transverse momentum PTcut i^s fixed value px, which is an argument to the 
differential cross section, is taken. 

5.5 HadCalc 

For the numerical evaluation of the cross sections presented in the following chap- 
ters a program called HadCalc was developed to facilitate this task. It is based 
on the established program packages FeynArts [9] and FormCalc [11, 12] which 
are used to generate the partonic cross sections. The main task of HadCalc then 
consists of the convolution with the PDFs that are taken from the PDFlib [92] 
or LHAPDF [93] library packages that include PDF fits from various groups. 

With this program it is possible to calculate both totally integrated and dif- 
ferential hadronic cross sections of processes with up to three particles in the final 
state. The latter ones can be differential with respect to the partonic center-of- 
mass energy, or the rapidity or the transverse momentum of one of the outgoing 
particles. Several cuts can be applied to the phase space. HadCalc operates 
either in batch mode, where the parameters are read from a file and the cross 
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sections are written back to disk, allowing for easy post-processing with e.g. a 
tool that generates plots. It can also be used in interactive mode where in- and 
output are done via keyboard and screen and which allows the user for example 
to tune the parameters most easily. 

A complete manual of HadCalc can be found in appendix F. The program 
code is available on request from the author^ 



"'^email: mraucMmppmu . mpg . de 



Chapter 6 

Associated Production of W H^ 



The discovery of a charged Higgs boson would be a clear signal of an extended 
Higgs sector and therefore of physics beyond the Standard Model. For relatively 
light charged Higgs bosons with a mass mn ^ f^t — rrib the main production 
process is tt-production via a subsequent decay sequence t -^ bH^ -^ br^Ur [94]. 
Both decay steps are enhanced by large Yukawa couplings. The experimental 
signature is an excess of rz/^ pairs in the detector. In the case of charged Higgs- 
boson masses above the top-quark mass the dominant production process is gb —>■ 
tH~ [95, 96, 97]. Afterwards the Higgs boson mainly decays into fet pairs with a 
branching ratio of at least 90%. The top-bottom-quark pairs lead to a detector 
signature which has a large QCD background at the LHC The detection of a 
heavy charged Higgs boson is therefore much more difficult. Later studies [98, 99, 
100] showed that the cross section is large enough so that the main decay channel 
can be ignored. It is sufficient to consider only the suppressed H~ -^ rv^ decay 
channel which has a clear detector signal while still yielding enough events. 

In this chapter we investigate another production mechanism, the production 
in association with a W boson. The leptonic decay modes of the W boson avoid 
large QCD backgrounds and can therefore provide an easier way of detecting a 
charged Higgs boson. 



6.1 The i/+Vl/- final state 

The production of a charged Higgs boson in association with a W boson was 
first studied in ref. [101]. This process proceeds either via bottom quark-anti- 
quark annihilation (Fig. 6.1) or via gluon fusion and an intermediate quark or 
squark loop (Fig. 6.2). The leptonic decays of the W boson could be used as 
a trigger for this process, thereby making charged Higgs boson detection easier. 
The calculation was updated in ref. [102] and triggered a detailed analysis [103] 
of the discovery potential at the LHC using this process. This paper concluded 
that an efficient separation of the signal process from the background processes 
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(a) s-channel contribution (b) t-channcl contribution 

Figure 6.1: Tree-level Feynman diagrams contributing to the dominant subpro- 
cess bb -^ H+W- 



such as top-quark pair production is difficult for semileptonic W boson decays 
including both low and high values of t^. The cross sections were evaluated at 
leading order for both production processes. 

The later studies of refs. [98, 99, 100] showed that a discovery is more likely 
in the main production channel gb — > tH~ where only the rare decay into a rUr 
pair is taken into account. Nevertheless the associated production of a charged 
Higgs with a W boson is an interesting process, especially when the existence 
of a charged Higgs boson has already been established before. If at that point 
no supersymmetric particles were detected, the question whether the H^ orig- 
inates from a Standard Model-like theory with an extended Higgs sector, like 
the Two-Higgs-Doublet Model (THDM), or from the MSSM remains open. In 
the latter case the cross section receives an additional contribution from virtual 
superpartners running in the loop. This can be used to tell the two models apart. 

As already mentioned earlier there are two important production processes for 
this final state in proton-proton collisions. The dominant one is the tree-level pro- 
duction (see Fig. 6.1) via bottom quark-ant i- quark annihilation. The s-channel 
diagrams shown in Fig. 6.1(a) are mediated by a virtual Higgs boson where all 
three neutral Higgs bosons of the MSSM (/i°, H^ and A^) can appear in the in- 
termediate state. The appearance of the massive particles in the s-channel leads 
to a propagator suppression of this diagram type. In the t-channel diagrams the 
exchange of a top quark occurs and yields the leading contribution to the bottom 
quark-ant i- quark annihilation process. As this class of diagrams contributes the 
most to the total H'^W" production rate one-loop corrections to this process are 
also important as they can modify the cross section significantly. Standard-QCD 
corrections of O (as) to this process were calculated in ref. [104]. Both the DR 
and the OS renormalization scheme were considered and good agreement between 
the two schemes could be found. The corrections are typically of O (15%) and 
can reach up to 30% in the small tj^ regime. The supersymmetric electroweak 
corrections, i.e. corrections where squarks together with charginos and neutrali- 

nos appear in the loop, of O (am^(.^^/m^ j and O [ctmtf^JTn^] were calculated 
in ref. [105]. 

The second parton process contributing to the H^W~ final state proceeds 
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(b) squark contribution 
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Figure 6.2: Leading-order types of Feynman diagrams contributing to the sub- 
process gg —^ H^W~ 



via gluon fusion and an intermediate loop as shown in Fig. 6.2. Since there is 
no tree-level process and the leading order contains a quark or squark loop it is 
suppressed by a factor a^ with respect to the bottom-quark annihilation process. 
The higher density of gluons in the proton partly compensates this effect, making 
both processes comparable in size. The contribution of quark loops was already 
included in the calculation of ref. [102]. The contribution of supersymmetric 
particles was calculated later on [106] and it was shown that they can reach up 
to 40%. Together with the QCD corrections this can raise the cross section for 
small tp so that it becomes comparable in size to the bottom-quark annihilation 
process which is not loop-suppressed. 

In this thesis the missing supersymmetric QCD corrections, i.e. corrections 
with squarks and gluinos running in the loop, to the leading bottom-quark anni- 
hilation process are considered. They are the last contribution of O {a^a^) to the 
associated production of a charged Higgs boson with a W boson in the MSSM 
which has not been calculated so far. 



6.2 SUSY-QCD corrections to hh -^ H+W' 

In this chapter the supersymmetric QCD corrections to the main production 
process bb -^ H~^W~ are calculated. As already shown in chapter 4.3 it is known 
that the coupling of the bottom quark to the Higgs bosons receives large one- 
loop corrections. These can be parametrized by introducing a correction term 
Anib to the bottom-quark Yukawa coupling. Yet other terms also give significant 
contributions, as we will see later. So it is necessary to not only use the tree- 
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(a) self-energy correction 
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Figure 6.3: Diagram types yielding SUSY-QCD corrections to the process bb 



level result with an effective bottom-quark Yukawa coupling but to perform a full 
one-loop calculation. 

The possible types of diagrams which appear in the calculation of SUSY-QCD 
corrections are depicted in Fig. 6.3. A SUSY-QCD self-energy contribution to 
the bottom-quark propagator of O (as) enters in the t-channel exchange diagram 
as shown in Fig. 6.3(a). Vertex corrections (Fig. 6.3(c)) appear in the s-channel 
diagrams in the vertex where the incoming bottom quark and anti-quark couple 
to the intermediate Higgs boson. The t-channel diagram receives vertex cor- 
rections at both the btW- and btH-Yeitices. Finally all four external particles 
can be connected via a box-shaped loop diagram (Fig. 6.3(e)). Additionally for 
the self-energy and vertex corrections appropriate counter-term diagrams appear 
(Fig. 6.3(b), (d)). 

The cross sections were calculated in both the OS and DR renormalization 
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schemes. Additionally a A?72fe-corrected tree-level cross section was calculated. 
As shown in chapter 4.3 the bottom-quark mass and the bottom-quark couplings 
to the Higgs fields receive large contributions from the one-loop SUSY-QCD 
corrections which are parametrized in the variable Ami,. To be able to compare 
the improved tree-level cross section with the full one-loop cross section it is 
necessary to use the same order in perturbation theory for both calculations. 
This means that one must use the non-resummed replacement eq. (4.21) 

rrif, -^ ruf, (1 — Anih) 

and the non-resummed correction to the bottom-quark Yukawa coupling as in 
eq. (4.27). A matrix element with this replacement must be treated as a one-loop 
matrix element. Let us recall that the standard way of computing a one- loop cross 
section is to add the interference terms {T*L + TL*) to the tree-level cross section 
\T\ and to discard the loop-squarei 
that for the squared matrix element 



T\ and to discard the loop-squared term \L\ which is a two-loop quantity so 



\Mf^\'^ = \T\'^ + 2Re{T*L) (6.1) 

is obtained. In this equation T denotes the tree-level amplitude and L the ampli- 
tude of the one-loop diagrams. In complete analogy the squared matrix element 
of the Am(,-corrected cross section is defined as 



^fiAn., 



\Tf + 2ReiT*LAnJ (6.2) 



with 

LAm,, = Tj\rni, — T . (6.3) 

Ti\mi denotes the tree-level cross section with the Ami, replacements of eq. (4.21) 
and eq. (4.27). Therefore L/^rm, contains the additional contribution which orig- 
inates from the correction terms. The corresponding cross section to 

is denoted by o"a in the following. 

In order to present the numerical results several relative corrections A using 
various contributions are defined. Firstly there is the relative correction in the 
OS scheme, 

OS _ OS 

Aos-'-^-oP^ . (6.4) 



^f^Am, 



The relative correction in the DR scheme is defined analogously as 



,DR _DR 
"^0 



^Dn= ' ' . (6.5) 



a?^ 



The third relation consists of the difference between the one-loop result and the 
A?72{,-corrected tree-level result which is calculated according to eq. (6.2). Hence 
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it signifies the true one-loop corrections. It is defined as 

Aa^,= ' os^ ■ (6-6) 

The subscript of the cross section a always denotes the order of the respective 
loop contribution, i.e. for the tree-level result, 1 for the one-loop result includ- 
ing the SUSY-QCD corrections, and A for the Am{,-corrected tree-level result. 
The superscript indicates the renormalization scheme in which the quantity is 
calculated. 



6.3 Numerical Results 

In this section the numerical results for the production process bb -^ H^W~ are 
presented. All quoted cross sections are given for the LHC with a proton-proton 
center-of-mass energy of 14 TeV. First of all we investigate the effect of varying 
the MSSM parameters on the SUSY-QCD corrections. To do so a parameter 
point is chosen for which the corrections of the Arrtb term are expected to have 
a large impact. To that effect the parameter point 

mH+ = 200 GeV 
t/3 = 30 

Mq = M(j = Mf) = -fx = m-g e [0, 10] TeV (6.7) 

is used as input, t^ of this point is fairly large to enhance the A?Tif, contribution. 
The soft SUSY-breaking mass terms in the squark sector, jointly denoted as 
-^susY = Mq = Mfj = Mfj, — /i and the gluino mass rrig all take the same value 
which is varied over a large mass range. Eq. 4.17 then predicts that the SUSY- 
QCD corrections in the OS scheme should be large and independent of the varied 
mass scale as we are in the limit where all these masses are equal and the mass 
dependence drops out. This is indeed the case as one can see in Fig. 6.4. A naive 
calculation in the OS scheme yields a correction of above 60%. After subtracting 
the Ami) contribution only the real one- loop corrections are left. Their size is of 
around 0.2%. In the DR scheme the corrections are equally small for small mass 
values and show a logarithmic rise with growing mass. This is an artefact of the 
mismatch between the renormalization scale and the masses of the squarks and 
gluinos appearing in the loop diagrams. The former one was fixed to the sum 
of the final-state particle masses hr = mu- + raw- Terms of the order In ^^^^ 
appear in the expression which originate from the dimensional regularization of 
the divergent loop integrals. Thus this logarithmic rise bears no physical meaning 
and will vanish if higher orders of perturbation theory are taken into account. 
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Figure 6.4: Hadronic cross section differences for the process bb -^ H~^W~ in the 
large tp regime. The soft SUSY-breaking mass terms in the squark sector Msusy, 
— /i and the gluino mass rrig are fixed to the same value and varied over a large 
mass range. The parameter set eq. (6.7) was used to obtain this plot. 
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Figure 6.5: Hadronic cross section differences for the process bb -^ H^W~ using 
the parameter set eq. (6.8). Only n is varied in this plot. 



In the next plot (Fig. 6.5) the squark and gluino masses are now fixed and 
only /i is left as a free parameter. The parameter set for this plot is 



mH+ = 200 GeV 



A, 



30 

At 



Mq = Mjj = M^ = 500 GeV 
m-g = 580 GeV . 



(6. 



As expected the cross section in the OS scheme grows linearly with /i. For fi 
values less than about —750 GeV the corrections even exceed 100%. Again, 
when the Arrib corrections are subtracted and only the true one-loop SUSY-QCD 
corrections remain the order of the corrections is below 1% and shows only a very 
low variation with /i. In the DR scheme the corrections are also much smaller 
than in the OS scheme and almost constant as a function of /i, as is expected 
from the remaining corrections appearing in this scheme. 

As next step the dependence of the hadronic cross section differences as a 
function of t/3 is investigated in Fig. 6.6. Here a data point with a smaller soft- 
supersymmetry breaking mass is chosen to emphasize the effect which will be 
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Figure 6.6: Hadronic cross section differences for the process bb -^ H~^W 
function of t^. This plot was calculated with the parameter set eq. (6.9). 



as a 



discussed below, namely 



mH+ = 200 GeV 
/i = -200 GeV 



Ma = M, 



M 



u 
m-g = 580 GeV 



D 



200 GeV 



(6.9) 



For large values of tp the respective corrections feature the behavior which was 
already observed in the previous plots. The Arrif, corrections are large, while the 
true one-loop corrections almost vanish. The corrections in the OS scheme rise 
linearly with tp as is expected from eq. (4.17). For t^ < 15 however the behavior 
changes. There the difference between the full one-loop computation and the 
Am^-corrected tree-level cross section can increase up to 10%. This contribution 
for small tp originates mainly from the diagram given in Fig. 6.7. The Yukawa 
coupling of the charged Higgs to the stop and sbottom is proportional to the 
top-quark mass if a right-handed stop couples to a left-handed sbottom. Another 
factor rrit appears in the trace over the fermion line where for the internal top 
quark line the mass term in the Dirac algebra is chosen. These two factors 
cancel the top-quark propagator which is dominated by the mass term and the 
top-quark mass dependence drops out. To get this left-right mixing term in the 
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Figure 6.7: Dominant one- loop contribution in the case of small tp. The crosses 
in this diagram denote a mass insertion, i.e. when calculating the trace over the 
fermion line the mass term in the Dirac algebra is chosen. 



Yukawa coupling also the mass term for the gluino appears during the calculation 
of the fermion trace, giving a factor rrig in the nominator. Accordingly this vertex 
correction to the tbH~^ vertex is proportional to 

-^fimgl{mi,mi,'mg) (6.10) 

where / was given in eq. (4.18) and is related to the three-point integral in the 
limit of vanishing external momenta. The expression is independent of tp. So for 
small values oft/3, where the bottom-quark Yukawa coupling is not enhanced, this 
gives an important contribution. When t^ takes larger values, the bottom-quark 
terms dominate. These terms have a factor of t^ appearing in the amplitude so 
the cross section increases quadratically with t'i. Hence the relative contribution 
of Fig. 6.7 is reduced. This plot underlines the importance of performing a 
full one-loop calculation because only in such a full calculation the non-universal 
corrections are included and a tree-level calculation with effective couplings would 
give wrong results in the low-t/3 regime. 

The variation of the cross section as a function of fi in the low-t^ regime is 
investigated in Fig. 6.8. For this plot t^ = 6 was chosen and the other parameters 
were left unchanged from the last plot. Again the rather small value of 200 GeV is 
chosen for the soft-supersymmetry breaking masses in the squark sector, so that 
the function /, which is proportional to -^ — , has a small denominator. One 
can clearly see two distinct properties. The one-loop corrections which cannot be 
absorbed into a redefinition of the bottom-quark Yukawa coupling now give a sig- 
nificant contribution. The Arrih corrections still yield an important contribution 
as can be seen from comparing the Aqs and AAm^ curves. Yet after subtracting 
the universal corrections to the bottom-quark Yukawa coupling the remaining 
non-universal corrections are large. For large values of /i they can reach more 
than 50%. They are equally present in the DR scheme where a numerical contri- 
bution close to the A?Tife-corrected result is obtained. Furthermore the correction 
to the cross section increases approximately linearly with the absolute value of fi 
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Figure 6.8: Hadronic cross section differences for the process bb -^ H~^W~ as a 
function of fi. For all other parameters the parameter set eq. (6.9) with tp = 6 
was used. 



as is expected from eq. (6.10). 

Furthermore, in Fig. 6.9 a plot where the complete SUSY mass spectrum, 
i.e. soft-supersymmetry breaking mass terms for the squarks, fi and the gluino 
mass, is fixed to the same value and run up to 10 TeV is presented for the low-t/j 
regime. The behavior as a function of the mass scale is the same as before in the 
case of large tf^. After having subtracted the Arrih corrections from the complete 
one-loop contributions there is still a correction of the order of 10% left which 
mainly originates from the diagram given in Fig. 6.7. For small mass values the 
relative correction slightly drops because in this region other diagrams also give 
a numerically significant contribution. 

The scale dependence of the SUSY-QCD corrections is given in Fig. 6.10. 
The factorization and renormalization scale of the process are fixed to the same 
value and varied between 0.1 and 10 times their basic value /i_R = /xf = 'm^w + 
rriH which is used for the other plots. Even for this large scale variations there 
is only a mild dependence for the corrections in the OS scheme and the Am^ 
corrections. This is due to the fact that the only scale-dependent parameter is 
the strong coupling constant ag and the PDFs. On-shell conditions render all 
other paramaters independent of the scale. In the DR scheme also the quark 
masses are scale dependent resulting in a much larger variation as a function of 
the scale. Including the Standard-QCD corrections which are already known for 
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Figure 6.9: Hadronic cross section differences for the process bb -^ H^W~ in the 
low-t/3 regime, i.e. t/5 = 6 is used. All other parameters take the values given in 
eq. (6.7). The soft SUSY-breaking mass terms in the squark sector, — // and the 
gluino mass are fixed to the same value and varied over a large mass range. 




|iR/(mw+mH) 

Figure 6.10: Hadronic cross section differences for the process hb -^ H^W~ as a 
function of the renormalization and factorization scale. The plot was calculated 
using the parameter set eq. (6.9) with additionally t/3 = 6. 
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Figure 6.11: Total hadronic cross sections for the process hh -^ H^W as a 
function of t/j using the parameter set eq. (6.9). 



this process [104] would reduce this variation, but implementing these additional 
contributions was beyond the scope of this dissertation. 

The last plot in Fig. 6.11 shows the total hadronic cross section as a function 
of t/3 in both OS and DR renormalization schemes. Various cross sections with 
different contributions taken into account are presented here. The same param- 
eter set as in eq. (6.9) is used for this plot. In all cases the total cross section 
rises quadratically with t^ in the region where t/3 is larger than about 15. This 
is the parameter space where the Yukawa coupling to the charged Higgs boson is 
dominated by the term proportional to the bottom-quark mass, which scales with 



AS 


Tree-level cross section 


UQ 


OS 


Am^-corrected tree-level result 


(Ta 


OS 


One-loop OS cross section 


ai 


OS 


One-loop OS cross section including 

resummed higher-order Amb corrections 


O-fuU 


DR 


One-loop DR cross section 


(Ti ^ 



Table 6.1: Key to the total hadronic cross sections of the process bb 
plotted in Fig. 6.11. 
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tp, and gives the leading contribution to the cross section. On the left-hand side 
of the plot, where t^ is small, in contrast the top-quark mass part is responsible 
for the overall behavior of the cross section and leads to a decrease with tf^. In the 
intermediate region both terms contribute equally much, leading to a minimum 
of the cross section for tp ^ 8. 

In total five different cross-section types including miscellaneous contributions 
are depicted. An overview is given in table 6.1. ao*"*^, the straight red line, denotes 
the tree-level contribution in the OS renormalization scheme. The short-dashed 
blue line, cq^^, is the one-loop cross section in the OS scheme without having 
used any further effective couplings, ctq^^, the long-dashed green line contains 
the Amfe-corrected tree-level result. As seen before in the plots of the relative 
corrections, this line agrees with the complete one-loop result in the case of large 
tjS as in this region only the universal Ami, corrections are relevant. In the 
small-t/3 regime these terms can only account for a part of the total corrections. 
There are also non-universal terms, mainly coming from the sub-diagram given 
in Fig. 6.10, which cannot be absorbed into an effective coupling. These ones do 
not, as observed before, contain any factors of tp and hence their effects diminish 
for higher t^ values as the total cross section scales with t^. For cTfuu^^, the dotted 
pink line, the Arrib corrections are included in the resummed version of eq. (4.16). 
Additionally the non-universal one-loop corrections are added. To avoid double- 
counting the non-resummed Arrih contribution must then be subtracted again, 
resulting in the following formula 

In this line also phase-space effects from the reduced bottom-quark mass are taken 
into account, leading to an additional shift with respect to the tree- level cross 
section. Nevertheless the curve again shows the expected behavior which can be 
deduced from the results given above with rising corrections for increasing tp. 
For large values of tf^ they largely exceed the one-loop result because of the Arub 
resummation where contributions are added that in a naive calculation would 
appear in higher-order diagrams of perturbation theory. This curve presents 
the current best estimate for the total cross section in the OS scheme where 
higher-order contributions are included as much as possible. Lastly the one- 
loop expression in the DR scheme is plotted as dash-dotted light-blue line and 
labeled cxi^^. It has a shape very similar to the previous curve because the Ami, 
corrections appear in the self-energy contributions to the bottom-quark mass. 
In this renormalization scheme they enter completely at one-loop order via the 
bottom-quark propagators and no further contributions at higher orders appear. 
Hence this corresponds to the resummed result in the OS scheme. The remaining 
difference between the two curves is a measure for the theoretical uncertainty of 
the calculation because of missing non-leading higher-order contributions from 
perturbation theory. 
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Tree-level cross section 


ao =2.684 fb 


Amfc-corrected tree-level cross section 


(TAm^ = 2.266 fb 


One-loop OS cross section 


(Ti =2.176 fb 


Relative one-loop OS correction 


Aos =-15.6 % 


Relative true one-loop OS correction 


AAm,= -4.0 % 



Table 6.2: Hadronic cross sections for the reference point SPSla', which is de- 
scribed in appendix A. 2. 



Finally the numerical results for the parameter point SPSla' [16] are given in 
table 6.2. This parameter point is described in appendix A. 2, it was chosen as a 
reference point for MSSM calculations. Because of the positive sign of // the one- 
loop cross section is now reduced with respect to the tree- level result. The tp value 
of 10 is in a region where the Arub corrections are already the dominant ones, but 
the non-universal corrections still yield a numerically significant contribution. 



Chapter 7 

Higgs-Boson Production via 
Vector Boson Fusion 



Proving the existence of a neutral Higgs boson is one of the main tasks of the 
LHC. Its main production processes for both SM and MSSM Higgs bosons include 
Higgs-boson production via vector-boson fusion (Fig. 7.1) [107, 108, 109]. Its rate 
is surpassed only by the gluon-fusion process gg —^ h^ [HO] shown in Fig. 7.2. 
This process has large NLO-QCD corrections with K-factors larger than 2 [111]. 
Even after including the NNLO-QCD corrections [112] theoretical uncertainties 
of O (10 — 20%) remain. They make the extraction of coupling constants from 
the gluon-fusion process difficult and lead to large errors. 

The Standard-QCD corrections to the vector-boson-fusion process were al- 
ready calculated before [113, 114]. At tree-level the process only consists of a 
t-channel exchange of a colorless gauge boson, which is why the contribution 
where an additional gluon connects the two quark lines has no interference term 
with the tree-level diagram. Only the quark-anti-quark-vector-boson vertex re- 
ceives corrections of O (a^) and hence the overall QCD corrections are relatively 
small and typically between 5 and 10%. Therefore, the extraction of Higgs cou- 
pling constants leads to much smaller errors than in the gluon-fusion case, and 
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Figure 7.1: Generic tree- level Feynman diagram of the vector-boson- fusion pro- 
cess qq -^ qqhP 
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Figure 7.2: Leading-order Feynman diagrams of the gluon-fusion process gg — > h^ 



the vector-boson-fusion production process is, despite its lower cross section, an 
ideal instrument to study the Higgs boson. 

This process possesses a clear experimental signature of two jets in the forward 
region and thus can be easily separated from background processes by applying 
appropriate cuts [115]. In this chapter the production of the lightest MSSM Higgs 
boson h^ via vector-boson fusion is investigated and the SUSY-QCD corrections 
to this process are calculated. 



7.1 The Partonic Process 



The partonic processes which contribute to the production of a Higgs boson via 
vector-boson fusion can be summarized in a single general Feynman diagram 
which is depicted in Fig. 7.1. It can be seen as scattering of two quarks which is 
mediated via a vector boson in the t-channel with a Higgs boson being radiated off 
the intermediate vector boson. This is why the process has a clear experimental 
signature of two jets in the forward region of the detector which allows one to 
easily distinguish the signal from background processes by using appropriate cuts. 
In a strict sense, the general diagram given in Fig. 7.1 is not the only one 
which contributes to this final state. When a quark-anti-quark pair of the same 
flavor appears in the initial state they can form a Zh^ pair via an intermediate 
virtual Z boson and the Z subsequently decays again into a quark-anti-quark 
pair. Hence these diagrams have exactly the same particle content in both the 
external and internal lines. There are also similar processes where an interme- 
diate W boson can appear in the same way. However, the vector-boson- fusion 
process has a very distinct signature of two jets in the forward region. Using 
only this particular phase-space region the interference between the two diagram 
types is strongly suppressed by the large momentum transfer appearing in the 
intermediate gauge bosons. Additionally a color suppression factor appears in 
the interference term [114]. Hence the effects from these additional diagrams can 
be safely neglected if appropriate cuts [115] to the phase space are applied. 
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An analysis of the statistical accuracies of the cross section which are achiev- 
able at the LHC was done in refs. [116, 117]. It could be shown that a measure- 
ment with an accuracy of 5 to 10% can be performed, also taking uncertainties 
in the decay branching ratios of the Higgs boson into account. So the size of the 
NLO-QCD corrections already matches the accuracy which is achievable in ex- 
periment. Theoretical uncertainties are not expected to limit the precision with 
which the cross sections can be measured. 

In the MSSM besides the Standard-QCD corrections also SUSY-QCD cor- 
rections appear which are of the same order O (as) or, in case of the pentagon 

diagrams, even enhanced and of (9 ( ^ 1 in the coupling constant. A full one- loop 

calculation of these corrections has not been done before and is presented in this 
thesis. 

7.2 SUSY-QCD Corrections 

In this chapter the SUSY-QCD corrections to /i°-production via vector-boson 
fusion are studied. If their size is larger than the experimental uncertainties one 
might be able to use this to tell the SM and the MSSM apart. In the limit of large 
tua the couplings of the h^ become SM-like. So if at the LHC only one Higgs 
boson with Standard-Model couplings and a mass below 140 GeV is found, the 
question arises whether a SM or a MSSM Higgs boson was seen in the detector. 
The SUSY-QCD corrections, which exist only in the case of a MSSM Higgs boson, 
could modify the Higgs boson coupling by an amount large enough and therefore 
make the distinction between the two models possible. Also if supersymmetry 
could be established by other means beforehand, these corrections give an indirect 
contribution to the coupling between the Higgs boson and two gauge bosons. To 
be able to extract the value from the experiment as precisely as possible it is 
necessary to include these higher-order corrections. 

The possible types of diagrams which appear in the SUSY-QCD corrections 
are depicted in Fig. 7.3. The quark-quark-gauge boson vertices receive cor- 
rections which are depicted in Fig. 7.3(a). Their divergencies are cancelled by 
appropriate counter terms which are shown in Fig. 7.3(b). The contribution of 
these diagrams was already investigated before [118] for the special case where all 
squarks have equal mass. Additionally one of the gauge bosons can be replaced 
by a box-shaped sparticle loop as shown in Fig. 7.3(c). Finally, all external par- 
ticles can be coupled via a pentagon-type loop as in Fig. 7.3(d). Because of the 
Majorana nature of the gluinos also the diagram displayed on the right-hand side 
of the figure, where two quark lines are connected, exists. 

All cross sections are calculated in the OS renormalization scheme. For the 
tree-level diagrams the leading-order parton distribution functions of ref. [119] 
were used. The one-loop cross sections were convoluted with the NLO-PDFs of 
the same group given in ref. [120]. In both cases the implementation from the 
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Figure 7.3: Diagram types contributing to SUSY-QCD corrections to h^ produc- 
tion via vector-boson fusion 
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program package LHAPDF [93] was used to obtain the numerical results which 
are presented in the next section. 

7.3 Numerical Results 

In this section numerical results for the SUSY-QCD corrections to /I'^-production 
via vector-boson fusion are presented. The hadronic cross sections for each indi- 
vidual process for the MSSM reference point SPSla', which is described in detail 
in appendix A. 2, are listed in the following table. The one- loop corrections are 
separated according to the loop type, where a^crtcx includes the contributions 
from the diagrams shown in Fig. 7.3(a) and (b), o"box from the ones of Fig. 7.3(c) 
and CTfive-pt. from the pentagon diagrams depicted in Fig. 7.3(d). 
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The quoted cross sections are hadronic ones, where the convolution with the 
PDFs has aheady been performed. They are given separately for each partonic 
subprocess to facihtate an easier analysis of the characteristic effects appearing 
in this process. Additionally the total hadronic cross section is stated, which is 
the sum of all subprocesses. For some partonic subprocesses no five-point loop 
diagrams exist at all. In this case the entry in the last column of the table is 
exactly zero. To exploit the unique characteristics of the final state of this process 
cuts were used to obtain the cross sections. A lower limit was placed on the 
transverse momentum pt and the pseudo-rapidity r] of the outgoing quarks and 
anti-quarks, so that the final-state jets are clearly separated from the beam pipe, 
but still in the forward region of the detector. Also a cut on the jet separation AR 
between each combination of outgoing particles was set to emulate the behavior 
of the jet-clustering algorithms used in experimental analyses and to be able to 
resolve the particles in the detector separately. Thus the apphed cuts were 



PT(g,g) > 40 GeV r]{q,q)>2 



AR 



m,qq,q<i,<iho,qho 



> 0.4 



(7.1^ 



The formal definition of these quantities was given in chapter 5.4. 

There is an interesting observation already at tree-level. The partonic sub- 
processes which enter via a Z-boson exchange are suppressed with respect to the 
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ones with a W boson as intermediate vector boson. Firstly, the couphng of the 
W boson to the quarks is enhanced by a factor of — , the inverse of the cosine 
of the electroweak mixing angle, with respect to the leading term of the qqZ- 
coupling. Secondly, the Z boson is heavier than the W boson, and the ratio of 
the two masses is also equal to ^. As the gauge-boson propagators are domi- 
nated by their mass terms this leads to an additional enhancement of -5- for each 
W boson propagator. So in total the amplitude oi aW boson-exchange diagram 
is enhanced by about -|- over one with a Z boson exchange. Accordingly, this 



1 _ 1 



amounts to a factor -^^ = ,„ „^^,ia — 4.8 for the tree-level cross section, which 
corresponds to the observed partonic cross-section ratios. Because of this effect 
and the large valence-quark densities of the up- and down-quarks in the proton, 
the partonic subprocess ud —>■ udh^ gives the leading contribution to the hadronic 
process. 

The vertex corrections all have the same size relative to the respective tree- 
level cross section and correspond to the results obtained in ref. [118]. Since the 
coupling of the W boson to the quarks is purely left-handed, they are largest if 
the off-diagonal elements in the squark mixing matrix are small and therefore left- 
and right-handed squarks have almost equal masses. Also for the intermediate 
gluino propagator only the momentum term survives when calculating the trace 
over the fermion line. It is proportional to the momentum transfer in the t- 
channel and thus small. For diagrams with Z-boson exchange the situation is 
more complicated, because the ggZ-coupling contains both left- and right-handed 
parts. Nevertheless, any subdiagrams involving a mixing of left- and right-handed 
squarks are proportional to the off-diagonal terms in the squark mixing matrix, 
which contain the Yukawa coupling of the corresponding quark and are hence 
small. Hence, the same effects as in the W boson case appear and lead to a 
similar relative size of the vertex corrections. 

In the case of box diagrams the relative size of the corrections shows a much 
wider range. For some partonic subprocesses they exceed the vertex corrections 
significantly, which underlines the importance of performing a full one-loop calcu- 
lation to take all effects into account. Yet for other subprocesses they are much 
smaller. This is due to the fact that for l^-boson exchange large cancellation 
is manifest. It occurs between the diagrams on the left-hand side of each row 
of Fig. 7.3(c), where the h^ couples to the squark with the flavor of the out- 
going quark, and the ones on the right-hand side which have an /I'^-coupling to 
the squark with the incoming-quark flavor. Because of the W boson one of the 
squarks is always up-type and the other one down-type. Their couplings to the 
h^ have a minus sign relative to each other. Any effects from CKM-mixing are 
neglected, therefore only squarks of the same generation can appear in a single 
diagram. This is also why only the superpartners of the four light quarks con- 
tribute at all. They all have very similar masses, so the absolute value of the 
ggl^-coupling is roughly the same everywhere. Hence the diagrams on the left- 
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and right-hand side of Fig. 7.3(c) almost exactly cancel for VT-boson exchange, 
leading to a strong suppression of this contribution. For Z-boson exchange no 
such effect occurs and the relative box-diagram contribution is larger by an or- 
der of magnitude. Since, as mentioned above, already the tree-level amplitude 
is smaller for this diagram type, the absolute value of the correction is small as 
well and these corrections cannot give a significant contribution to the total cross 
section. Again the correction is maximized in the case where the off-diagonal 
elements in the squarx mixing matrix are small. Also the change of the coupling 
constant between tree-level, which is proportional to Sa+/3, to the one- loop one 
of SjS-a cannot give an important effect, because the ratio of the two is always 
above 0.9 and approaches 1 in the decoupling scenario, where the additional Higgs 
bosons of the MSSM are heavy and the /i°-coupling becomes SM-like. 

For the five-point diagrams a cancellation similar to the box diagrams occurs 
when the corresponding tree-level process is mediated by a IV-boson exchange. 
Additionally at least twice a left-right mixing term in the squark sector appears 
in the amplitude. For this reason a term proportional to the Yukawa coupling 
of the four light quarks enters the expression and leads to a suppression of the 
one-loop correction. The choice of parameters which yield the biggest one-loop 
corrections are in this case large terms in the off-diagonal entries of the squark 
mixing matrices and therefore a larger mass splitting in the squark sector. Ad- 
ditionally for both types the larger masses of the gluinos and squarks, where 
experimental limits require that they are heavier than the W or Z boson, lead 
to a further reduction of the cross section. Yet, there is also an enhancement 
factor. Except for the Higgs coupling, all other four couplings of the pentagon 
diagrams are proportional to the strong coupling constant, while for the trian- 
gular and box-type diagrams two of the couplings are strong and two of them 
are electroweak. Hence this type of diagram contains an enhancement factor of 
— ~ 14. This is however not sufficient to give a significant contribution to the 
cross section. 



7.4 /i^-Production with External Gluons 

Additionally the production of an h^ with one or two gluons in the initial state was 
considered. The corresponding Feynman diagrams are shown in Fig. 7.4. Since 
there is no tree-level coupling of the gluons to the Higgs bosons, these processes 
occur at the one-loop level in leading order. This leads to an additional factor 
of as in the cross section so the total amplitude is of O (af a) , while the vector- 
boson- fusion diagram has a factor a^. Given that a^ is of the same numerical size 
as a these contributions might prove important. Additionally, the gluon densities 
in the proton are much higher than those of the quarks at typical LHC energies 
and this further enhances this type of diagram. 

The numerical results for these processes with external gluons for the reference 
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(b) box diagrams 
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Figure 7.4: Leading-order diagram types for /i°-production with one external 
gluon in the initial and final state. The Feynman diagrams with two gluons in 
the initial state have the same topology. They are obtained from these ones by 
taking both gluons as incoming particles and changing the incoming quark to an 
outgoing anti-quark. 
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Table 7.2: One- loop hadronic cross sections for the subprocesses with one or two 
gluons in the initial state for the MSSM reference point SPSla'. 



point SPSla' are given in table 7.2. In this case the convolution with the PDFs 
is also already included in the numbers for the cross section. The same cuts as 
for the vector-boson-fusion process, given in eq. (7.1), were used. Also the same 
cuts were applied to the final-state gluons as to the quarks and anti-quarks. 

Formally this process type constitutes a background to the previously con- 
sidered process of /i°-pro duct ion via vector-boson fusion. Therefore one wants to 
pursue the question, how large the total contribution of these diagrams is and if 
there are cuts which reduce their size with respect to the signal process. 

The processes with one gluon in the initial and final state have a momentum 
distribution similar to the vector-boson- fusion one. The gluon densities in the 
proton are very large for small x, but rapidly diminish for larger x. In contrast 
the sea-quark densities fall off much slower and the valence-quark densities have 
their maximum at about |. So the most favorable configuration is the one where 
the energy to produce the final state mostly originates from the quark. Thus the 
hadronic center-of-mass frame is strongly boosted with respect to the partonic 
one, which leads to jets in the forward region of the detector. This would mimic 
the signature of a vector-boson-fusion process and produce events which cannot 
be eliminated by cuts. In contrast for processes with two gluons in the initial 
state the momentum configuration which maximizes the hadronic cross section 
is the one where both gluons have similar values of x. This leads to more central 
jets which are suppressed by the applied cuts. 

As one can see from the cross sections in table 7.2, the total contribution of 
these diagrams is of O (10-^) and therefore well below the experimental uncer- 
tainties, which are in the range of 5 to 10%. As the total contribution of these 
background processes is below the statistical uncertainties which can be reached 
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in a measurement of the vector-boson-fusion cross section, these background pro- 
cesses do not affect the experimental determination of the /i°-production rate via 
vector-boson fusion. 



Chapter 8 

Higgs-Boson Production in 
Association with Heavy Quarks 



The coupling of Higgs bosons to fermions is of the Yukawa type and therefore 
proportional to the mass of the fermion. The four light quarks, u, d, s and 
c, all have a mass below or of about 1 GeV. This mass should be compared 
to the Higgs vacuum expectation value v, the scale of electroweak symmetry 
breaking, to obtain the strength of their respective Yukawa couplings, which 
are therefore small. In contrast the top-quark mass is of the same order as the 
electroweak symmetry-breaking scale, making the top-Higgs coupling numerically 
sizable. The bottom-quark mass of a few GeV also leads to a rather weak coupling 
to the Higgs boson in the Standard Model. In the MSSM, the coupling to the h^ is 
enhanced by a factor tp, so for large values of t/j its size can become comparable to 
the top-quark Yukawa coupling. These large Yukawa couplings make Higgs-boson 
production in association with heavy quarks [121, 122] a phenomenologically 
interesting process. 

In this chapter the production of the lightest CP-even neutral MSSM-Higgs 
boson h^ in association with a bottom or top quark-anti-quark pair is studied. 
The top quarks decay rapidly into mainly bW and the outgoing bottom quarks 
can be identified in the detector via 6-quark tagging. Therefore both processes 
form distinct final states which neither interfere with each other nor with the 
Higgs-boson production via vector-boson fusion presented in chapter 7. First the 
peculiarities of each of the two processes are discussed separately. Then the one- 
loop SUSY-QCD corrections for both processes are described. Since the same 
basic Feynman diagrams appear in both cases, this task is done jointly. The 
Standard-QCD corrections have already been calculated in refs. [123, 124] for 
66/i°-production and refs. [125, 126, 127] for tt/i°-production. Finally in the last 
two sections the numerical results for both processes are shown. 
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(a) s-channcl contribution 
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(b) t-channcl contribution 

Figure 8.1: Tree- level Feynman diagrams of the process pp 



bbh^ 



8.1 The bbh^ Final State 

The production of a Higgs boson in association with bottom quarks in the Stan- 
dard Model was intensively studied in the literature [117, 128, 129]. At tree-level 
it originates from the annihilation of a quark-anti-quark pair or from a gluon fu- 
sion process, where the final-state 66-pair is produced via an intermediate gluon, 
and the Higgs boson radiates off from one of the bottom quarks (Fig. 8.1(a)). 
Besides these s-channel diagrams the partonic gluon-fusion process also proceeds 
via a t-channel diagram shown in Fig. 8.1(b), where the Higgs boson can be 
emitted from both the internal and external bottom-quark lines. The analysis 
was soon extended [83, 94, 130] to include the lightest MSSM-Higgs boson h^. 
The diagram types are exactly the same as in the Standard Model case. Only 
the bottom-quark-Higgs coupling is changed to its supersymmetric counterpart, 
resulting in 



ctmssm {pp -^ bbh^) 



-— osm ipp^ bbH) 



l.l\ 



where ——is the ratio of the bottom-quark coupling to the MSSM h^ boson and 
to the one of the SM Higgs boson H . 

The Standard-QCD corrections [123, 124] to this process are also known and 
reduce the dependence of the cross section on the factorization and renormaliza- 
tion scales. 

However, there are subtleties when making a theoretical prediction for total 
integrated /i°-production via this process, i.e. when the final-state bottom quarks 
are not explicitly detected. In a four-flavor-number scheme, where only gluons 
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and the four light quarks, but no bottom quarks appear in the initial state, large 
logarithms of (9 ( In ^ j in the total cross section emerge, where Q is of the order 
of the Higgs-boson mass. They arise from the kinematical configuration where a 
gluon splits into a 66-pair and the bottom quarks are collinear to the gluon. These 
logarithms can be resummed using bottom-quark parton densities, thereby using 
a five-flavor-number scheme. The bottom-quark densities in the proton originate 
purely from such splitting gluons. So for every bottom quark which appears 
in a partonic process another bottom (anti-) quark exists in the hadronic final 
state. This scheme uses the approximation that the outgoing bottom quarks have 
small transverse momentum and they are assigned zero transverse momentum at 
leading order. In the five-fiavor-number scheme the leading-order partonic process 
is then 66 — > /i°. gg ^ hhh^ only appears at NNLO together with the two- loop 
corrections to this process [131]. 

In our case though these large logarithms are avoided by requiring bottom- 
quark jets with high transverse momenta and a tagging of the final-state bottom 
quarks in the detector. The additional cuts reduce the cross section by one or two 
orders of magnitude, but also greatly reduce the background and make this ap- 
proach more interesting. The existence of bottom-quark jets with large transverse 
momenta also guarantees that the Higgs boson was emitted from a bottom quark 
and is therefore proportional to the bottom-quark Yukawa coupling, allowing its 
precise measurement. 

The SUSY-QCD corrections to this process were partly calculated in ref. [132]. 
There an effective 66/i°-coupling was used which includes the one-loop squark 
and gluino contributions, but no box-type or pentagon diagrams were added in 
their analysis. In this dissertation a full one-loop calculation of the SUSY-QCD 
corrections is performed. 



8.2 The tth^ Final State 

The production of a Higgs boson in association with a top quark-anti-quark 
pair [122, 133, 134] proceeds in the same way as the one with a bottom quark- 
anti-quark pair discussed in the previous section and the same diagrams as in 
Fig. 8.1 appear. Since the mass of the top quark is of the same order as the 
electroweak symmetry-breaking scale, its Yukawa coupling gives a sizable contri- 
bution and the process is an important channel for Higgs-boson production in 
the mass region below 125 GeV [135]. Furthermore, this process can be used to 
measure the top-quark Yukawa coupling precisely [136]. The extension of the SM 
tree-level calculations to the MSSM, where the Higgs boson is an /i°, is again 
straightforward and amounts to a replacement of the Yukawa coupling such that 

ctmssm [pp -^ tth^) = f — ) cTsM {pp -^ tiH) , (8.2) 
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where — is the ratio of the top-quark couphng to the MSSM h^ boson and to the 
Higgs boson H of the SM. Thus the total cross section in the MSSM is reduced 
with respect to the SM one by approximately a factor of tt- 

The Standard-QCD corrections for this process are available in the litera- 
ture [125, 126, 127]. Their numerical size is of O (20% - 40%) and leads to a 
stable prediction of total and differential cross sections with respect to variation 
of the renormalization and factorization scales. 

The large mass of the top quark also reduces the size of the collinear loga- 
rithms to (9 (in ^j. Now the argument of the logarithm is close to 1. So the 

higher-order corrections are small and no resummation of these terms needs to 
be performed. Hence one can safely use the four-flavor-number scheme for this 
process and need not apply any additional cuts to the final-state top quarks in 
this case. 

A calculation of the SUSY-QCD corrections was performed quite recently in 
ref. [137]. As the figures of this article include both the Standard-QCD and the 
SUSY-QCD one-loop contributions a direct comparison of the numerical results 
is difficult. As far as the principal behavior with respect to a variation of the 
MSSM parameters is concerned, agreement could be found. 



8.3 SUSY-QCD Corrections 

In this section the SUSY-QCD corrections to /i°-production in association with a 
heavy quark-anti-quark pair are described. In the Feynman diagrams the heavy 
quark is denoted by a Q, which represents a b for the bottom-quark and a t for 
the top-quark final state. Correspondingly, in part of the diagrams the supersym- 
metric partners to the heavy quark appears, which is marked by Q, specifying 
b and t, respectively. A small q on the other hand signifies that all squarks can 
be inserted in the propagator. The q' in the qq diagrams denotes the super- 
partner to the initial-state quark. In Fig. 8.2(a)-(f) the basic types of Feynman 
diagrams which appear as one-loop SUSY-QCD corrections to /i°-production via 
gluon fusion are depicted. Self-energy corrections (Fig. 8.2(a)) enter either via 
a squark or gluino loop which is inserted into the intermediate gluon propaga- 
tor, or a combined squark-gluino two-point loop inserted into the heavy-quark 
line. In Fig. 8.2(b) the possible vertex corrections are displayed. Squark loops 
induce an effective gluon-Higgs coupling appearing in the two diagrams on the 
left-hand side of the first two rows. The diagrams on the right-hand side of the 
first two rows contain a correction to the coupling between h^ and the heavy 
quark. The third and fourth row feature corrections to the gluon-quark inter- 
action and the last row an additional contribution to the triple-gluon vertex. 
Diagrams where four particles are connected via a sparticle loop are presented 
in Fig. 8.2(c). Finally, all five external particles can be joined by a squark- 
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(a) self-energy corrections 
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Figure 8.2: Types of Feynman diagrams contributing to SUSY-QCD corrections 
to /i^-production in association with heavy quarks 
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(c) box diagrams 
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(c) self-energy counter terms 
Figure 8.2: (continued) 
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(f) vertex counter terms 
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(g) additional diagrams for quark-anti-quark annihilation 

Figure 8.2: (continued) 
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gluino loop as shown in Fig. 8.2(d). The emerging divergences are cancelled by 
counter-term diagrams shown in Fig. 8.2(e) for the self-energy contributions and 
(f) for the vertex corrections. All gluon-fusion s-channel diagrams, where the 
two initial-state gluons couple to a further, intermediate gluon, also exist in the 
quark-anti-quark-annihilation subprocesses. The only change is the replacement 
of the two incoming gluons by a quark-anti-quark pair. The additional diagrams 
which appear for this type of subprocess are depicted in Fig. 8.2(g). They are 
corrections to the quark-anti-quark-gluon coupling together with the associated 
counter-term diagram as shown in the first row. Secondly, the coupling of the 
incoming light quark to the Higgs boson, which is neglected at tree-level, appears 
at one-loop order as indicated by the diagram in the second row. Finally in the 
last row the additional box and pentagon diagrams are shown. 

In the remaining sections of this chapter the numerical results of /i°-production 
in association with a heavy quark-ant i- quark pair are presented. In analogy to 
chapter 6 several cross-section differences are defined to illustrate the results. All 
calculations in this chapter were performed in the OS renormalization scheme, so 
the label indicating the renormalization scheme will in the following be dropped 
for all items. The relative one-loop correction is defined as 

A, = ^i^^ , (8.3) 

where ctq denotes the tree-level and cri the one-loop cross section. For the calcu- 
lation of hadronic cross sections the PDF set of ref. [120] was used. Additionally, 
a A?7ife^i-corrected tree-level cross section o"a was calculated in a similar way as 
already described in chapter 6, i.e. using the non-resummed version of eq. (4.21) 
and treating the Anib^t term as a one-loop contribution. Additionally, the contri- 
bution to the vertex from the term proportional to the second mixing angle in the 
MSSM-Higgs sector, a, was included in a a according to eqs. (4.27) and (4.33). 
The relative correction using only these contributions is defined as 

A, = "-^^ . (8,4) 

Finally, the difference between the Arrib, t-corrected tree-level cross section and 
the full one-loop result, which denotes the true one- loop corrections, is given by 

A.„,. = ^^i^i^ . (8.5) 

The renormalization of the strong coupling constant a^ was performed as de- 
scribed in chapter 4.2.3. 

8.4 Numerical Results for bbh^ 

In this section the numerical results for the process pp — > bbh^ are presented. First 
the total hadronic cross section for the MSSM reference point SPSla' is given in 
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Partonic subprocess 


(To [fb] 


(Ti [fb] 


Ai [%] 


Ag [%] 


dd -^ bbh^ 


0.107 


0.104 


-2.48 


-1.95 


uu -^ bbh^ 


0.168 


0.164 


-2.56 


-1.95 


ss -^ bW 


0.028 


0.028 


-2.26 


-1.95 


cc -^ bbh^ 


0.013 


0.012 


-2.20 


-1.95 


gg ^ bW 


35.647 


33.734 


-5.37 


-1.95 


E [PP ^ bbh') 


35.963 


34.042 


-5.34 


-1.95 



Table 8.1: Hadronic cross sections for 66/i'^-production at the parameter point 
SPSla' (see appendix A. 2). 



table 8.1. It is also given separately for each partonic subprocess. As described 
before, the outgoing bottom-quark jets are required to have a high transverse 
momentum, so that large logarithms are avoided and the background processes, 
where the Higgs boson does not radiate off a bottom quark, are reduced. To this 
end a cut on the bottom quarks. 



PT{b,b) >20 GeV 



(8.6) 



was applied to obtain these results. The same cut will also be used for all other 
cross sections of this section. 

As one can see in the table, the dominant contribution originates from the 
gluon-fusion process and the quark-anti-quark-annihilation processes are sup- 
pressed by two orders of magnitude. Hence their contribution is negligible and 
in the following analysis only the gluon-fusion subprocess is considered. This 
large difference in the cross sections is due to the fact that the quark-anti-quark- 
annihilation process can only proceed via the s-channel diagram shown on the 
left-hand side of Fig. 8.1(a). It contains a propagator suppression from the in- 
termediate gluon which must carry at least the energy to procuce the final-state 
Higgs boson and the two bottom quarks. In contrast the gluon-fusion subprocess 
also contains a t-channel diagram (Fig. 8.1(b)) which does not suffer from such 
a suppression. In fact, if one takes only the s-channel diagram on the right-hand 
side of Fig. 8.1(a) into account, the cross section of the gluon-fusion contribution 
is of comparable size (cri = 1.103 fb) to the one of quark-anti-quark annihilation. 

In the following plots the effect of varying MSSM parameters on the SUSY- 
QCD contributions is investigated. To that end a parameter point with a fairly 



90 



Chapter 8. Higgs-Boson Production in Association with Heavy Quarks 






10 







-10 



-20 



-30 



-40 




400 600 800 1000 1200 1400 1600 1800 2000 



M 



SUSY 



[GeV] 



Figure 8.3: Partonic cross section differences for the process gg -^ bbh^, using 
t/3 = 30, as a function of Msusy = Mq = Mjy = Mfj. All other parameters take 
the values given in eq. (8.7). 



light SUSY spectrum was chosen, namely 

tua = 200 GeV 
/i = 300 GeV 



Msusy = Mq = M^ = M^ = 250 GeV 
m-g = 400 GeV . 



^.7) 



The MSSM parameters were then varied around this point. The renormalization 
scale, which appears in a^ (see eq. (4.9)), was set to fi^i = 2mi, + m/jO. As the 
contribution of the quark-anti-quark annihilation diagrams is negligible compared 
to the gluon-fusion subprocess, only the latter one is considered in the following. 
Also for simplicity the quoted cross section differences are partonic ones with a 
center-of-mass energy of vl = 500 GeV. 

For the first plots tp is set to the large value 30. In the plot given in Fig. 8.3 
a common mass scale Mgusy, where all soft-supersymmetry breaking masses in 
the squark sector take the same value, is chosen. According to chapter 4.3, where 
the bottom-quark Yukawa coupling was studied, one would expect the universal 
corrections, which are parametrized in Am;,, to give the dominant contribution. 
This is indeed the case for almost all MsusY-values. Also the decrease of the 
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Figure 8.4: Partonic cross section differences for 6&/i''-production via gluon fusion 
in tfie large t^-regime (t^ = 30) as a function of /i. For the value of the other 
parameters see eq. (8.7). 



corrections with growing SUSY mass scale, which is predicted by eq. (4.17) to 
fall off as TTT — , can be seen in the plot. Only for rather small values of Msusy a 



"-'SUSY 



deviation from this behavior occurs. Other terms contribute significantly in this 
region and lead to smaller cross-section differences than one would expect from 
the Arub terms alone. 

The numbers for the second plot (Fig. 8.4) are also calculated in the regime 
of large t/3, but now /i is varied. For small values of /i, the A?72;,-corrected tree- 
level result and the full one-loop cross section coincide and show the expected 
linear rise with /i. When /i becomes large, and thus the off-diagonal elements in 
the sbottom mixing matrix lead to a larger split between the lighter and heavier 
sbottom, this behavior changes and leads to a decelerated increase with /i. Also 
other terms begin to contribute to the cross section in a significant way and 
induce a deviation of the full one-loop result from the Anib corrections by up to 
20%. 

The effect of varying tp is studied in Fig. 8.5. The ATTi^-corrected tree-level 
result grows linearly with tp as predicted from eq. (4.17). In the small t/3-regime 
it approximates the full one-loop result rather well with a deviation of only about 
one percent. For larger values of t^ the complete one- loop corrections begin to 
deviate and additional contributions lead to a slower rise. Finally, the absolute 
value of the full corrections slightly decreases again. This is the same effect which 
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Figure 8.5: t/3- dependence of the partonic cross section differences for gg 
The values of all other parameters are given in eq. (8.7). 



bbh^. 



was already observed on the left-hand side of Fig. 8.3. As the common value of 
the soft supersymmetry-breaking masses MgusY was chosen to be 250 GeV we are 
exactly in this regime. A higher value for MgusY leads to a one- loop cross- section 
difference which coincides with the Amt-corrected tree-level one over the whole 
range of t/3. To verify this a value of Msusy = 400 GeV was chosen to obtain 
Fig. 8.6. Additionally, the gluino mass was set to nig = 640 GeV such that the 
ratio of the two masses is the same as in the parameter set eq. (8.7). In this 
case the discrepancy AAm^ between Ai and Aq stays below one percent for all 
t/3- values. 

In the final two plots the behavior of the cross-section differences in the small 
t^- regime, namely for t^ = 6, is studied. Firstly, in Fig. 8.7 the common mass 
scale MgusY of the soft supersymmetry-breaking masses appearing in the squark 
sector is varied. The A?Tib-corrected tree-level result is a good approximation of 
the full one-loop result over the whole mass range. The difference is about one 
percent for small Msusy and rapidly vanishes for larger values. 

The last plot (Fig. 8.8) depicts the dependence of the cross-section differences 
on yU for t^ = 6. It shows a similar behavior as the same plot for large t^ 
(Fig. 8.4). For small values of /i the A?72b-corrected tree-level result and the full 
one-loop calculation coincide, while for larger ones significant deviations occur. 
In the small t/3-regime the absolute value of the Aq corrections even grows slightly 
larger than the linear behavior of eq. (4.17), which is due to additional effects 
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Figure 8.6: t/3-dependence of the partonic cross section differences for gg 
For this plot a shghtly higher Msusy = 400 GeV and m^ 
than the one of eq. (8.7). 
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Figure 8.7: Partonic cross section differences for 66/i''-production via gluon fusion 
for t/j = 6 as function of a common mass MgusY for the soft supersymmetry- 
breaking terms. All other parameters take the values given in eq. (8.7). 
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200 400 600 800 1000 1200 1400 1600 1800 2000 

1^ [GeV] 

Figure 8.8: /i-dependence of the partonic cross section differences for bbh^- 
production via gluon fusion in the small t^-regime {tf^ = 6). The values of all 
other parameters is given in eq. (8.7). 

originating from the growing mass splitting between the two sbottoms. The one- 
loop corrections in contrast decrease again, once the absolute value has reached 
a maximum at about 1500 GeV, and even change sign and become positive. So 
in this parameter region true one-loop corrections contribute significantly. 



8.5 Numerical Results for tth^ 

The numerical results for the second process of Higgs production in association 
with heavy quarks, pp -^ tih^, are presented in this section. In table 8.2 the 
hadronic cross section for the MSSM reference point SPSla' is denoted. It is 
given separately for each partonic subprocess which contributes to the tt/i°-final 
state. 

In this case the quark-ant i- quark annihilation diagrams give a contribution 
which is of comparable size to the gluon-fusion ones. On the partonic level the 
same analysis as in the previous section for /i°-production in association with 
a bottom quark-anti-quark pair holds. The quark-anti-quark-annihilation pro- 
cesses are suppressed because there only a propagator-suppressed s-channel di- 
agram exists, while the gluon-fusion subprocess also proceeds via a t-channel 
diagram which does not suffer from such a suppression. After the convolution 
with the parton distribution functions the situation however changes. The gluon 
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Partonic subprocess 


o-Q [fb] 


ai [fb] 


Ai [%] 


dd -^ tihP 


42.7 


37.6 


-11.77 


uu -^ tihP 


71.9 


63.4 


-11.81 


ss -^ tih^ 


7.5 


6.6 


-11.58 


cc -^ tih^ 


2.8 


2.5 


-11.53 


99 -^ tW 


273.7 


264.7 


-3.30 


E {PP - tth^) 


399.0 


374.8 


-5.96 



Table 8.2: Hadronic cross sections for tt/i'^-production at the parameter point 
SPSla', which is defined in appendix A. 2. 



densities in the proton show a much steeper fall with growing parton-momentum 
fraction x than the sea-quark ones. As the top quarks are much heavier than 
the bottom quarks, also the energy and thus the x of the incoming partons must 
be larger to be above the threshold for tt/i°-pro duct ion. For this final state 
it approximately compensates the effect from the propagator suppression. The 
gluon-fusion process is still the dominant production mode, but all processes need 
to be taken into account for a complete analysis. 

In the following plots the effect of varying MSSM parameters on the SUSY- 
QCD contributions is investigated. To that end the same parameter point as in 
the previous section with a fairly light SUSY spectrum was chosen, namely 



ruA = 200 GeV 
t/j = 6 
fx = 300 GeV 

Mq = Mfj = Mi) = 250 GeV 
m-g = 400 GeV . (8.8) 

The MSSM parameters were then varied around this point, tp = 6 was kept 
fixed for all plots of this section. The renormalization scale, which enters as via 
eq. (4.9), and the factorization scale were set to fiR = hf = '2mt + m/jo. The 
hadronic cross-section calculations were performed for the LHC with a proton- 
proton center-of-mass energy of 14 TeV. 

First a common mass scale Msusy, where the soft supersymmetry-breaking 
squark mass terms all take the same value, is introduced and varied between 
200 and 2000 GeV, as shown in Fig. 8.9. The difference between the tree-level 
cross section and the Amf-corrected one falls off as ytt — , as expected from the 

form of the Arrit term given in eq. (4.30). The total one-loop contributions show 
a similar decrease, but with a larger coefficient which leads to a much steeper 
descent. This originates from the fact that the Arrit term only includes vertex 
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Figure 8.9: Hadronic cross section differences for the process pp — > tih'^ as a 
function of MgusY = Mq = Mjy = Mfj. For the value of the other parameters 
see eq. (8.8). 



corrections to the tih^ vertex. Yet there are many other one-loop diagrams which 
also contribute and lead to the modified behavior. In contrast to Arrih, which 
is enhanced by a factor t^, the Anit corrections are suppressed by j- and their 
numerical effect is expected to be smaller. For small SUSY masses threshold 
effects of the squark masses induce a deviation from the scaling with 



Mi 



In the second plot (Fig. 8.10) the dependence of the relative corrections on /i, 
the mass parameter mixing the two Higgs doublets, is presented. Also in this case 
the Am^-corrected tree-level cross section is not a good approximation. Whereas 
the term rises with growing /i, the full one-loop correction decreases in this case. 
The slope is constant over a large range of fi. Only for bigger values, when the 
off-diagonal entries in the squark mixing matrices become very large and yield 
an additional contribution, also the gradient increases. 

Finally, in Fig. 8.11 tp is varied. The corrections to the Amt-corrected 
tree- level result fall off with growing tjj. This is again the expected behavior 
of eq. (4.30). The full one-loop corrections are significantly larger in size. They 
show a mild dependence on this parameter, with a maximum at around t/? = 15. 
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Figure 8.10: Hadronic cross section differences for tt/i^-pro duct ion as a function 
of /i. The values of all other parameters are given in eq. (8.8). 
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Figure 8.11: t/j-dependence of the hadronic cross section differences for pp 
All other parameters take the values given in eq. (8.8). 



tth^. 



Chapter 9 

Quartic Higgs Coupling at 
Hadron Colliders 



In this chapter another possibihty is investigated to test the means of electroweak 
symmetry breaking. This is achieved by measuring the quartic Higgs couphng 
and hence fully determining the Higgs potential [22]. 

After the discovery of a light Higgs boson the next step will be to study its 
properties, including its couplings to other particles. At the planned Interna- 
tional Linear Collider (ILC) measuring these couplings with high precision will 
be possible for all Standard Model bosons and fermions [138]. Furthermore, if 
super symmetric particles are found, the coupling of the Higgs to charginos and 
neutralinos can be measured precisely [139]. To fully understand electroweak 
symmetry breaking it is important to measure the Higgs self-couplings and to 
thereby determine the parameters of the Higgs potential. 

9.1 Higgs potential 

The Higgs potential of the Standard Model was already given in eq. (2.2). In 
this model the trilinear (A3) and quartic (A4) Higgs self-coupling are related to 
the Higgs mass via 



-Simlr , —Simlr Aq 

I2' 



I'M X -ovrnjj A3 



A3 = ^ A4 = ^ = ^, (9.1) 



where v is the vacuum expectation value of the Higgs field. 

In models with more than one Higgs field, like the MSSM with its two fun- 
damental Higgs doublets, the relations between the trilinear and quartic Higgs 
coupling can change significantly. For the lighter CP-even MSSM Higgs boson 
h^ the ratio of the self-couplings is 

^ = ^'J±^^ (9.2) 

A4h0 C2a 
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where v = \/v\ + v\ and V\ and f 2 are the vacuum expectation values of the two 
Higgs fields. If the parameter vtij^ is sufficiently large there is a mass splitting 
between /i*^ and the remaining Higgs sector. Additionally, the angles a and /? 
are related via s^ — —Cjs in this limit. Therefore ^^^^ approaches 1 and the 
/i^-coupling becomes Standard Model-like. 

In this chapter we will not refer to the MSSM as our underlying theory. 
Instead we use an effective theory whose particle content is the same as the one 
of the Standard Model. Its Higgs sector also contains one doublet but the trilinear 
and quartic couplings are left as free parameters of the theory. In this way we 
are not restricted on a specific model but can accommodate for many different 
ones. Such deviations from the Standard Model couplings can for example be 
generated when higher-dimensional powers of the Higgs doublet are added to the 
potential as shown in chapter 2.2.2. Taking the first two higher-order terms into 
account the Higgs self-couplings become 



A, =A, 1 






A.^A.(l + S^ + i^ . (9.3) 

The additional terms are suppressed by A which is the scale where new physics 
sets in. Both self-couplings receive different contributions from the additional 
terms. In general, the self-couplings may even become negative. The stability 
of the Higgs potential is guaranteed if the highest non-vanishing term in the 
potential has a positive sign. All other terms can have arbitrary values as long 
as the ground state has a non-vanishing vacuum expectation value to break the 
electroweak symmetry. 

9.2 Trilinear Higgs coupling 

As we will see below it is essential for the measurement of the quartic Higgs 
coupling to know the value of the trilinear Higgs coupling as precisely as possible. 
For a Higgs boson with a mass larger than 150 GeV this coupling can be extracted 
at the LHC [140, 141, 142]. At least two Higgs bosons must be produced to 
measure the three-Higgs coupling. At hadron colliders this is performed via a 
gluon fusion process. In this process two distinct types of diagrams appear, as 
shown in Fig. 9.1. Either (a) an intermediate Higgs boson is produced via a three- 
point top-quark loop diagram which couples to the two final-state Higgs bosons 
and contains the required trilinear coupling, or (b) the particles couple via a four- 
point box-type top-quark loop. For the detection of the Higgs bosons the decay 
channel into two W^W~ boson pairs is analyzed [142]. Two or three of the four 
W bosons are required to decay leptonically into a lepton and a neutrino to have 
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Figure 9.1: Leading-order types of Feynman diagrams contributing to the process 
gg^HH 



a clear detector signal and the other two or one, respectively, decay hadronically 
into two jets. 

Additional information can be obtained from the kinematic distributions of 
the Higgs bosons. Not only the total cross section carries information on the 
trilinear Higgs coupling but also the differential hadronic distribution with respect 
to the invariant mass of the final state. This can only be calculated correctly if 
the top-quark loop is fully taken into account. Using the infinite top-quark mass 
limit and an effective gluon-gluon-Higgs coupling will yield completely incorrect 
results, as was shown in ref. [142]. 

Using this information it was evaluated [142] with which precision the trilinear 
Higgs coupling for Higgs bosons heavier than 150 GeV can be measured. In the 
beginning of LHC a non-zero value can be established with a confidence level 
of 95% after having accumulated a luminosity of 300 fb~^. After a luminosity 
upgrade the self-coupling can be measured with a precision of up to 20% at 
the 95% confidence level using an integrated luminosity of 3 ab~ . At a future 
high-energy Very Large Hadron Collider (VLHC) with a hadronic center-of-mass 
energy of 200 GeV, the measurement of A3 can be performed with an uncertainty 
of about 10% at a confidence level of 95% after having accumulated a luminosity 
of 1 ab-^ 

In contrast a linear collider can measure the three-Higgs coupling for small 
Higgs masses of around 120 GeV, but possibly not for higher ones [143, 144]. 
The mass region below 140 GeV is more difficult for hadron colliders, because 
the dominant decay mode for the Higgs boson is a bottom quark-anti-quark 
pair which has large QCD backgrounds. Anyway, it is still accessible using rare 
decays [145]. The creation of the two Higgs bosons at an e^e~ linear collider 
happens as double-Higgs strahlung in association with a Z boson or as WW 
double- Higgs fusion in association with a uu-paii as shown in Fig. 9.2 and Fig. 9.3, 
respectively. For a Higgs mass of 120 GeV the planned International Linear Col- 
lider (ILC) with a center-of-mass energy of 500 GeV will be able to measure the 
trilinear Higgs coupling with an accuracy of 20% within one standard deviation 
after having accumulated a luminosity of 1 ab~^. A proposed Compact Linear 
Collider (CLIC) with a center-of-mass energy of 3 TeV finally could measure a 
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Figure 9.2: Leading-order types of Feynman diagrams contributing to the process 
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Figure 9.3: Leading-order types of Feynman diagrams contributing to the process 
e^e~ -^ vvHH 



180 GeV Higgs boson with a precision of 8% within one standard deviation using 
5 ab^ of integrated luminosity. 

In general the proposed future generation of hadron (VLHC) and linear (CLIC) 
colliders will be able to measure the trilinear Higgs coupling with an accuracy 
of O (10%). A combination of both collider types thereby covers the whole mass 
range where a Standard Model Higgs boson is expected to be found, as derived 
from electroweak precision analyses [20]. 



9.3 Quartic Higgs coupling 

As the production of two Higgs bosons is needed for measuring the trilinear Higgs 
coupling, three final-state Higgs bosons are necessary for a measurement of the 
quartic Higgs coupling. 

Three-Higgs production at linear colliders has already been studied in ref. [146]. 
It was found that even at CLIC the cross section is too low, with only about five 
three-Higgs events per year being produced there at a center-of-mass energy of 
10 TeV. Hence it will be impossible to determine the quartic Higgs coupling at 
the next generations of linear colliders. 

At hadron colliders the three Higgs bosons are dominantly produced via gluon 
fusion and an intermediate top-quark loop, like in the two-Higgs case. Four 
distinct topologies appear as shown in Fig. 9.4: (a) continuum production of 
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Figure 9.4: Leading-order types of Feynman diagrams contributing to the process 
gg ^ HHH 



three Higgs bosons via a five-point top-quark loop, (b) production of two Higgs 
bosons via a box-type loop and subsequent decay of one of the Higgs bosons via 
the trilinear self-coupling into two Higgs bosons, and finally the production of one 
intermediate Higgs boson via a three-point loop. This can either (c) decay via 
a chain of two three-Higgs couplings or (d) through one quartic Higgs coupling. 
Only the last diagram type contains the coupling we want to measure. 

Looking at the diagrams it is clear that a precise knowledge of the trilinear 
self-coupling is necessary to obtain results on the quartic self-coupling. The 
process is also very sensitive to the top-quark Yukawa coupling which must be 
known very well. In the numerical analysis we have also included the diagrams 
where the top-quark loops are replaced by bottom-quark loops. The contribution 
of these diagrams is however less than one percent. 

The total cross section as function of the Higgs mass is shown in Fig. 9.5 for 
the (a) LHC and a (b) 200 TeV VLHC In the following the Higgs boson mass is 
set to 120 GeV for all cross sections and plots. 

In Fig. 9.6 the dependence of the total cross section on the trilinear and quartic 
Higgs coupling is shown. The values of the trilinear and quartic self-couplings 
are varied between minus and plus two times the Standard Model value. One can 
clearly see the strong dependence on A3. The variation on A4 is much smaller, 
as one can see in more detail in Fig. 9.7. For positive values of A3 and A4 the 
variation of the cross section stays below 20% at both the LHC and the VLHC. 
Including negative values of A4 induces changes in the cross section of up to a 
factor 2. For negative values of A3 the absolute variation as a function of A4 stays 
at the same order of magnitude. While the total cross section depends strongly 
on A3, the relative variation with A4 is heavily suppressed. 

As the quartic Higgs coupling contributes only to the single diagram Fig. 9.4(d) 
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Figure 9.5: Total hadronic cross section for the triple-Higgs production process 
via gluon fusion with Standard Model couplings as a function of the Higgs boson 
mass. 



this behavior is expected. It is also reflected in the partial contributions from the 
different diagram types. Taking into account all diagrams leads to a Standard 
Model cross section of 6.25 ■ 10~^ fb at the LHC. Using Standard Model couplings 
the five-(Fig. 9.4(a)), four-(Fig. 9.4(b)) and three-point (Fig. 9.4(c,d)) loop dia- 
grams alone yield a cross section of (17.07, 8.20, 0.46) ■ 10^^ fb, respectively. The 
small size of the triangle-type diagram results from a suppression factor by the 
intermediate Higgs propagator. If one only takes this diagram type into account 
and sets either A4 or A3 to zero, a cross section of 0.17 ■ 10~^ fb from the trilinear 
self-coupling and of 0.08 ■ 10~^ fb from the quartic self-coupling is obtained. For 
the VLHC the partial cross sections have similar ratios. So the diagram which 
contains the quartic self-coupling is almost two orders of magnitude smaller than 
the total cross section. 

As one would expect from the evaluation of the trilinear Higgs coupling in 
two-Higgs production [142] the interference between the pentagon and the box 
diagrams is indeed destructive for positive values of A3. This results in the large 
increase of the cross section when A3 gets smaller and therefore the contribution 
from the box diagram diminishes as one can see in Fig. 9.6. For negative values of 
A3 the single three-Higgs vertex in the box diagrams changes sign and therefore 
makes this interference constructive, leading to the sharp rise of the cross section. 
The trilinear and pentagon diagrams interfere constructively, but because the 
triangle and box contributions have a more similar kinematic configuration, the 
destructive interference between those two results in a slight decrease of the cross 
section with growing A4. Only for A3 = 0, where the box diagrams do not 
contribute any longer, the behavior of the cross section reverses and rises with 
increasing A4. The relative signs of the different topologies can be understood 
analytically by using the low-energy theorem for the leading form factors [147]. 
An expansion is performed in the ratio ^^^ with a partonic center of mass energy 
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Figure 9.6: Total hadronic cross section for triple-Higgs production via gluon 
fusion as a function of the trilinear and quartic Higgs coupling normalized to the 
Standard Model values. A color of green denotes the Standard Model value. A 
deviation of plus and minus 20 % is signified by red and blue color, respectively. 
The maximum values obtained in the scanning interval are colored white and 
black, respectively, using a linear color gradient for intermediate values. The 
Standard Model point is additionally marked by a blue cross. A Higgs boson 
mass of 120 GeV was used to obtain this plot. 
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Figure 9.7: Variation of the hadronic cross section for iJiJiJ-production as a 
function of A4 normalized to the value where A4 takes its Standard Model value. 
The Higgs mass was fixed to 120 GeV. 
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s ~ ml. The top-quark mass and the top-quark Yukawa couphng are both 
denoted by nit. These form factors are basically the squared matrix element 
without any couplings or additional propagators which are not part of the loop 
integral. They are obtained starting from the top loop in the gluon self-energy. 
An additional Higgs boson can be attached to the loop by using the following 
recursion relation: 

F{n+i)H = mt . (9.4) 



onit nit 




triangle 'J 


\ nif 



This relation yields 

-'^pentagon -'^box -'^triangle ^7 I 2 J ' l'^"^/ 



Therefore the structure of constructive and destructive interferences in the dia- 
gram types is explained by the relative minus sign in front of the box-type term. 

In the case of two-Higgs production the information on the value of A3 was not 
only encoded in the total cross section. Also the differential hadronic cross section 
with respect to the invariant mass carries information on A3. The same is true 
for three-Higgs production. In Fig. 9.8 the normalized cross section as a function 
of the partonic center-of-mass energy is shown. The trilinear and quartic self- 
coupling both take the values 0, 1 and 2 times their respective Standard Model 
value. When varying A3 the position of the peak changes significantly and an 
extraction of this coupling is possible, as was already found in the analysis of 
two-Higgs production [142]. When changing A4, and keeping A3 constant, the 
size of the shift is about an order of magnitude smaller. Additionally for A3 = 
the order of the A4/ASM = 0, 1, 2 peaks is inverted which is due to the different 
sign in the interference as explained above. 

The total hadronic cross section for triple-Higgs production via gluon fusion 
using Standard Model Higgs self-couplings is 6.25 ■ 10~^ fb at the LHC for a 
Higgs boson with a mass of 120 GeV. This rate is too low to be measurable even 
for the high-luminosity mode of the LHC. At the VLHC with a center-of-mass 
energy of 200 GeV the cross section is 9.45 fb so three-Higgs production might 
be observable at this future collider. The rather strong dependence of the total 
cross section on A3, especially for values smaller than the Standard Model value, 
allows to extract the value of this coupling, thereby possibly confirming the two- 
Higgs production result. The variation in A4 is much smaller, typically below 
20%. Hence the extraction of this coupling is much harder. If one takes into 
account the theoretical uncertainties from missing higher-order corrections and 
the experimental error on the measurements of A3 and the top-quark mass, the 
chances to be able to extract the quartic Higgs coupling are tiny. 

Also in the differential cross section there is a clear effect on the peak position 
when varying A3. This shift will be the mode to extract the trilinear coupling 
in double-Higgs production. The size of the shift for a variation of A4 is much 
smaller. If the errors on the measurements of A3 and the top-quark mass are 
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Figure 9.8: Differential hadronic cross section for three-Higgs production witli 
respect to tlie partonic center-of-mass energy, normalized to the respective total 
cross section. The mass of the Higgs boson was set to 120 GeV. 
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again taken into account, the extraction of the quartic Higgs self-couphng looks 
challenging. 



Chapter 10 
Conclusions 



In this thesis production processes for Higgs bosons at hadron coUiders were 
considered. In order to facilitate the computation of hadronic cross sections from 
a large number of complicated parton processes, a computer code was developed. 
The calculation of hadronic cross sections, in particular for the production of 
supersymmetric Higgs bosons at the LHC in various processes was examined in 
the first part of this thesis. One-loop SUSY-QCD corrections, i.e. corrections 
with squarks and gluinos running in the loop, were calculated and the numerical 
results discussed. In the second part triple-Higgs production in an effective theory 
was examined and the question whether the quartic Higgs self-coupling can be 
measured at hadron colliders, pursued. 

The calculation of cross sections for processes which contain hundreds of sin- 
gle Feynman diagrams is not possible without the help of automated tools. For 
partonic cross sections there are already programs, like the packages FeynArts, 
FormCalc and LoopTools. The latter one was extended, so that the five-point 
loop integrals are now included which makes a calculation of 2 — i> 3 processes 
with these tools possible. Additionally, the numerical stability of the loop inte- 
grals was improved. To obtain hadronic cross sections, which are the observables 
which will be measured at the LHC, the partonic cross sections must be con- 
voluted with the parton distribution functions. Therefore a computer program, 
called HadCalc, was written which performs this task. Using this program it is 
possible to calculate both integrated and differential partonic and hadronic cross 
sections of processes which were generated by FormCalc beforehand. The cross 
sections can be differential with respect to the invariant mass of the final state, or 
rapidity or transverse momentum of one of the outgoing particles. Additionally, 
the possibility to apply cuts on the rapidity, the transverse momentum, and the 
jet separation of the final-state particles was implemented. 

This program was then used to calculate SUSY-QCD corrections to various 
Higgs-boson-production processes at the LHC. The first considered process was 
the production of a charged Higgs boson in association with a W boson via bot- 
tom quark-anti-quark annihilation. It is known that the bottom-quark Yukawa 
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coupling receives large one- loop corrections for large tp. They are universal and 
can be parametrized via the variable Aruf, and summed up to all orders in per- 
turbation theory. The numerical analysis showed that in the large-tp regime this 
term indeed represents the dominant contribution and the one-loop cross section 
is well approximated by the A?n6-corrected tree- level result. For small tf3, also 
other terms play an important role. The leading subdiagram of this region was 
identified and its analytical behavior studied. These contributions can reach a 
significant size, yielding corrections of up to 50% for certain parameter combina- 
tions. It is a true one-loop result and cannot be taken into account by an effective 
tree-level coupling. 

Furthermore, the production of an h^ via vector-boson fusion was calculated. 
This process has a clear experimental signature of two jets in the forward re- 
gion of the detector. In the theoretical analysis this phase-space region was 
selected by applying corresponding cuts. The numerical size of the SUSY-QCD 
corrections to this process is of O{10~^), and therefore significantly below the 
experimental uncertainty which LHC will be able to reach. This smallness of the 
corrections could be explained by cancellations between different one-loop Feyn- 
man diagrams and the appearance of suppression factors. Hence, they do not 
induce a sufficiently large modification of the cross section, which would allow 
to distinguish between the SM and the MSSM in the Higgs sector in this way. 
Additionally, /i°-pro duct ion with two final-state jets, where one or two of the in- 
coming partons are gluons, was considered. This constitutes a background to the 
above-mentioned vector-boson-fusion process. Thus for the calculation the same 
cuts were applied. The total contribution of these processes is smaller than the 
total vector-boson-fusion cross section by more than four orders of magnitude, 
so they can be safely neglected in an analysis of /i°-production via vector-boson 
fusion. 

Moreover, the associated production of an h^ with a heavy, i.e. bottom or top, 
quark-ant i- quark pair was studied. These processes also form discovery channels 
for the Higgs boson. Additionally, they can be used to measure the respective 
Yukawa couplings. In the case of bottom quarks the final-state jets are required 
to have large transverse momenta to avoid the appearance of large logarithms 
and corresponding cuts on the phase space were applied. For the top quarks 
no such cuts are necessary and therefore the calculation of the cross section was 
conducted over the full phase space. Both times the SUSY-QCD corrections 
provide a significant contribution to the total cross section, which must be taken 
into account. For final-state bottom-quarks the A?Ti(,-corrected tree-level cross 
sections give a good estimate of the full one-loop result in large regions of the 
MSSM parameter space. Only when the off-diagonal elements of the sbottom 
mixing matrix become large, does this approximation break down. Here other 
terms also give significant contributions, which lead to a sizable change in the 
numerical result. As the sign of the two terms differs, the overall size of the one- 
loop corrections is reduced. Nevertheless, they provide a significant contribution 
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to the total cross section, which can reach up to 40% for certain combinations of 
MSSM parameters. In the top-quark case the full one-loop calculation is never 
well approximated by the Ami-corrected tree-level cross section. As the Arrit 
term contains a suppression of j- instead of the t^-enhancement of Arrif,, this 

behavior is expected. The total size of the corrections to tt/i'^-production is of 
the order of several percent. 

In the last part of this dissertation the production of three Higgs bosons at 
hadron colliders was considered. This process can be used to extract the quartic 
Higgs self-coupling, thereby determining the Higgs potential. In this calculation 
not the MSSM was used as the underlying model, but an effective theory deduced 
from the Standard Model, where the three- and four-Higgs self-couplings were 
left as free parameters. The numerical analysis showed that the cross section 
is too small to be measured at the LHC A future Very Large Hadron Collider 
with a projected center-of-mass energy of 200 TeV would produce enough three- 
Higgs events. However, because of the interference structures of the different 
diagrams which contribute to this final state, the extraction of the quartic Higgs 
coupling from the invariant-mass distribution will be seriously challenging. This 
is especially true if the theoretical uncertainties and the statistical errors on the 
measurement of the trilinear self-coupling and the top-quark mass are taken into 
account. 



Appendix A 

Choice of Parameters 

A.l Standard Model Parameters 

The parameters of the Standard Model have been measured by various experi- 
ments. Their current best average is [17]: 

• Masses: 
Quarks: 



mu = 53.8 MeV 


TTic = 1.5 GeV 


mt 


= 178 GeV 


ma = 53.8 MeV 


rris = 150 MeV 


rrib 


= 4.7 GeV 


Leptons: 








me = 510.999 keV 


m^ = 105.658 MeV 


m^ 


= 1.777 GeV 


m,y^ = 


rriu^ =0 


rriy^ 


= 


Gauge bosons: 








mw = 80.45 GeV 


mz = 91.1875 GeV 


m^ 


= nig = 


Couphng constants 








a-i (g2 = 0) = 137.0359895 «f ^(Q^ = 


= m|) = 


= 0.1172 . 



The masses of the first-generation quarks rriu and rrid are effective parameters. 
They were chosen in a way that the vacuum polarization of the photon, which 
was determined from experimental data via a dispersion relation and is known 
more exactly, is correctly reproduced [148]. 

The calculations in this thesis were performed with the parameters quoted 
above. Only the masses of the light quarks m„, m^, rric and m^ were set to zero 
exactly as their contribution is negligible. Also the lepton masses are mentioned 
here only for completeness, as they do not enter any diagram. 

115 



116 Appendix A. Choice of Parameters 

A.2 SPA scenario of the MSSM 

Even after adding experimental bounds and eliminating regions of the parameter 
space which are disfavored by theoretical arguments the possible choices of the 
miscellaneous MSSM parameters are still plenty. In order to unify the conventions 
used in the calculations and to allow for a comparison of the results of different 
working groups, an effort was made to apply the same conventions and use certain 
points in the MSSM parameter space as benchmark scenarios. For this the SPA 
(Supersymmetry Parameter Analysis) conventions [16] were established. 

One of the favored reference points of the MSSM parameter space in the 
SPA conventions is called SPSla'. It is defined in the minimal supergravity 
(mSUGRA) scenario. This scenario assumes that all parameters are real, a uni- 
fication of the gauge couplings happens at the GUT scale and the soft super- 
symmetry-breaking terms are universal at the high-energy scale. Therefore the 
number of parameters is greatly reduced to only four. They are a common scalar 
mass Mq, a common gaugino mass M1/2, a common trilinear coupling Aq and the 
ratio of the Higgs vacuum expectation values tp. Additionally the sign of fi is 
not fixed and can be chosen freely. For SPSla' these variables take the following 
values 

Mo = 70 GeV M1/2 = 250 GeV Ao = -300 GeV 

tp{M) = 10 sign(/i) = +1 . (A.l) 

This parameter point was chosen in such a way that it is compatible with all 
current experimental bounds. The mass parameters are defined at the GUT 
scale and then evolved via renormalization group equations (RGEs) to the SPA 
scale M = 1 TeV which is also the scale where t^ is specified. 

Using this procedure the MSSM parameters take the following values at the 
SPA scale [16, 80] 

t(3 =10 fi =402.87 GeV m^ =431.02 GeV 

Ml =103.22 GeV M2 =193.31 GeV M3 =572.33 GeV 

A^'2 = - 784.7 GeV A^ = - 1025.7 GeV Al'^ = - 449.0 GeV 

Al = - 535.4 GeV A^ = - 938.5 GeV A^ = - 445.5 GeV 

M^'^ =526.9 GeV Mi'^ =507.7 GeV mJ^'^ =505.5 GeV 

Mj =471.3 GeV M? =384.6 GeV M| =501.3 GeV 

mI'^ =181.3 GeV Mg'^ =115.6 GeV 

L/ Jrt 

Ml =179.5 GeV M| =109.8 GeV . (A.2) 

The RGE evolution leads to different values for the soft-breaking mass parameters 
and trilinear couplings of the first two and the third generation. In the above- 
mentioned table the generation index is denoted by the superscript. 
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These parameters must be interpreted as parameters in the DR renormal- 
ization scheme. As most of the calculations in this thesis are done using OS 
renormalization a further conversion step is necessary. Using these DR parame- 
ters the masses of all supersymmetric particles are computed at the one-loop level. 
Then a minimal set of these masses is chosen and the OS parameters are calcu- 
lated using the tree-level relations between parameters and masses. With this 
procedure the physical masses of the particles in both renormalization schemes 
are equal. Thus a meaningful comparison of the results of calculations in both 
schemes is possible. 

This conversion procedure yields the following MSSM parameters in the OS 
scheme [80]: 

tp =10 /i =399.26 GeV m^ =431.02 GeV 

Ml =100.11 GeV Ma =197.55 GeV M3 =612.85 GeV 

A^'2 = - 784.7 GeV A]'^ = - 1025.7 GeV A]'"^ = - 449.0 GeV 

Al = - 535.4 GeV A^ = - 938.5 GeV A^ = - 445.5 GeV 

M^ =565.97 GeV M^ =546.78 GeV M^ =544.95 GeV 

M? =565.91 GeV M? =546.84 GeV M? =544.97 GeV 

M| =453.05 GeV M? =460.52 GeV M? =538.13 GeV 

M} =184.12 GeV Ml =118.02 GeV 

Ml =184.11 GeV M| =117.99 GeV 

Ml =182.18 GeV M| =111.29 GeV . (A.3) 

These are the parameters used for the calculations in this thesis when the SPSla' 
parameter set is referred to. 



Appendix B 

Basic Principles of 
Supersymmetry 

B.l Poincare group 

Every point in four-dimensional space-time of Minkowski space is characterized 
by a contravariant vector which is defined by 

x^' = (x°, x\x^, x^) = {t, x) (B.l) 

as generahzed space coordinate. With the metric tensor 

^M- = ^'^" = diag (1,-1, -1,-1) (B.2) 

a covariant vector 

Xf, = g/^yx" = (t, -x) (B.3) 

can also be defined. 

Here and in the following Greek indices run from to 3 and Latin ones from 
1 to 3 except where denoted otherwise. Einstein's sum convention is implicitly 
assumed, i.e. indices which appear once as covariant and once as contravariant 
index are summed over. Additionally natural units are used where h = c = 1. 

Derivatives with respect to generalized space coordinates can be abbreviated 
as 

The momentum four-vector is defined as 

p^' = i^^' = {E,p) . (B.5) 
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B.2 Spinors 

B.2.1 Weyl spinors 

Two-component anticommuting objects ^a which transform under a matrix M of 

SL{2X) as 



c-^{M-');e r-(M-^)*%f (b.6) 



/9 
13'-' 



are called Weyl spinors. The spinor indices a, a, (3 and (3 can take the values 1 
and 2. Additionally the relations 

L = C ia = ea^e (B.7) 

hold. On the right-hand side the two-dimensional totally antisymmetric tensor 

1 for even permutations of {1,2} 

-1 0, 



e"^ = -ea/3 = ( °T i ) = ^ -1 for odd permutations of {1, 2} (B.8) 



else 
has been introduced. Along the same lines also three-dimensional 

{1 for even permutations of {1, 2, 3} 
-1 for odd permutations of {1, 2, 3} (B.9) 

else 

and four- dimensional versions 

{1 for even permutations of {0, 1, 2, 3} 
— 1 for odd permutations of {0, 1, 2, 3} (B.IO) 

else 

can be defined, which are needed later. 

If one generalizes the definition of the Pauli matrices a 

^'K?J) ^'K°"oO ^-G-0 '^-"^ 

to four dimensions via 

a>^ = {l,a) a'^ = (l,-a) , (B.12) 

then the Dirac equation can be written in two-component notation as follows: 

{^,pT^^P = mr {<y,v')^pf = ^ia ■ (B.13) 
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B.2.2 Dirac and Majorana spinors 

Out of two two- component Weyl spinors a four-component spinor 

^ = (P) (B.14) 

can be constructed. lit] = C,, ip is called a Majorana spinor, else ip is called Dirac 
spinor. The four-component equivalent to the Pauli matrices are the 7 matrices 
which are defined via the Clifford algebra 

{7/^, Y} = 7^7^^ + 7^7^ = 20^'" . (B. 15) 

Their hermitian conjugates are 

^Mt = YyY ^ (B.16) 

so 7*^ is hermitian and the 7* are antihermitian. There are different representa- 
tions of the 7 matrices which all fulfill eq. (B.15). The one corresponding to the 
form of the spinors in eq. (B.14) is called chiral representation and in that one 
the 7 matrices have the following form 

Additionally one defines 75 as 

,5 _ „• ,0„,1,,2„,3 chiral /—I 



representation \ U -11 / 

for which the following relations hold: 

{7^7'} = (75)' = . (B.19) 

Then the projection operators 

Pl^\{1~ 75) Pr = ^ (1 + 75) (B.20) 

yield the left- and right-handed part of a Dirac spinor, respectively 

i^L = Pii^ = (^q") ^r = Pr^ = (-°J . (B.21) 

The 4x4 tensor matrices a^^ are constructed from the 7 matrices via 

a'^^ = ^[7^71 • (B.22) 
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The spinor ijj which is adjoint to ip is defined as 



(B.23) 



Starting from eq. (B.13) the Dirac equation can now also be written in a four- 
component notation 



ilfj.P^ — rn)il) ^ if — m) ip = ■ 



[B.24:) 



In the following a few contraction identities and traces over 7 matrices are 
collected which are needed for the calculation of cross sections: 



# 


= a ■ b — ia^ya^h" 


Yi,i 


= 4 


Tr(l) 


= 4 


Tr(7^) 


= 


Tr (7^) 


= 


Tr (7'^7^) 


= 


Tr(7'^7^) 


= Ag^^ 


Tr (7^7'^7') 


= 


Tr (7^7'^7''7") 


= 4((7^V'^-(7^V" + ^^V 


Tr (7''7^7''7"7^) 


= -4,ie^'^p^ 


Tr(/7^..7V) 


= . 


odd number of 7's 





(B.25) 



B.3 Grassmann variables 

For the definition of the supersymmetric anticommutation relations in chapter 3.2 
anticommuting numbers, so-called Grassmann numbers, were introduced. The 
basic relation between two such numbers 9 and ^ is 



which is equivalent to 



oi = -ie 



{^,e} = o 



(B.26) 



(B.27) 



The last relation shows the fermionic character of these numbers. In particular, 
the square of any Grassmann number is zero. 

Grassmann numbers form an Abelian group under the operation of addition. 
The multiplication with ordinary complex numbers has the same properties as 
scalar multiplication of a vector, in particular the distributive law holds. 
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The integral over Grassmann numbers is defined as 

I d9{A + B9) = B , (B.28) 

where A and B are complex numbers. This leads to the expression for differen- 
tiation with respect to Grassmann numbers 

^{A + Be) = B . (B.29) 

Complex Grassmann numbers can be built out of real and imaginary parts in 
the same way as for ordinary complex numbers. It is convenient to define the 
complex conjugation in such a way that the order of a product is reversed, as is 
done in hermitian conjugation of matrices: 

(e^y = ^*e* = -e*c ■ (b.30) 

In the integral over complex Grassmann numbers 6 and 6* are treated as inde- 
pendent variables as the real and the imaginary part are independent of each 
other, so 

dMr e*e = i . (b.si) 
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Phase-space parametrization 



In this appendix the parametrization of the phase space for 2—^2 and 2 — > 3 
processes, as it was used for the calculations of this thesis, is presented. It 
is the same parametrization which is also used in FormCalc [11, 12, 13]. The 
parametrization is performed in the center-of-mass system of the two incoming 
particles, which define the beam axis and carry a center-of-mass energy of ^/s. 
For each final-state particle an integral over its three-momentum k occurs in the 
calculation of integrated cross sections. The energy fc^ of the particle is fixed by 



the on-shell condition k^ = y | A;^ -|- m^, where m denotes the mass of the particle. 
Four of these integrals are eliminated by global energy-momentum conservation. 
In the following sections the parametrizations of the two- and three-particle phase 
space are shown. 

C.l Two-particle phase space 

With two particles in the final state, labeled by the subscripts 3 and 4 in the 
following, the phase-space integral can be written in terms of two angles. They 
are the azimuth angle (p and the polar angle 6 with respect to the beam axis. 
Because of rotational invariance around the beam axis the integration over (p is 
trivial and amounts to a factor of 27r. So the integral over the two-particle phase 
space is given by 

where 

ik-^r = Ikif = i^ ^-^ (C.2) 

4s 

denotes the squared absolute value of the three-momentum of the final-state 
particles, ms and m,4 are their respective masses, and y/s specifies the center-of- 
mass energy of the incoming particles. 
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Figure CI: Graphical representation of the variables used in the parametrization 
of the 2^3 phase space. The figure is taken from ref. [13]. 



C.2 Three-particle phase space 

For the three-particle phase space, where the outgoing particles are labeled by 
the indices 3, 4 and 5, five independent integration variables remain after global 
energy-momentum conservation has been applied. They are the energies fcg and 
/cg, the azimuth angle (p and the polar angle 6 of the fifth particle with respect to 
the beam axis, and the angle fj which rotates particle 3 out of the plane defined 
by particle 5 and the beam axis. A graphical representation of the angles is given 
in Fig. C.l. 

The four-momenta of the outgoing particles have the following explicit form 



h = (^3' l^sles) 
h = ikl, \h\e^) 



ki = {^/s - /cg - k^, -h - h 



with 
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CeCfjS^ + SgCj: 



es 




(C.3) 



(C.4) 



The angle 9, which is also plotted in the figure, is defined over 



ce 



{^s-kl-k 



0N2 



m^ 



1^3 



Ifcfi 



2|fellA;fi 



(C.5) 



C.2. Three-particle phase space 
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rrii again denotes the mass of the respective particle i and ^/s is the center-of-mass 
energy of the initial-state particles. Due to axial symmetry the trivial integration 
over (j) can be carried out immediately and yields a factor of 27r. 

Then the parametrization of the three-particle phase space takes the following 
form 



dr. 



(fcO)" 



dkl 



ms 



(fcO)" 
(fcO)-i 



< 



1 /•2pi 

dc0 / di7 

1 Jo 



647r3 



(C.6) 



where the integration limits are given by 



o^max Vs {ms + m^f - m\ 



(^0)n.ax ^ ^ 



2x/i 



(C.7) 



and 



(^3) 



O^inaXjinin 



1 

27 



0"(r + raj^raS) ± \kr\\KT — 77i^)(r — my) 



(Ci 



usmg 



a 



V~s-kl 



T = a 



\k,? 



m± = m^ ± m^ 



(C.9) 



Appendix D 
Loop Integrals 



When calculating quantum corrections to physical processes Feynman diagrams 
appear which contain loops. The rules for evaluating Feynman diagrams state 
that for every closed loop an integral over the loop momentum appears in the 
expression for the amplitude. In this thesis one-loop corrections are calculated 
and thus we are concerned only with one-loop integrals. The general one-loop 
integral which corresponds to the general N-point one-loop diagram depicted in 
Fig. D.l is given by 
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[q^ - mj] [{q + ki) - ml] . . . [{q + kN-i) - m%] 

(D.l) 
The notation and conventions used in this chapter correspond to that used in 
ref. [15] . The loop momentum is denoted by q, the momenta Pi are the momenta of 
the external propagators and the momenta ki, which appear in the denominator, 
are related to the former via 



Pi = h, 


P2 = k2- ki, 


Pn = ^N 


ki = pi, 


k2=Pl+P2, 


kN =Pi 
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Figure D.l: General one-loop diagram. The arrows denote the momentum flow. 
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It is assumed implicitly that all propagators in the denominator have an infinitely 
small positive imaginary part e' which is set to zero only after the integration has 
been performed. Explicitly this is achieved by replacing the masses mf with 
mj — ie' everywhere. 

The expression given in eq. (D.l) is valid for both dimensional regular izat ion 
and dimensional reduction (see chapter 4.1). Z^ = 4 — 2e is the dimension of the 
integral, where e is a small positive number which will be sent to zero at the end 
of the calculation. The regularization parameter /i has the dimension of a mass 
and is introduced to keep the dimension of the whole expression fixed when going 
from 4 to D dimensions. 

Following ref. [149] the loop integrals are denoted by capital letters in ascend- 
ing order, resulting in 

T^ = A one-point loop integral, 

T^ = B two-point loop integral, 

T^ = C three-point loop integral, 

T'^ = D four-point loop integral, 

T^ = E five-point loop integral, .... 

Scalar integrals, which do not have an index //«, are denoted with an index 0, e.g. 
Ao. 

The scalar one-point integral Aq can be calculated analytically and reads 

Ao{m)=m''(/^-\n'\ + l\+0{e) (D.3) 



where 



A = - -7£;-hln47r (D.4) 



contains the divergent part of the loop integral with the Euler-Mascheroni con- 
stant 

7^-0.577215664901532... . (D.5) 

This one and expressions for the two-, three- and four-point loop integrals were 
first given in ref. [149] and further improved later [150, 151]. Scalar loop integrals 
with five and more internal propagators can be reduced to four-point ones [78, 
152]. 

The tensor integrals can be decomposed into linear combinations of Lorentz- 
covariant tensors [153]. They consist of a basis which is formed of linearly inde- 
pendent momenta and the metric tensor g^^ , and coefficients which are Lorentz 
scalars. This decomposition is not unique. Here the momenta are chosen as the 
momenta ki appearing in the denominator. In this basis the coefficient functions 
are totally symmetric in their indices. The decomposition for tensors up to rank 
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four reads 

AT-l 



a=l 

N-1 

a,b=l 

N-1 

A'iM2M3 ~ / y \9tJ.iH2^afj,i ~T~ 9^J.2^J.3^alJ.l ' 9^la^J.l^afl2) -'ooa 

a=l 

+ 2_^ kafj,ikbfj,2kcii3T^bc (D-8) 

a,fe,c=l 
A'lM2M3M4 ~ \9lJ.ltl2 9tJ.3IJ,4 ' 9/11(13 9 tl2IJ.4 ' 9/11(14 9(12113) "^ 0000 

Af-1 
~r / J \9/iili2^a/i^'^b/i4 ~T~ 9liili3'^a.fi2'^hfi4 ~T~ 9/ii/i4'^a/i2^b/iz 



a,b=l 

'9/i2/i3'^<ifiii^b/i4 + 9/12/14'^a/ii'^bfii, + 9/i3/i4'^afii'^bfi2) -'-OOab 
N-l 

+ _^ kafj_-^kbfj_2kcfj.^kdf^^T^f,cd ■ (D-9) 

Furthermore the coefficient functions of the tensorial loop integrals can be 
written as functions of the scalar integrals [153]. This is known as Passarino- 
Veltman reduction scheme. A complete set of equations for reducing loop inte- 
grals up to point rank A^ = 4 and up to tensor dimension A^ + 1 can be found for 
example in ref. [152]. These are all loop integrals which can appear in processes 
with up to four external legs. 

In the reduction the inverse of the Gram matrix occurs. The Gram matrix Z 
is a matrix which is built up from the momenta ki by Zij = 2kikj. This matrix 
can become singular. For up to four-point loop integrals this happens only at 
the borders of phase space. Care has to be taken when calculating phase space 
points close to the edges, as the computation can become numerically unstable. 
A technique to improve stability will be presented in appendix E.l. 

For loop integrals with five and more internal propagators the Gram matrix 
can become singular also at points inside the phase space. However the loop mo- 
menta are four-component Lorentz vectors, so a linear independent combination 
of four of them spans the whole Minkowski space. This allows one to eliminate 
the inverse Gram matrix. The reduction is done in such a way that A^-point loop 
integrals are reduced to a combination of (A^ — l)-point loop integrals with the 
tensor rank increased by one [152]. Recently a slightly different scheme has been 
found which even reduces the tensor rank by one in the decomposition [154]. 

For phase space points close to those where the Gram matrix becomes singular 
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also expansions around vanishing determinant of the Gram matrix or methods 
which use numerical integration can be applied. An overview of possible tech- 
niques with explicit formulae for loop functions with point rank up to six was 
given in ref. [155]. 

For the numerical evaluation of the loop integrals the Loop Tools package [11, 
14], which is based on FF [156], was used. In this package the stability of calculat- 
ing the Passarino-Veltman reduction was improved numerically with the method 
of Gaussian elimination, which will be described in appendix E.l. Addition- 
ally the five-point functions up to tensor rank four were implemented based on 
ref. [152] and the scalar four-point function amended according to ref. [151] so it 
is valid for all cases. The numerical results of the code were compared to those 
of an independent code from Dittmaier [157] and very good agreement could be 
found. Moreover a Passarino-Veltman reduction of the five-point tensor integrals 
was implemented. Also here a comparison yielded excellent agreement except for 
points very close to the edges of phase space where the decomposition algorithm 
is known to become numerically unstable. 

The explicit formulae for this decomposition are given below. To shorten the 
notation some abbreviations are introduced: 

S,=l-S., = l' ^'''' = \ ^, = l' ^'''^' . (D.IO) 

I 1 for z 7^ j I ^ ~ -'- ^'^^ ^ > J 

A number in brackets behind the loop integral denotes that the term in the 
denominator with a mass with this index is left out from the integrand, e.g. the 
normal scalar four-point integral 

Uq — 



m"^ 



D, 



d^g 



[g2 - ml] [{q + kif - mf\ [{q + fca)^ - rnf\ [{q + k^f - mj] 

(D.ll) 



where the second propagator is left out, becomes 



D„(2)=^.2^fi% 1 



ivr2 J [g2 _ ^2] ^(^q ^ ^^)2 _ ^2j ^^^ ^ ^^)2 _ ^2j 



(D.12) 



The resulting integral is a loop integral with the point-rank reduced by one, 
as one can see in the example above, which is a three-point loop integral. For 
integrals where the first propagator is eliminated the integration momentum must 
be shifted by g — > g — /c2, so that the standard form eq. (D.l) is again obtained. 
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With this the tensorial coefficients of the five-point loop integrals are: 

^.=E(^'*'),;'si 
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+ 2(Z(^))-;(Doo(n+l)-Doo(l)) 
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+ 2 (Z(^))^/ (Doo,„ {n + 1) 4, - D,,, (1)) 
+ 2 (^('));7 (/^oofc„ (n + 1) 4fc - Dook (1)) 

and the Gram matrix 

[kiki kik2 kik^ kik4\ 
kik2 k2k2 k2k'i k2kA 

klks ^2^3 ^^3^3 ^3^4 
\k1k4 ^2^4 ^3^4 k^k^^J 

All coefficients which are multiplied by the metric tensor, i.e. have a 00 in the 
index, vanish identically 

Eqq = Eooi = -Eqooo = Eooij = 0. (D.23) 

The decomposition into coefficients multiplied by the metric tensor and such 
multiplied by the momenta is not unique as the four linearly independent mo- 
menta ki span the whole Minkowski space and are related to the metric tensor 

by 



(D.21) 
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The formulae given in ref. [152] use a different decomposition in order to avoid 
inverse Gram matrices. This leads to a non-vanishing £'00 which is compensated 
by different Eij and the same happens for higher-tensor ranks. 

For a numerical comparison the two expressions must be transformed into 
each other by exploiting the relation eq. (D.24): 



E^=ES + 2{Z^%ES, (D.25) 

ESI =ES, + 2 {Z^%' Eg, + 2 (Z(^))-; Eg,, + 2 (Z^^);^ E,% (D.26) 

E5Z =E^u + 2 [Z^^^T^ <oH + 2 [Z^^^y^ <o.. + 2 [Z^%^ <o., 
+ 2 {Z^%1 Eg, + 2 (Z(^)):/ <, + 2 (Z(^))^/ Eo-,. 

+ 4 Uz^y^ {Z^^r: + (^^^^).' (^^^0;/ + izn7 iz'%^) E'o 



^0000 



(D.27) 

In the equations above the expressions from [152] are denoted by a superscript 
D and the ones from the Passarino-Veltman reduction by a superscript PV. 



Appendix E 
Numerical Methods 



The numerical calculation of cross sections is only possible with the aid of com- 
puter programs. Computers can do floating point operations only with a finite 
precision so rounding errors occur inevitably in many steps of the program. Ex- 
pressions which are still valid analytically might give a numerical result which is 
utter nonsense. Therefore not only the analytical correctness must be checked 
when implementing algorithms, but also that the code is numerically stable. 



E.l Gaussian Elimination 

When calculating one-loop integrals, Gram matrices occur which contain scalar 
products of the momenta, as was shown in chapter D . For loop integrals with up 
to four external legs the inverse of the Gram matrix has to be calculated. This 
matrix can become singular at the edges of phase space, i.e. for forward scattering 
or at the production threshold. Already close to the edge naive matrix inversion 
can become unstable as will be shown below. 

As is often the case in problems where the inverse of a matrix appears in the 
analytical expression the inverse of the Gram matrix itself actually is not needed 
in the Passarino-Veltman reduction scheme. This can be seen easily for example 
from the expression for the five-point loop integrals given in chapter D. Left- 
multiplying eqs. (D.13)-(D.16) with the Gram matrix Z^^^ reduces the problem 
to the problem of solving a system of linear equations of the form 

Ax = b, (E.l) 

where A is an n x n matrix and x and b are vectors of dimension n. A and b are 
input parameters and x is the solution vector of the system of equations. We are 
going to consider the most general case here, that all components of the matrix 
and the vectors can be complex numbers and 6 is a generic vector. The physical 
case of real matrices A then follows directly from this as a special case. 
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Appendix E. Numerical Methods 



Let us first consider the analytic way of multiplying the input vector with the 
inverse matrix. Then eq. (E.l) can be written as a function / 
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A^'b 



(E.2) 



Calculating the solution vector x in this manner however leads to numerical 
instabilities. Let us first decompose / into the two partial steps 



g: A 
h: A-^ 
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(E.3) 
(E.4) 



and assume the ideal case that both partial steps can be calculated in a numer- 
ically stable way. Then a stability analysis [158] yields, that the error on x is 
proportional to k{A), the condition of A. The condition is defined as 
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where || ■ || denotes a matrix norm. In a geometrical interpretation this is the 
distance of A to a singular matrix for which eq. (E.l) has no longer a unique 
solution. 

The numerical evaluation of mathematical expressions was done in double 
precision in this thesis, i.e. floating point arithmetics with double precision as 
defined in ref . [159] . These numbers offer about sixteen valid digits, so that for a 
condition k{A) > 10^^ the matrix cannot be distinguished any longer numerically 
from a singular matrix. One also says that the matrix is numerically singular. 
When calculating one-loop integrals this case occurs at the edges of phase space. 

Even earlier the error Ax of the solution vector x increases, 
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and the result becomes inaccurate. Such a behavior could indeed be observed 
while doing the numerical evaluation of the vector boson fusion processes (chap- 
ter 7). 

To avoid this problem, one decomposes the matrix A into a unipotent, i.e. 
whose diagonal contains 1, lower triangular matrix L and a non-singular upper 
triangular matrix R 
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* here is a place holder for an arbitrary complex number. The algorithm partitions 
the matrices such that 



V Aj~\w Lj [o R^ 



(E.8) 



where 

aeC; u,v,we C"-^ A„L,,R,e C^'-^^^^""!) (E.9) 

so 

w = - (E.IO) 

a 

L^R^=A^-wu^. (E.ll) 

L^, and R^ are again a unipotent lower and non-singular upper triangular matrix, 
respectively. In the next step of the iteration the matrix L^R^ has to be par- 
titioned in this way. After applying this procedure recursively one obtains the 
desired decomposition. 

If a happens to be zero, the algorithm breaks down. Therefore, not A is 
decomposed, but PA, where P is an n x n permutation matrix and chosen in such 
a way, that in every iteration step the first row and the row whose first column 
contains the largest element by absolute value are swapped. This method is called 
partial pivoting. It can be shown that with partial pivoting every non-singular 
matrix A can be decomposed into L and R and the decomposition is unique. The 
system of linear equations now has the form 

LRx = p-^b. (E.12) 

To solve it one first solves the system 

Ly = P~^b (E.13) 

with an auxiliary vector y. Because of the triangular structure of L this can be 
done very easily by a recursive forward substitution 

y, = {p-^b),-J2y, 1 = 1 ^n . (E.14) 

The arrow denotes the order in which the yi must be calculated. Finally, 

Rx = y (E.15) 

is solved by backward substitution and one obtains the solution vector x 

n 

Xi = Hi — y^ Xj i = n ^>- I . (E.16) 

j=i+i 
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The algorithm presented here is known in the hterature by the name of Gaussian 
ehmination with partial pivoting [158, 160]. 

An error analysis yields that the error on x is determined by 

< ||L|| ■ IIL^^II =k(L) (E.18) 

If one chooses the norm as the maximum norm, the inequality \Lij\ < 1 holds 
because L is unipotent and partial pivoting was used. Therefore 

\\L\\oo<n (E.19) 

||^"'||oo<2"-^ [161], (E.20) 

so 

fi;<n-2""i . (E.21) 

When calculating loop integrals only matrices of dimension n < A occur, so the 
calculation of x in this manner is absolutely stable. 



Appendix F 

Manual of the HadCalc Program 



For the calculation of hadronic cross sections a computer code, called HadCalc, 
was written (see chapter 5.5). In this appendix the manual of the program is 
presented. 

F.l Prerequisites and Compilation 

F.1.1 Prerequisites 

The following programs are required for compiling and running HadCalc and 
must be installed: 

• a Fortran compiler compliant with the Fortran?? standard, 

• a C compiler conforming to ANSI-C, 

• GNU make, 

• FormCalc 4 [11], 

• one of the two following packages that include sets of parton distribution 
functions from various groups 

- PDFLIB (CERN Computer Program Library entry W5051) [92], or 

- LHAPDF 



Additionally, support for the following two programs is integrated into HadCalc 

• FeynHiggs 2.1beta or newer [162], 

• Condor workload management system for compute-intensive jobs. 
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PDG flavor code Particle 

gluon g 

1 down quark d 

2 up quark u 

3 strange quark s 

4 charm quark c 

5 bottom quark b 

6 top quark t 

-1 down anti-quark d 

-2 up anti-quark u 

-3 strange anti-quark s 

-4 charm anti-quark c 

-5 bottom anti-quark b 

-6 top anti-quark i 

Table F.l: PDG flavor codes 

F.1.2 Configuration and Compilation 

First the partonic subprocess must be generated and prepared by following the 
instructions in the FormCalc4 manual. Especially the deflnitions in process. h 
have to be updated correctly as HadCalc relies on those. It is not necessary to 
flU in correct MSSM parameters or tune integration parameters, however. 

Then the distribution flle HadCalc-0 . 5 . tar . gz should be unpacked. As next 
step change into its subdirectory and run configure from there. The following 
conflgure options are mandatory: 

— with-partonprocess=DIR This is the location of the FormCalc- 

generated partonic output. 

— with-processtype=mn By this option the processtype is flxed, 

specifled by the number of incoming par- 
ticles m and the number of outgoing par- 
ticles n. Note that m and n form a sin- 
gle number, i.e. for a 2 ^ 2 process 
one would write — with-processtype=22. 
Currently, 2 ^ 1, 2 ^ 2 and 2 ^ 3 is 
implemented and can be entered here. 

— with-partonl='i The type of the flrst parton is specifled 

by i, given as the PDG flavor code [17] 
(see table F.l). 

— with-parton2=z Similarly, this is the PDG flavor code for 

the second parton. 



F.J. Prerequisites and Compilation 
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Additionally the following options are recognized by configure and enable optional 

features: 
— enable- antiprotonl 



— enable- antiproton2 
— with-condor[=DIR] 



— witli-feynliiggs[=DIR] 



Hadron 1 is an anti-proton instead of a 
proton. 

Hadron 2 is an anti-proton instead of a 
proton. 

Link the final code with the Condor 
workload management system libraries. 
If the binary is not in the standard search 
path of your shell, its location can be 
specified with the optional DIR argu- 
ment. 

Link the final code with the Feyn- 
Higgs library. This is mandatory if the 
FormCalc option to compute the MSSM- 
Higgs sector via FeynHiggs is chosen. 
The optional DIR specifies the location 
of the FeynHiggs library libFH.a, if it is 
not in the standard search path of the 
compiler. 

If LoopTools is not in the standard 
search path of the compiler, its location 
can be specified here. 
Use LHAPDF for the parton distribution 
functions. If the LHAPDF library is not 
in the standard search path, its location 
can be given by the optional DIR argu- 
ment. The PDF data is assumed to be 
found at the same place. 
Use PDFlib for the parton distribution 
functions. If the PDFlib library is not in 
the standard search path and the CERN- 
lib environment variables $CERN and 
$CERN_LEVEL are not set, the DIR 
argument designates where it can be 
found. 

Only one of the last two options can be given on the command line. If neither 
— with-lhapdf nor — with-pdflib was given, configure first tries to find LHAPDF 
and, if this fails, probes the existence of PDFlib. 

After having run configure, a call to make compiles the program. When 
it successfully finishes, a binary called HadCalc has been created in the current 
path. 



-with-looptools=DIR 
-with-lhapdf[=DIR] 



-with-pdflib [=DIR] 
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F.2 Running the program 

The program is simply started by running ./HadCalc. It will then present a 
menu which allows one to tune various settings and start the calculation of cross 
sections. The following subsections describe the possible settings in detail. An 
item is chosen by typing the number shown in brackets before the item and 
pressing "Enter". In every menu "(0)" exits the submenu or, for the top level 
menu, quits the program. Invalid input is ignored and an error message is written 
on the screen. 



F.2.1 Physics parameters 

This submenu sets the parameters of the MSSM and related things and is divided 
in three further submenus. 



MSSM parameters 

Here all values which correspond to parameters of the MSSM are set. 

First let us look at menu item 16. This decides whether the program should 
use a common mass Msusy in the sfermion sector, or if individual values for 
the left-handed squarks and sleptons and the right-handed sups, sdowns and 
selectrons are allowed. Depending on this flag either the MS* variables cannot 
be set (because they are fixed at Msusy) ot Msusy itself cannot be set (because 
it is irrevelant and not used in the computation). When choosing a common 
SUSY mass scale, the settings in the MS* variables are retained and restored 
when deselecting this option. 

All other parameter settings can be in two states. They can either have a 
fixed setting, then this value is used for all calculations. Or their value can be 
running. In this case a lower and upper bound and the number of intermediate 
intervals must be chosen. Then the computation of the cross section is done 
("intermediate intervals" + 1) times, with the value of the parameter increasing 
from lower bound to upper bound^. The distance between two values is equal for 
the setting "linear" and exponentially increasing for the setting "logarithmic", 
i.e. the values are closer at the lower bound and they have equal distance again 
when plotting them on a logarithmic scale. A behavior vice versa with values 
closer at the upper bound can be easily achieved by exchanging lower and upper 
bound. If more than one parameter is chosen to be running, the iteration loops 
are nested, with the first parameter varying fastest. 



^Despite its name, the lower bound can be mathematically larger than the upper bound; 
then the value of the parameter is decreasing. 
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Kinematic parameters 

In this menu all parameters are set which are related to kinematic variables of 
the process. 

The underlying parameter of items 3 and 4 depends on the type of process. 
For processes with two particles in the final state, this is the angle 6 between the 
two outgoing particles, for those with three final particles, it denotes the energy 
^5 of the fifth particle, which is the third final-state particle. The menu items 8, 
12, 14 and 15, which refer to the fifth particle, are ignored for 2 — > 2 processes 
and cannot be changed. 

The settings of the parameters are possible in the same way as already de- 
scribed in the previous item. 

Scale parameters 

This menu sets the renormalization and factorization scale of the process in the 
same way as described above. A negative number for the renormalization scale 
has a special meaning. Then the sum of the masses of the final-state particles 
is taken, multiplied with the absolute value of the setting, and this number is 
taken as the renormalization scale. Additionally it can be chosen that both 
renormalization and factorization scale are always set to the same value. 

Show ModelDigest (FormCalc) 

Finally this choice invokes FormCalc's ModelDigest function, which takes the pa- 
rameters as an input and calculates the physical masses of the particles. Thereby 
it applies lower bounds on the masses established by experiment and refuses the 
calculation if these bounds are violated. The calculated cross section will also 
be zero in that case. The FormCalc manual contains a more detailed explana- 
tion of this function. There it is also described how the check for the violation 
of experimental bounds can be switched off by flipping a switch in FormCalc's 
process.h. 

F.2.2 PDF parameters 

The set used for the calculation of the parton distribution functions is chosen by 
this submenu. The layout and choices presented depends on whether LHAPDF 
or PDFlib is used. For PDFlib three numbers must be entered. The first de- 
notes the type of parton distribution functions and is 1 for proton PDFs. The 
second number specifies the respective group which has performed the fit to the 
experimental data and the third number chooses a specific PDF set. When using 
LHAPDF a string must be entered that directly specifies the filename of the PDF 
set in the LHAPDF subdirectory. 



144 Appendix F. Manual of the HadCalc Program 

F.2.3 Integration parameters 

This submenu chooses the integration routine and sets its parameters. Currently 
there are six integration routines available: 

GAUSSKR One-dimensional Gauss-Kronrod algorithm 

GAUSSAD One-dimensional adaptive Gauss algorithm 

DCUHRE Multi-dimensional adaptive Gauss algorithm 

VEGAS Monte Carlo integration algorithm 

SUAVE Subregion adaptive Monte Carlo integration algorithm 

DIVONNE Monte Carlo integration via stratified sampling 

The last four algorithms are part of the CUBA library [163]. In the following 
only a short description of the possible parameter settings is given. The technical 
details of these algorithms and the precise impact of the variables are described 
in the CUBA manual and shall not be repeated here. 

The GAUSSKR and GAUSSAD algorithms can only handle one- dimensional 
integrals. If multi-dimensional integrals are attempted to be computed, VEGAS 
is used as a fallback. In contrast the DCUHRE and DIVONNE algorithms can- 
not handle one-dimensional integrals. There the GAUSSKR algorithm is used 
instead. In both cases a warning is printed on the screen. 

All integration routines share these two variables: 

• relative error: the desired relative error 

• absolute error: the desired absolute error 

Additionally, the following variables are available for one or more of the routines. 
Which ones these are is denoted in square brackets after the entry. 

• maximum 7^ of points: the maximum number of function evaluations used 
[GAUSSAD, DCUHRE, VEGAS, SUAVE, DIVONNE] 

• 7^ of points for starting: the initial number of points per iteration [VEGAS] 

• increase in ^ of points: the number of points the previous number is incre- 
mented for the next iteration [VEGAS] 

• # of points for subdivision: the number of points used to sample a subdi- 
vision [SUAVE] 

• flatness ^ for splits: the type of norm used to compute the fluctuation of 
a sample [SUAVE] 

• # of passes: the number of passes after which the partitioning phase ter- 
minates [DIVONNE] 
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• key 1: determination of sampling in the partitioning phase [DIVONNE] 

• key 2: determination of samphng in the final integration phase [DIVONNE] 

• key 3: sets the strategy for the refinement phase [DIVONNE] 

• maximum x^ for subregion: the maximum x^ value a single subregion is 
allowed to have in the final integration phase [DIVONNE] 

• minimum deviation for split: a bound which determines whether it is worth- 
while to further examine a region that failed the x^ test [DIVONNE] 

F.2.4 Amplitude switches 

This submenu sets the type of diagrams used for the computation and how the 
cuts should be applied. The value of the cuts is set in the parameter section and 
was already described there. 

The first entry decides whether the tree-level and the one-loop result shall 
be computed in one go or only one of them. Possible choices are "Tree only", 
"Tree+Loop" and "Loop only". Which way is better depends on the concrete 
circumstances and features of the problem. Computing both at the same time 
might save computation time, but the integration routine has to optimize its 
choices for both at the same time, which might lead to sub-optimal performance. 
On the other hand it is not too likely that there are problematic regions in the 
tree-level cross section which are no longer present in the one-loop computation, 
so normally this procedure gives satisfactory results. If only one cross section is 
computed, the value of the other one is set to zero. 

The remaining entries decide if and how the cuts on rapidity, transverse mo- 
mentum and jet separation should be applied. It is either possible to have the 
particle, or a pair of particles in case of the jet separation, fulfill a cut, violate it 
or ignore the cut altogether. Since HadCalc relies on FormCalc for the partonic 
process and implementation details, for the cuts for particle three in the 2 — »• 2 
case and particle five in the 2-^3 case it cannot be chosen that the rapidity and 
transverse momentum cut is violated, but they always have to be fulfilled. They 
can, however, be switched off by setting the relevant entry in the parameters 
section to zero. 

F.2.5 Input /Output options 

This submenu allows one to read in a set of parameters from a file and specify 
where and how to write the calculation output. 

To read in a set of parameters a parameter specification must have been 
written into a file and this filename then has to be entered here. All possible 
variables which can be set in such a file are given in section F.3. There are three 
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basic types of variables. Those which specify a parameter can either take four 
comma-separated values that are the lower and upper bound, the behavior with 
respect to increments, i.e. linear or logarithmic, and the number of intermediate 
intervals, or a single number denoting its fixed value. The ones of type boolean 
turn on a certain switch and take no arguments. All remaining ones take a single 
argument and the variable is assigned to this parameter. 

In the following also a formal definition of the parsing rules is given: 

• The file is read line by line. 

• White space at the beginning of a line is ignored. 

• Empty lines are ignored. 

• Lines starting with the character "7^" (after optional white space) are com- 
ment lines and ignored. 

• The first token which is separated by white space from the rest of the line 
is extracted. This token has to be a token from the list of valid tokens in 
section F.3. 

• If the token type is boolean, its associated parameter is set. 

• If the token type is integer, an attempt to read an integer value is made 
and if it succeeds, this is assigned to the associated parameter. 

• If the token type is double, an attempt to read a double precision floating- 
point number is made and if it succeeds, this is assigned to the associated 
parameter. 

• If the token type is string, the second token is assigned to the associated 
parameter. 

• If the token type is parameter, the following actions happen: 

— An attempt to read four comma-separated double precision floating- 
point numbers is made. 

— If this attempt succeeds, the four numbers are assigned to lower bound, 
upper bound, log and number of intermediate intervals of the param- 
eter, log means linear increase if this variable is zero and exponential 
one otherwise. 

— If this does not succeed, an attempt to read a single double precision 
floating point number is made. 

— If this succeeds, this number is the constant value of the parameter. 

— If this also does not succeed, the line is flagged as not parsable. 
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• For lines not parsable by the rules above a warning message is printed and 
their content is ignored. 

Furthermore some integration routines offer the possibility to write interme- 
diate results or progress report to the screen. This is turned on with Verbose 
integration output. For hadronic cross sections this also enables writing PDFlib 
statistics on the screen at the end of the calculation. 

Finally one can choose whether the calculation results will be written to the 
screen or into a file. In the latter case the variable outputstring describes which 
elements should be written to the output file. The form of this variable together 
with the valid tokens is described in section F.4. The output-file format starts 
with a "#" -quoted header with a file identification and the content of output- 
string. Then, each on a line by itself, for every scanned parameter point the 
values defined in outputstring are written, separated from each other by a space. 



F.2.6 Amplitude calculation 

This submenu finally allows one to choose the cross section one wants to compute 
and does the calculation. During the following integration the process may be 
interrupted with "Ctrl-C", after which it aborts the current calculation and jumps 
back into the main menu. Due to restrictions imposed by Condor this feature 
is not available if HadCalc was configured with the option — with-condor. Here 
pressing "Ctrl-C" aborts the whole HadCalc program. 



F.3 Allowed tokens in input files 

The following list shows all token names that may appear in an input file together 
with its associated type. The tokens are not case-sensitive. Thereby parameter 
means that the variable can either be followed by four comma-separated values 
that denote the lower and upper bound, whether the increase is linear or loga- 
rithmic, and the number of intermediate intervals, or a single number that is the 
fixed value of this parameter, boolean means that a specific behavior is switched 
on. There is a corresponding separate token that switches the same behavior off 
again, double and integer tokens take a single double-precision or integer value, 
respectively, as input, string assigns the remainder of the line to the parameter. 
Finally preselected takes special values as an argument. The possible choices 
for each of these ones were discussed during the description of the menus given 
in section F.2. Any settings referring to particle 5 are relevant only for 2 — i> 3 
processes and will be silently ignored otherwise. 
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token 


: type 


description 


MAO 


: parameter 


mass of the CP-odd Higgs boson 


TB 


: parameter 


ratio of the Higgs vacuum expecta- 
tion values 


MUE 


: parameter 


/i parameter in the Higgs sector 


MSusy 


: parameter 


common SUSY mass scale 


MSQ 


: parameter 


mass parameter of the left-handed 
squarks 


MSU 


: parameter 


mass parameter of the right-handed 
sup-like squarks 


MSD 


: parameter 


mass parameter of the right-handed 
sdown-like squarks 


MSL 


: parameter 


mass parameter of the left-handed 
sleptons 


MSE 


: parameter 


mass parameter of the right-handed 
selectron-like sleptons 


A_t 


: parameter 


trilinear coupling of the sup-like 

squarks 


A_b 


: parameter 


trilinear coupling of the sdown-like 
squarks 


A_tau 


: parameter 


trilinear coupling of the selectron- 
like sleptons 


M_l 


: parameter 


U(l)y gaugino mass 


M_2 


: parameter 


SU(2)l gaugino mass 


MGl 


: parameter 


SU(3)c gaugino mass 


SQRTS 


: parameter 


square root of the hadronic center 
of mass energy 


SQRTSHAT 


: parameter 


square root of the partonic center of 
mass energy 


THETA 2 


: parameter 


angle between the two outgoing 
particles (in degrees) ^ 


THETACUT ^ 


: parameter 


cut on the angle between the two 
outgoing particles (in degrees) 


K50 3 


: parameter 


energy of the third outgoing particle 


K50CUT 


: parameter 


cut on the energy of the third out- 
going particle 


PTRANS 3 


: parameter 


transverse momentum 


PTRANS3CUT 


: parameter 


cut on the transverse momentum of 
particle 3 


PTRANS4CUT 


: parameter 


cut on the transverse momentum of 
particle 4 



^only for 2 — s- 2 processes 
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token 



type 



description 



PTRANS5CUT : parameter 

RAPID ^ : parameter 

RAPID3CUT : parameter 

RAPID4CUT : parameter 

RAPID5CUT : parameter 

DELTAR34CUT : parameter 

DELTAR35CUT : parameter 

DELTAR45CUT : parameter 

RENSCALE : parameter 

FACTSCALE : parameter 

CommonSUSYMassScale : boolean 

NoCommonSUSYMassScale : boolean 

CommonRenFactScale : boolean 

NoCommonRenPactScale : boolean 



AMPLITUDE 
Ptrans3>cut 

Ptrans3<cut 



: preselected 
: boolean 

: boolean 



Ptrans3nocut 


: boolean 


Rapid3>cut 


: boolean 


Rapid3<cut 


: boolean 


Rapid3nocut 


: boolean 


Ptrans4>cut 


: boolean 



cut on the transverse momentum of 
particle 5 
rapidity 

cut on the rapidity of particle 3 
cut on the rapidity of particle 4 
cut on the rapidity of particle 5 
cut on the distance between parti- 
cles 3 and 4 

cut on the distance between parti- 
cles 3 and 5 

cut on the distance between parti- 
cles 4 and 5 
renormalization scale 
factorization scale 
choose a common SUSY mass scale 
do not choose a common SUSY 
mass scale 

choose a common remormalization 
and factorization scale 
do not choose a common remormal- 
ization and factorization scale 
choose which amplitude (s) to calcu- 
late 

require the transverse momentum 
of particle 3 to be larger than the 
cut 

require the transverse momentum 
of particle 3 to be smaller than the 
cut 

disable cut on the transverse mo- 
mentum of particle 3 
require the rapidity of particle 3 to 
be larger than the cut 
require the rapidity of particle 3 to 
be smaller than the cut 
disable cut on the rapidity of parti- 
cle 3 

require the transverse momentum 
of particle 4 to be larger than the 
cut 



^only for differential cross sections 
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token 



type 



description 



Ptrans4<cut 


: boolean 


Ptrans4nocut 


: boolean 


Rapid4>cut 


: boolean 


Rapid4<cut 


: boolean 


Rapid4nocut 


: boolean 


DeltaR34>cut 


: boolean 


DeltaR34<cut 


: boolean 


DeltaR34nocut 


: boolean 


DeltaR35>cut 


: boolean 



DeltaR35<cut 



boolean 



DeltaR35nocut 


boolean 


DeltaR45>cut 


boolean 


DeltaR45<cut 


boolean 


DeltaR45nocut 


boolean 


INTMETHOD 


preselected 


EPSABS 


double 


EPSREL 


double 


MAXPTS 


integer 


VSTARTPTS 


integer 


VINCREASE 


integer 


SNNEW 


integer 



require the transverse momentum 
of particle 4 to be smaller than the 
cut 

disable cut on the transverse mo- 
mentum of particle 4 
require the rapidity of particle 4 to 
be larger than the cut 
require the rapidity of particle 4 to 
be smaller than the cut 
disable cut on the rapidity of parti- 
cle 4 

require the jet separation between 
particles 3 and 4 to be larger than 
the cut 

require the jet separation between 
particles 3 and 4 to be smaller than 
the cut 

disable the cut on the jet separation 
between particles 3 and 4 
require the jet separation between 
particles 3 and 5 to be larger than 
the cut 

require the jet separation between 
particles 3 and 5 to be smaller than 
the cut 

disable the cut on the jet separation 
between particles 3 and 5 
require the jet separation between 
particles 4 and 5 to be larger than 
the cut 

require the jet separation between 
particles 4 and 5 to be smaller than 
the cut 

disable the cut on the jet separation 
between particles 4 and 5 
choose the integration routine 
absolute integration error 
relative integration error 
maximum number of points 
number of points for starting 
increase in number of points 
number of points for subdivision 
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token 



type 



description 



SFLATNESS 


integer 


MAXDPASS 


integer 


DKEYl 


integer 


DKEY2 


integer 


DKEY3 


integer 


DBORDER 


double 


MAXDCHISq 


double 


MINDDEV 


double 


VERBOSITY 


integer 


PDFTYPE 


double 


PDFGROUP 


double 


PDFSET 


double 


PDFPATH 


string 


PDFNAME 


string 


ScreenOutput 


boolean 


OUTPUTFILE 


string 


OUTPUTSTRING 


string 



flatness number for splits 

number of passes in partitioning 

phase 

Divonne key 1 

Divonne key 2 

Divonne key 3 

border of the integration region 

maximum x^ for subregion 

minimum deviation for split 

verbosity of integration output 

type of the PDF [PDFhb] 

group of the PDF [PDFlib] 

set of the PDF [PDFlib] 

path where the PDF files are 

[LHAPDF] 

name of the PDF [LHAPDF] 

print output on the screen 

print output into file 

parameters to print in output (see 

section F.4) 



F.4 Allowed variable names for outputstring 

The following list shows all variable names that may appear in outputstring. The 
individual entries are separated from each other by spaces. Variables with the 
dimension of a mass are output in GeV. Note that these names are case-sensitive. 



Name 



Parameter description 



MAO 


niA, 


TB 


tan/3 


MUE 


/i 


MSusy 


"f^SUSY 


MSQ 


rriq 


MSU 


Ulfi 



MSD 

MSL 
MSE 

A_t 



: mj 

: mj 
■ me 

: A 



mass of the CP-odd Higgs boson 

ratio of the Higgs vacuum expectation values 

/i parameter in the Higgs sector 

common SUSY mass scale 

mass parameter of the left-handed squarks 

mass parameter of the right-handed sup-like 

squarks 

mass parameter of the right-handed sdown-like 

squarks 

mass parameter of the left-handed sleptons 

mass parameter of the right-handed selectron-like 

sleptons 

trilinear coupling of the sup-like squarks 
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Name 



Parameter description 



A_b 

A_tau 

M_l 

M_2 

MGl 

SQRTS ^ 

SQRTSHAT ^ 

THETA 6 

THETACUT ^ 

K50 ^ 

K50CUT ■^ 

PTRANS 

PTRANS3CUT 

PTRANS4CUT 

PTRANS5CUT 

RAPID 

RAPID3CUT 

RAPID4CUT 

RAPID5CUT 

DELTAR34CUT 

DELTAR35CUT 

DELTAR45CUT 

RENSCALE 

FACTSCALE 

MhO 

MHO 

MHpm 

MCha(l) 

MCha(2) 

MNeu(l) 

MNeu(2) 

MNeu(3) 

MNeu(4) 

MGl 

MSn(l) 



Ab 

Ar 

Ml 
M2 



■.e, 



cut 

K 

'^5 cut 
Ptrans 

Ptranscut 

4 

Ptrans cut 
Ptrans cut 

V 

Vcut 
4 

Vcut 
Vcut 

^R cut 
^R cut 
^R cut 

fip 

rriho 
niHO 
mH± 



m 



"Xi 
X2 



Xi 
X2 

X3 
rrig 



trilinear coupling of the sdown-like squarks 
trilinear coupling of the selectron-like sleptons 
U(l)y gaugino mass 
SU(2)i gaugino mass 
gluino mass 

square root of the hadronic center-of-mass energy 
square root of the partonic center-of-mass energy 
angle between the two outgoing particles (in de- 
grees) 

cut on the angle between the two outgoing particles 
(in degrees) 

energy of the third outgoing particle 
cut on the energy of the third outgoing particle 
transverse momentum 

cut on the transverse momentum of particle 3 
cut on the transverse momentum of particle 4 
cut on the transverse momentum of particle 5 
rapidity 

cut on the rapidity of particle 3 
cut on the rapidity of particle 4 
cut on the rapidity of particle 5 
cut on the distance between particles 3 and 4 
cut on the distance between particles 3 and 5 
cut on the distance between particles 4 and 5 
renormalization scale 
factorization scale 

mass of the lighter CP-even Higgs boson 
mass of the heavier CP-even Higgs boson 
mass of the charged Higgs boson 
mass of the lighter chargino 
mass of the heavier chargino 
mass of the lightest neutralino 
mass of the second-lightest neutralino 
mass of the second-heaviest neutralino 
mass of the heaviest neutralino 
mass of the gluino 
mass of the electron sneutrino 



^only relevant for the computation of hadronic cross sections 

^only relevant for the computation of partonic cross sections and differential hadronic cross 
section with respect to invariant mass 
^only for 2^2 processes 
^only for 2^3 processes 
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Name 


Parameter desci 


ription 


MSn(2) 


mi>^ 


mass 


of the muon sneutrino 


MSn(3) 


rriu. 


mass 


of the tau sneutrino 


MSl(l) 


rrie^ 


mass 


of the lighter selectron 


MSI (2) 


rufi. 


mass 


of the lighter smuon 


MS1(3) 


nif^ 


mass 


of the lighter stau 


MSL(l) 


rrie^ 


mass 


of the heavier selectron 


MSL(2) 


m-^2 


mass 


of the heavier smuon 


MSL(3) 


mf2 


mass 


of the heavier stau 


MSu(l) 


rriu^ 


mass 


of the lighter sup 


MSu(2) 


rric. 


mass 


of the lighter scharm 


MSu(3) 


m^ 


mass 


of the lighter stop 


MSU(l) 


rriu^ 


mass 


of the heavier sup 


MSU(2) 


mc2 


mass 


of the heavier scharm 


MSU(3) 


^h 


mass 


of the heavier stop 


MSd(l) 


^d. 


mass 


of the lighter sdown 


MSd(2) 


m-s^ 


mass 


of the lighter sstrange 


MSd(3) 


^l. 


mass 


of the lighter sbottom 


MSD(l) 


^d2 


mass 


of the heavier sdown 


MSD(2) 


ms2 


mass 


of the heavier sstrange 


MSD(3) 


^b2 


mass 


of the heavier sbottom 


TREE 


a' 


tree-level cross section 


LOOP 


a^ 


one-loop cross section 


TREEERR 


o-(c^o) 


integ 


ration error of the tree-level cross section 


LOOPERR 


a{ax) 


integ 


ration error of the one-loop cross section 


TREEPROB 


xHcriao)) 


probability of the integration error of the tree- level 



LOOPPROB 

NREGIONS^ 

NEVAL^ 

FAIL^ 



cross section 
• X^('^(ci)) probability of the integration error of the one- loop 
cross section 

number of regions used for integration 
number of function evaluations used for integration 
a non-zero value indicates that the desired accuracy 
could not be reached 



only relevant for some integration routines 



Bibliography 



[1] S. L. Glashow. Nucl. Phys., 22:579-588, 1961; S. Weinberg. Phys. 
Rev. Lett, 19:1264-1266, 1967; A. Salam. WEAK AND ELECTRO- 
MAGNETIC INTERACTIONS, Originally printed in *Svartliolm: Elemen- 
tary Particle Theory, Proceedings Of The Nobel Symposium Held 1968 At 
Lerum, Sweden*, Stockholm, 367-377, 1968. 

[2] S. L. Glashow, J. Iliopoulos, and L. Maiani. Phys. Rev., D2:1285-1292, 
1970. 

[3] P. W. Higgs. Phys. Lett, 12:132-133, 1964; P. W. Higgs. Phys. Rev. 
Lett, 13:508-509, 1964; P. W. Higgs. Phys. Rev., 145:1156-1163, 1966; 
F. Englert and R. Brout. Phys. Rev. Lett, 13:321-322, 1964; T. W. B. 
Kibble. Phys. Rev., 155:1554-1561, 1967. 

[4] A. Heister et al. Phys. Lett, B526:191-205, 2002. 

[5] J. Wess and B. Zumino. Nucl. Phys., B78:l, 1974. 

[6] H. E. Haber and G. L. Kane. Phys. Rept, 117:75, 1985. 

[7] M. Carena, D. Garcia, U. Nierste, and C. E. M. Wagner. Nucl. Phys., 
B577:88-120, 2000. 

[8] J. Guasch, W. HoUik, and S. Penaranda. Phys. Lett, B515:367-374, 2001. 

[9] J. Kublbeck, M. Bohm, and A. Denner. Comput Phys. Commun., 60:165- 
180, 1990; T. Hahn. Comput Phys. Commun., 140:418-431, 2001; 
T. Hahn and C. Schappacher. Com,put Phys. Com,m,un., 143:54-68, 2002. 

[10] T. Hahn. FeynArts 3.2 User's Cuide, October 2004. 

http : //www . f eynarts . de/looptools/LT2Guide . ps . gz. 

[11] T. Hahn and M. Perez- Victoria. Comput. Phys. Com,m,un., 118:153-165, 
1999. 

[12] T. Hahn. New developm,ents in FormCalc 4-1, hep-ph/0506201, 2005. 

155 



156 Bibliography 

[13] T. Hahn. FormCalc 4 User's Guide, March 2005. 

http : //www . f eynarts . de/f ormcalc/FC4Guide . ps . gz. 

[14] T. Hahn and M. Ranch. News from FormCalc and LoopTools, MPP-2006-7, 
2006. 

[15] T. Hahn. LoopTools User's Guide, June 2004. 

http : //www . f eynarts . de/looptools/LT2Guide . ps . gz. 

[16] J. A. Aguilar-Saavedra et ah Supersymmetry Parameter Analysis: SPA 
Gonvention and Project, hep-ph/0511344, 2005. 

[17] S. Eidelman et ah Phys. Lett, B592:l, 2004. 

[18] Y. Nambu. Phys. Rev. Lett, 4:380-382, 1960; J. Goldstone. Nuovo Gim., 
19:154-164, 1961; J. Goldstone, A. Salam, and S. Weinberg. Phys. Rev., 
127:965-970, 1962. 

[19] R. Barate et al. Phys. Lett, B565:61-75, 2003. 

[20] The LEP Electroweak Working Group. A combination of preliminary 
electroweak measurements and constraints on the standard model, hep- 
ex/0511027, 2005. 

[21] V. Barger, T. Han, P. Langacker, B. McElrath, and P. Zerwas. Phys. Rev., 
D67:115001, 2003; C. Grojean, G. Servant, and J. D. Wells. Phys. Rev., 
D71:036001, 2005; S. Kanemura, Y. Okada, E. Senaha, and C. P. Yuan. 
Phys. Rev., D70: 115002, 2004. 

[22] T. Plehn and M. Ranch. Phys. Rev., D72:053008, 2005. 

[23] G. W. Bennett et al. Phys. Rev. Lett, 89:101804, 2002. 

[24] M. Passera. Status of the standard model prediction of the muon g-2, hep- 
ph/0509372, 2005. 

[25] G. Bertone, D. Hooper, and J. Silk. Phys. Rept, 405:279-390, 2005. 

[26] D. N. Spergel et al. Astrophys. J. SuppL, 148:175, 2003; C. J. MacTavish 
et al. Gosmological parameters from the 2003 flight of BOOMERANG, 
astro-ph/0507503, 2005. 

[27] H. Georgi and S. L. Glashow. Phys. Rev. Lett, 32:438-441, 1974; S. Di- 
mopoulos and H. Georgi. Nucl. Phys., B193:150, 1981; N. Sakai. Zeit. 
Phys., Cll:153, 1981. 

[28] U. Amaldi, W. de Boer, and H. Furstenau. Phys. Lett., B260:44 7-455, 1991. 



Bibliography 157 

[29] H. Fritzsch and P. Minkowski. Ann. Phys., 93:193-266, 1975. 

[30] M. Drees. An introduction to super symmetry., hep-ph/9611409, 1996. 

[31] M. B. Smy. Five years of neutrino physics with Super- Kamiokande, hep- 
ex/0206016, 2002. 

[32] S. R. Coleman and J. Mandula. Phys. Rev., 159:1251-1256, 1967. 

[33] R. Haag, J. T. Lopuszanski, and M. Sohnius. Nucl. Phys., B88:257, 1975. 

[34] J. Wess and J. Bagger. Supersymmetry and Supergravitation. Princeton 
Series in Physics. Princeton University Press, 1992. 

[35] S. P. Martin. A supersymmetry primer, hep-ph/9709356, 1997. 

[36] M. F. Sohnius. Phys. Kept, 128:39-204, 1985. 

[37] P. Fayet and J. Ihopoulos. Phys. Lett, B51:461-464, 1974. 

[38] L. O'Raifeartaigh. Nucl. Phys., B96:331, 1975. 

[39] H. E. Haber. Introductory low- energy supersymmetry, hep-ph/9306207, 
1993. 

[40] H. P. Miles. Phys. Kept, 110:1, 1984. 

[41] L. Girardello and M. T. Grisaru. Nucl. Phys., B194:65, 1982. 

[42] S. L. Adler and W. A. Bardeen. Phys. Rev., 182:1517-1536, 1969; J. S. 
Bell and R. Jackiw. Nuovo Cim., A60:47-61, 1969. 

[43] J. Rosiek. Phys. Rev., D41:3464, 1990; J. Rosiek. Complete set of Feynman 
rules for the MSSM - ERRATUM, hep-ph/9511250, 1995. 

[44] L. D. Faddeev and V. N. Popov. Phys. Lett, B25:29-30, 1967. 

[45] G. R. Farrar and P. Fayet. Phys. Lett, B76:575-579, 1978. 

[46] J. R. Elhs, G. Ridolfi, and F. Zwirner. Phys. Lett, B257:83-91, 1991; 
Y. Okada, M. Yamaguchi, and T. Yanagida. Prog. Theor. Phys., 85:1-6, 
1991; H. E. Haber and R. Hempfiing. Phys. Rev. Lett, 66:1815-1818, 1991. 

[47] A. Brignole. Phys. Lett, B281:284-294, 1992. 

[48] P. H. Chankowski, S. Pokorski, and J. Rosiek. Phys. Lett, B286:307-314, 
1992; P. Chankowski, S. Pokorski, and J. Rosiek. Nucl. Phys., B423:437- 
496, 1994. 



158 Bibliography 



[49] A. Dabelstein. Nud. Phys., B456:25-56, 1995; A. Dabelstein. Z. Phys., 
C67:495-512, 1995. 

[50] G. Degrassi, S. Heineineyer, W. Hollik, P. Slavich, and G. Weiglein. Eur. 
Phys. J., C28:133-143, 2003. 

[51] R. Hempfling and A. H. Hoang. Phys. Lett, B331:99-106, 1994. 

[52] S. Heinemeyer, W. Hollik, and G. Weiglein. Phys. Rev., D58:091701, 1998; 
S. Heinemeyer, W. Hollik, and G. Weiglein. Phys. Lett, B440:296-304, 
1998. 

[53] S. Heinemeyer, W. Hollik, and G. Weiglein. Eur. Phys. J., C9:343-366, 
1999. 

[54] R.-J. Zhang. Phys. Lett, B447:89-97, 1999; J. R. Espinosa and R.-J. 
Zhang. JHEP, 03:026, 2000. 

[55] G. Degrassi, R Slavich, and F. Zwirner. Nucl. Phys., B611:403-422, 2001. 

[56] J. R. Espinosa and R.-J. Zhang. Nucl. Phys., B586:3-38, 2000. 

[57] A. Brignole, G. Degrassi, P. Slavich, and F. Zwirner. Nucl. Phys., B631:195- 
218, 2002. 

[58] A. Brignole, G. Degrassi, P. Slavich, and F. Zwirner. Nucl. Phys., B643:79- 
92, 2002. 

[59] S. Heinemeyer, W. Hollik, H. Rzehak, and G. Weiglein. Eur. Phys. J., 
C39:465-481, 2005. 

[60] A. Dedes, G. Degrassi, and P. Slavich. Nucl. Phys., B672: 144-162, 2003. 

[61] S. P. Martin. Phys. Rev., D65: 116003, 2002; S. P. Martin. Phys. Rev., 
D66:096001, 2002; S. P. Martin. Phys. Rev., D67:095012, 2003; S. P. Mar- 
tin. Phys. Rev., D68:075002, 2003; S. P. Martin. Phys. Rev., D70:016005, 
2004. 

[62] S. P. Martin. Phys. Rev., D71:016012, 2005. 

[63] S. Heinemeyer, W. Hollik, and G. Weiglein. Phys. Lett, B455: 179-191, 
1999. 

[64] N. Cabibbo and M. J. G. Veltman. Weak interactions, CERN-65-30, 1965; 
M. Kobayashi and T. Maskawa. Prog. Theor. Phys., 49:652-657, 1973. 

[65] F. Gabbiani, E. Gabrielli, A. Masiero, and L. Silvestrini. Nucl. Phys., 
B477:321-352, 1996. 



Bibliography 159 



66 

67 

68 
69 

70 
71 

72 

73 

74 

75 
76 

77 

78 
79 

80 



82 



W. Pauli and F. Villars. Rev. Mod. Phys., 21:434-444, 1949. 

C. G. BoUini and J. J. Giambiagi. Nuovo Cim., B12:20-25, 1972; J. F. 
Ashmore. Lett. Nuovo Cim., 4:289-290, 1972; G. 't Hooft and M. J. G. 
Veltman. Nucl. Phys., B44:189-213, 1972. 

W. Siegel. Phys. Lett., B84:193, 1979. 

D. M. Capper, D. R. T. Jones, and P. van Nieuwenhuizen. Nucl. Phys., 
B167:479, 1980. 

W. Siegel. Phys. Lett., B94:37, 1980. 

D. Stockinger. JHEP, 03:076, 2005. 

G. 't Hooft. Nucl. Phys., B33:173-199, 1971; G. 't Hooft. Nucl. Phys., 
B35:167-188, 1971. 

G. 't Hooft. Nucl. Phys., B61:455-468, 1973. 

W. A. Bardeen, A. J. Buras, D. W. Duke, and T. Muta. Phys. Rev., 
D18:3998, 1978; A. J. Buras. Rev. Mod. Phys., 52:199, 1980. 

G. Passarino and M. J. G. Veltman. Phys. Lett., B237:537, 1990. 

W. J. Marciano and A. Sirlin. Phys. Rev. Lett., 46:163, 1981; A. Sirlin. 
Phys. Lett., B232:123, 1989. 

M. Bohm, H. Spiesberger, and W. Hollik. Fortsch. Phys., 34:687-751, 1986; 
W. Hollik. Renormalization of the Standard Model, MPI-PH-93-21, 1993. 

A. Denner. Fortschr. Phys., 41:307-420, 1993. 

M. Frank, S. Heinemeyer, W. Hollik, and G. Weiglein. FeynHiggsl.2: Hy- 
brid MS-bar / on-shell renormalization for the CP-even Higgs boson sector 
in the MSSM, hep-ph/0202166, 2002. 

T. Fritzsche. Berechnung von Observablen zur supersymmetrischen Teilch- 
enerzeugung an Hochenergie-Collidern unter Einschlufi hoherer Ordnungen. 
PhD thesis, Universitat Karlsruhe, 2005. 

S. Alam, K. Hagiwara, S. Matsunioto, K. Hagiwara, and S. Matsumoto. 
Phys. Rev., D55:1307-1315, 1997; Z. Sullivan. Phys. Rev., D56:451-457, 
1997. 

S. Berge. Supersymmetrische Einschleifen- QCD- Strahlungskorrekturen zur 
Top-Quark-Produktion am Fermilab Tevatron. Diplom thesis, Universitat 
Karlsruhe, 1999. 



160 Bibliography 

[83] M. Carena, S. Mrenna, and C. E. M. Wagner. Phys. Rev., D60:075010, 1999; 
M. Carena, S. Mrenna, and C. E. M. Wagner. Phys. Rev., D62:055008, 2000. 



84 
85 

86 



89 

90 
91 

92 

93 

94 
95 



96 



97 



A. Dobado, M. J. Herrero, and D. Temes. Phys. Rev., D65:075023, 2002. 

D. J. Gross and F. Wilczek. Phys. Rev. Lett, 30:1343-1346, 1973; H. D. 
Politzer. Phys. Rev. Lett, 30:1346-1349, 1973. 

R. Brock et al. Rev. Mod. Phys., 67:157-248, 1995; R. Brock et al. Hand- 
book of perturbative QCD; Version 1.1: September 1994, FERMILAB-PUB- 
94-316, 1994. 

J. D. Bjorken and E. A. Paschos. Phys. Rev., 185:1975-1982, 1969; R. P. 
Feynman. Phys. Rev. Lett, 23:1415-1417, 1969. 

E. D. Bloom et al. Phys. Rev. Lett, 23:930-934, 1969; M. Breidenbach 
et al. Phys. Rev. Lett, 23:935-939, 1969; J. I. Friedman and H. W. Kendall. 
Ann. Rev. Nucl. Part Sci., 22:203-254, 1972. 

J. C Collins, D. E. Soper, and G. Sterman. Adv. Ser. Direct. High Energy 
Phys., 5:1-91, 1988. 

G. Altarelli and G. Parisi. Nucl. Phys., B126:298, 1977. 

C. Meier. ^Z-Produktion in P-P-Kollisionen mit elektroschwachen 1- 
Schleifen-Korrekturen. Diplom thesis, Universitat Karlsruhe, July 2001. 

H. Plothow-Besch. PDFLIB - User's Manual Version 8.O4. CERN- 
ETT/TT 2000.04.17. 

W. Giele and M. Whalley. LHAPDF version 4 Users Guide, April 2002. 
http : //durpdg . dur . ac . uk/lhapdf /manual . htm. 

Z. Kunszt and F. Zwirner. Nucl. Phys., B385:3-75, 1992. 

J. F. Gunion, H. E. Haber, F. E. Paige, W.-K. Tung, and S. S. D. Wil- 
lenbrock. Nucl. Phys., B294:621, 1987; R. M. Barnett, H. E. Haber, and 

D. E. Soper. Nucl. Phys., B306:697, 1988; F. I. Olness and W.-K. Tung. 
Nucl. Phys., B308:813, 1988. 

V. D. Barger, R. J. N. Phillips, and D. P. Roy. Phys. Lett, B324:236-240, 
1994. 

C. S. Huang and S.-H. Zhu. Phys. Rev., D60:075012, 1999; L. G. Jin, C. S. 
Li, R. J. Oakes, and S. H. Zhu. Eur. Phys. J., C14:91-101, 2000; L. G. 
Jin, C. S. Li, R. J. Oakes, and S. H. Zhu. Phys. Rev., D62:053008, 2000. 



Bibliography 161 

[98 
[99 



[100 
[101 

[102 

[103 
[104 
[105 

[106 
[107 

[108 
[109 
[110 

[111 
[112 



[113 
[114 
[115 



K. A. Assamagan et al. The Higgs working group: Summary report 2003, 
hep-ph/0406152, 2004. 



J. G. Branson et al. High transverse momentum physics at the Large Hadron 
ColUder, hep-ph/0110021, 2001. 

V. Buscher and K. Jakobs. Int. J. Mod. Phys., A20:2523-2602, 2005. 

D. A. Dicus, J. L. Hewett, C. Kao, and T. G. Rizzo. Phys. Rev., D40:787, 
1989. 

A. A. Barrientos Bendezu and B. A. Kniehl. Phys. Rev., D59:015009, 1999; 
A. A. Barrientos Bendezu and B. A. Kniehl. Phys. Rev., D61:097701, 2000. 

S. Moretti and K. Odagiri. Phys. Rev., D59:055008, 1999. 

W. Hollik and S.-h. Zhu. Phys. Rev., D65:075015, 2002. 

Y.-S. Yang, C.-S. Li, L.-G. Jin, and S. H. Zhu. Phys. Rev., D62:095012, 
2000. 

O. Brein, W. Hollik, and S. Kanemura. Phys. Rev., D63:095001, 2001. 

R. N. Cahn and S. Dawson. Phys. Lett, B136:196, 1984; Erratum-ibid. 
Phys. Lett, B138:464, 1984. 

D. A. Dicus and S. S. D. Willenbrock. Phys. Rev., D32:1642, 1985. 

G. Altarelh, B. Mele, and F. Pitolh. Nucl. Phys., B287:205-224, 1987. 

H. M. Georgi, S. L. Glashow, M. E. Machacek, and D. V. Nanopoulos. 
Phys. Rev. Lett, 40:692, 1978. 

D. Graudenz, M. Spira, and P. M. Zerwas. Phys. Rev. Lett., 70:1372-1375, 
1993; M. Spira, A. Djouadi, D. Graudenz, and P. M. Zerwas. Nucl. Phys., 
B453:17-82, 1995; M. Spira. Fortsch. Phys., 46:203-284, 1998. 

S. Catani, D. de Florian, and M. Grazzini. JHEP, 05:025, 2001; R. V. Har- 
lander and W. B. Kilgore. Phys. Rev., D64:013015, 2001; R. V. Harlander 
and W. B. Kilgore. Phys. Rev. Lett, 88:201801, 2002; C. Anastasiou and 
K. Melnikov. Nucl. Phys., B646:220-256, 2002; V. Ravindran, J. Smith, 
and W. L. van Neerven. Nucl. Phys., B665:325-366, 2003. 

T. Han and S. Willenbrock. Phys. Lett, B273:167-172, 1991. 

T. Figy, C. Oleari, and D. Zeppenfeld. Phys. Rev., D68:073005, 2003. 

E. L. Berger and J. Campbell. Phys. Rev., D70:073011, 2004. 



162 Bibliography 

[116] D. Zeppenfeld, R. Kinnunen, A. Nikitenko, and E. Richter-Was. Phys. Rev., 
D62:013009, 2000; D. Zeppenfeld. eConf, C010630:P123, 2001. 

[117] ATLAS Collaboration. ATLAS Detector and Physics Performance - Tech- 
nical Design Report, CERN/LHCC 99-14, 5 1999. 

[118] A. Djouadi and M. Spira. Phys. Rev., D62:014004, 2000. 

[119] A. D. Martin, R. G. Roberts, W. J. Stirling, and R. S. Thorne. Phys. Lett, 
B531:216-224, 2002. 

[120] A. D. Martin, R. G. Roberts, W. J. Stirling, and R. S. Thorne. Eur. Phys. 
J., C28:455-473, 2003. 

[121] R. Raitio and W. W. Wada. Phys. Rev., D19:941, 1979; A. S. Bagdasaryan, 
R. S. Egorian, S. G. Grigorian, and S. G. Matinyan. Sov. J. Nucl. Phys., 
46:315, 1987; J. N. Ng and R Zakarauskas. Phys. Rev., D29:876, 1984. 

[122] Z. Kunszt. Nucl. Phys., B247:339, 1984. 

[123] C. Balazs, H.-J. He, and C. R Yuan. Phys. Rev., D60: 114001, 1999; D. Di- 
cus, T. Stelzer, Z. Sullivan, and S. Willenbrock. Phys. Rev., D59:094016, 
1999. 

[124] S. Dittmaier, M. Kramer, and M. Spira. Phys. Rev., D70:074010, 2004; 
S. Dawson, C. B. Jackson, L. Reina, and D. Wackeroth. Phys. Rev., 
D69:074027, 2004; S. Dawson, C. B. Jackson, L. Reina, and D. Wackeroth. 
Higgs production in association with bottom quarks at hadron colliders, hep- 
ph/0508293, 2005. 

[125] W. Beenakker et al. Phys. Rev. Lett, 87:201805, 2001. 

[126] L. Reina and S. Dawson. Phys. Rev. Lett, 87:201804, 2001; S. Dawson, 
L. H. Orr, L. Reina, and D. Wackeroth. Phys. Rev., D67:071503, 2003. 

[127] W. Beenakker et al. Nucl. Phys., B653:151-203, 2003. 

[128] D. A. Dicus and S. Willenbrock. Phys. Rev., D39:751, 1989. 

[129] J. Campbell et al. Higgs boson production in association with bottom quarks, 
hep-ph/0405302, 2004; J. L. Diaz-Cruz, H.-J. He, T. Tait, and C. R Yuan. 
Phys. Rev. Lett, 80:4641-4644, 1998; C. Balazs, J. L. Diaz-Cruz, H. J. He, 
T. Tait, and C. R Yuan. Phys. Rev., D59:055016, 1999. 

[130] J. Dai, J. F. Gunion, and R. Vega. Phys. Lett, B345:29-35, 1995; J. Dai, 
J. F. Gunion, and R. Vega. Phys. Lett, B387:801-803, 1996; E. Richter- 
Was and D. Froidevaux. Z. Phys., C76:665-676, 1997. 



Bibliography 163 

[131] R. V. Harlander and W. B. Kilgore. Phys. Rev., D68:013001, 2003. 

[132] G. Gao, R. J. Oakes, and J. M. Yang. Phys. Rev., D71:095005, 2005. 

[133] W. J. Marciano and F. E. Paige. Phys. Rev. Lett, 66:2433-2435, 1991; 
J. F. Gunion. Phys. Lett, B261:510-517, 1991. 

[134] J. Goldstein et al. Phys. Rev. Lett, 86:1694-1697, 2001. 

[135] E. Richter-Was and M. Sapinski. Acta Phys. Polon., B30:1001-1040, 1999; 
V. Drollinger, T. Muller, and D. Denegri. Searching for Higgs bosons in 
association with top quark pairs in the HO -^ b anti-b decay mode, hep- 
ph/0111312, 2001. 

[136] F. Maltoni, D. L. Rainwater, and S. Willenbrock. Phys. Rev., D66:034022, 
2002; A. Belyaev and L. Reina. JHEP, 08:041, 2002. 

[137] W. Peng et al. NLO supersymmetric QCD corrections to t anti-t hO asso- 
ciated production at hadron colliders, hep-ph/0505086, 2005. 

[138] K. Desch and M. Battaglia. Determination of the Higgs profile: HFITTER, 
Prepared for 5th International Linear Collider Workshop (LCWS 2000), 
Fermilab, Batavia, Illinois, 24-28 Oct 2000, 2000. 

[139] W. Kilian, T. Plehn, P. Richardson, and E. Schmidt. Eur. Phys. J., 
039:229-243, 2005. 

[140] T. Plehn, M. Spira, and P. M. Zerwas. Nucl. Phys., B479:46-64, 1996; 
T. Plehn, M. Spira, and P. M. Zerwas. Nucl. Phys., B531:655, 1998. 

[141] D. A. Dicus, 0. Kao, and S. S. D. Willenbrock. Phys. Lett, B203:457, 
1988; E. W. N. Glover and J. J. van der Bij. Nucl. Phys., B309:282, 1988; 
A. Djouadi, W. Kilian, M. Muhlleitner, and P. M. Zerwas. Eur. Phys. J., 
010:45-49, 1999; S. Dawson, S. Dittmaier, and M. Spira. Phys. Rev., 
D58:115012, 1998. 

[142] U. Baur, T. Plehn, and D. L. Rainwater. Phys. Rev. Lett., 89:151801, 2002; 
U. Baur, T. Plehn, and D. L. Rainwater. Phys. Rev., D67:033003, 2003; 
A. Dahlhoff. Prospects to measure the Higgs boson properties in ATLAS, 
hep-ex/0505022, 2005; F. Gianotti et al. Eur. Phys. J., 039:293-333, 2005. 

[143] A. Djouadi, W. Kilian, M. Muhlleitner, and P. M. Zerwas. Eur. Phys. J., 
010:27-43, 1999. 

[144] 0. J. Gounaris, D. Schildknecht, and F. M. Renard. Phys. Lett, B83:191, 
1979; V. D. Barger, T. Han, and R. J. N. Phillips. Phys. Rev., D38:2766, 
1988; F. Boudjema and E. Ohopin. Z. Phys., 073:85-110, 1996; D. J. 



164 Bibliography 

Miller and S. Moretti. Eur. Phys. J., 013:459-470, 2000; M. Battaglia, 

E. Boos, and W.-M. Yao. eConf, C010630:E3016, 2001; C. Castanier, 
P. Gay, P. Lutz, and J. Orloff. Higgs self coupling measurement in e+ 
e- collisions at center- of -mass energy of 500-GeV, liep-ex/0101028, 2001; 
U. Baur, T. Plehn, and D. L. Rainwater. Phys. Rev., D68:033001, 2003. 

[145] U. Baur, T. Plehn, and D. L. Rainwater. Phys. Rev., D69:053004, 2004. 

[146] E. Accomando et al. Physics at the CLIC multi- Te V linear collider, liep- 
ph/0412251, 2004. 

[147] J. R. Ellis, M. K. Gaillard, and D. V. Nanopoulos. Nucl. Phys., B106:292, 
1976; M. A. Shifman, A. I. Vainshtein, M. B. Voloshin, and V. I. Zakharov. 
Sov. J. Nucl. Phys., 30:711-716, 1979; B. A. Kniehl and M. Spira. Z. Phys., 
C69:77-88, 1995. 

[148] F. Jegerlehner. Nucl. Phys. Proc. Suppl, 51C:131-141, 1996; 

F. Jegerlehner. Hadronic vacuum polarization effects in alpha(em)(M(Z)), 
hep-ph/0308117, 2003. 

[149] G. 't Hooft and M. J. G. Veltman. Nucl. Phys., B153:365-401, 1979. 

[150] W. Beenakker and A. Denner. Nucl. Phys., B338:349-370, 1990. 

[151] A. Denner, U. Nierste, and R. Scharf. Nucl. Phys., B367:637-656, 1991. 

[152] A. Denner and S. Dittmaier. Nucl. Phys., B658:175-202, 2003. 

[153] G. Passarino and M. J. G. Veltman. Nucl. Phys., B160:151, 1979. 

[154] T. Binoth, J. P. Guillet, G. Heinrich, E. Pilon, and C. Schubert. JHEP, 
10:015, 2005. 

[155] A. Denner and S. Dittmaier. Nucl. Phys., B734:62-115, 2006. 

[156] G. J. van Oldenborgh. Comput. Phys. Commun., 66:1-15, 1991. 

[157] S. Dittmaier. Private communication. 

[158] F. Bornemann and C. Lasser. Numerische Mathematik, WS 2004/05. Lec- 
ture and tutorials at TU Miinchen. 

[159] I. C. Society. IEEE Standard for Binary Floating-Point Arithmetic. Tech- 
nical report, IEEE Std 754, 1985. 

[160] P. Deuflhard and A. Hohmann. Numerische Mathematik: eine algorith- 
misch orientierte Einfiihrung. de Gruyter, 1993. 



Bibliography 165 



[161] J. H. Wilkinson. J. ACM, 8(3):281-330, 1961. 

[162] S. Heinemeyer, W. Hollik, and G. Weiglein. Comput. Phys. Commun., 
124:76-89, 2000. 

[163] T. Hahn. Comput. Phys. Commun., 168:78-95, 2005. 



Acknowledgments 



First of all, I would like to cordially thank my PhD advisor, Prof. Wolfgang HoUik, 
for his continuous support and the possibility to do my PhD thesis in his working 
group. Despite his many obligations as one of the directors of the Max-Planck- 
Institut fiir Physik, he has always found the time to answer my questions, discuss 
new results, and give many further valuable suggestions and advice. 

I would also like to thank Prof. Manuel Drees for initially taking over the 
formal supervision of my PhD. 

I am also indebted to Tilman Plehn for his support and for the many fruitful 
discussions about Higgs physics and collider phenomenology. The work on our 
joint paper has always been much fun. A big thanks as well to Heidi Rzehak for 
answering all my questions about renormalization in different schemes, for many 
valuable discussions and for suggesting to me to join such a great working group. 
Also Pietro Slavich has taught me many details about quantum field theory and in 
particular about the universal corrections to the bottom-quark Yukawa coupling. 
Thomas Hahn has always had an answer for even the trickiest challenges in pro- 
gramming, regardless whether the problems occurred in Mathematica, Fortran, 
UNIX-shell or ETeX. 

Stefan Dittmaier has been of great help during the implementation of the 
five-point loop integrals and has also kindly provided his code for a numeri- 
cal comparison. The development of HadCalc was greatly influenced by Oliver 
Brein, who has told me all details about calculating hadronic cross sections and 
whose private computer code for the calculation of hadronic cross sections has 
triggered my work on HadCalc. I would also like to thank Martin Fiirst and 
Christian Schappacher for their feedback on the program, the bug reports and 
the helpful suggestions to make the program more useful and user-friendly. Also 
Cailin Farrell, Axel Bredenstein and Thomas Fritzsche have invested much time 
in answering all my questions about their respective work, from effective field 
theories over Monte Carlo generators to the renormalization of the MSSM. The 
great atmosphere at the institute and all other colleagues, which I have not men- 
tioned so far and which would go beyond the scope of these acknowledgments, 
have contributed significantly to the success of this thesis as well. 

Performing numerical calculations would not have been possible without the 
excellent maintenance of the computer system by Thomas Hahn, Peter Breiten- 
lohner and all other system administrators of the MPI and the Rechenzentrum 
Garching der MPG, where parts of the calculations were carried out. 

Many thanks also to Heidi Rzehak, Ulrich Meier, Michael Pliimacher, Christoph 
ReiBer and Axel Bredenstein for proof-reading several chapters of my thesis each 
and contributing many useful corrections and suggestions, and especially to Arne 
Weber, who has not only carefully read large parts of this dissertation and given 
helpful remarks, but also greatly improved its linguistic aspect. 

Finally, I would like to thank my parents for their support. 



